a Эмєz_evalf..zNot Enough Initial ConditionsrщNщ cSsg|]}t|dГСqS)rrrrrr r)r)Zannihilator┌order┌parent┌base┌ get_field┌_euler┌_rk4Z listofpoly┌append┌new┌rep┌range┌y0┌len┌ TypeError┌x0┌ enumerater)┌funcZpointsZderivatives┌method┌ann┌R┌K┌meth┌j┌redrr┌solr rrr ┌_evalf s* r*c Cs╕t|Гаtjб}t|Гаtjб}ddД|DГ}||}|ddЕ} d} t|ГD].}| tt|||ddНГаtjб||7} qN| а| бg}t|ГD]}|а|||| |бqФ|S)zs Euler's method for numerical integration. From x0 to x1 with initial values given at x0 as vector y0. cSsg|]}t|ГаtjбСqSrйr┌ _to_mpmathr┌precrrrr r6rz_euler..rNrTйZmpmйrr,rr-rrr) r(r┌x1rr┌A┌B┌y_0┌h┌f_0┌f_0_nr r)rrr r.s, rc sHt|Гаtjб}t|Гаtjб}ddД|DГЙ||Йd}d}d} d} ИddЕЙt|ГD].}|tt|||ddНГаtjбИ|7}qZИа|бЗЗЗfddДtd|ГDГЙt|ГD]F}|tt|||Иd ddНГаtjбИ|И|Иd 7}q╕Иа|бЗЗЗfd dДtd|ГDГЙt|ГD]H}| tt|||Иd ddНГаtjбИ|И|Иd 7} Рq.Иа| бЗЗЗfddДtd|ГDГ}t|ГD]@}| tt|||ИddНГаtjбИ|И|И7} Рqж|а| бg} t|ГD]D}| аИ|ИИ|d И|d И|||dбРq■| S) z1 Runge-Kutta 4th order numerical method. cSsg|]}t|ГаtjбСqSrr+rrrr rMrz_rk4..rrNTr.cs$g|]}И|И|ИdСqSйщrr)r5r4r3rr rZrr8cs$g|]}И|И|ИdСqSr7rr)┌f_1r4r3rr r_rcs g|]}И|И|ИСqSrrr)┌f_2r4r3rr rdrщr/)r(rr0rrr1r2r6Zf_1_nZf_2_nZf_3_nr Zf_3r)r)r5r9r:r4r3r rFs8, D F > BrN)Fr) ┌__doc__Zsympy.core.sympifyrZsympy.holonomic.holonomicrZmpmathrr*rrrrrr ┌s %