g| ddlmZmZmZmZddlmZmZmZm Z ddl m Z ddl m Z ddlmZmZddlmZmZddlmZmZddlmZdd lmZdd lmZmZmZmZdd l m!Z!dd l"m#Z#dd l$m%Z%ddl&m'Z'ddl(m)Z)m*Z*m+Z+ddl,m-Z-ddl.m/Z/m0Z0m1Z1m2Z2m3Z3ddl4m5Z5ddl6m7Z7ddl8m9Z9ddl:m;Z;mZ>ddl?m@Z@ed\ ZAZBZCZDZEZFZGZHZIZJZKedd\ZLZMZNZOedZPedZQedZRdZSdZTd ZUe=d!ZVd"ZWd#ZXd$ZYd%ZZd&Z[d'Z\e;d(Z]d)Z^d*Z_d+Z`d,Zad-Zbd.Zcd/Zdd0Zed1Zfd2Zgd3Zhd4Zid5Zje=d6Zkd7Zld8Zmd9Znd:Zod;Zpd<Zqd=Zrd>Zsd?Ztd@Zue>dAdBZvdCS)D) DerivativeFunctionSubsdiff)EIRationalpi)Eq)S)Symbolsymbols)imre)explog)acosh)sqrt)atan2cossintan)Integral)Poly)O)collect) classify_odehomogeneous_orderdsolve) checkodesol)classify_sysodeconstant_renumber constantsimpget_numbered_constants solve_ics) _undetermined_coefficients_match)LinearCoefficients) ode_order)XFAILraisesslowSKIP) filldedentzC0:11zu,x:zTrealfghc~tt5tdddddS#1swxYwYdSN)r* ValueErrorr$\/var/www/html/ai-engine/env/lib/python3.11/site-packages/sympy/solvers/ode/tests/test_ode.pytest_get_numbered_constantsr90s   %%t$$$%%%%%%%%%%%%%%%%%%s 266c L ttt}t|d}|d}|jjdd}|d}|jjdjdd}idtttttdtzdttttd ttttd ttttd ttttd ttttd ttttdttttdttttdttttdtttttdtzdttttdtttdttd|z|tttdttttdttttdtttttd|z |tttz fzdttttttd|ttfttdtzttttttttd d dd}||ksJt|dttttksJdS)Nallhintseparable_Integralr11st_homogeneous_coeff_subs_dep_div_indep_IntegralBernoulli_Integral1st_homogeneous_coeff_best Bernoulli nth_algebraicnth_linear_euler_eq_homogeneous%nth_linear_constant_coeff_homogeneous separable(1st_homogeneous_coeff_subs_indep_div_depnth_algebraic_Integral 1st_linear1st_linear_Integral 1st_exact1st_exact_Integral lie_group(1st_homogeneous_coeff_subs_dep_div_indep1st_power_series)r>11st_homogeneous_coeff_subs_indep_div_dep_Integralbest best_hintdefaultorderrS) r0xrrlhsargsrhsr C1rrr)eqoutputsol1_y_u1expecteds r8test_dsolve_all_hintrb5sI 11B BU # # #F & 'D q ! B E FD (-   "1 %C$b1rHQNN/B&C&C$b1rllR!b\\ AaDD" *2addB<<  0AaDD"  R!b\\  3BqttRLL !"QqTT2,, b1rll r!A$$Xa^^(;<< R!b\\ bhq!nnx2&FAaDD!Q!QSUVV R!b\\ 3BqttRLL ) rMr@rRrKrArlrSrTrUrmrErVr>rVrSrUrmrTrLrGrBrOrHrJrKrMr@rRr>c&tdS)N notarealhintr<rr\sr8z%test_dsolve_options..svb~>>>r7c&tdS)N Liouviller<rprqsr8rrz%test_dsolve_options..svb{;;;r7r?)rWr0rrsortedkeysr(r r[hasrr*r5rr)abcrw Integral_keysr\s @r8test_dsolve_optionsr}as 1Q4499Q<<!A$$ BrAr...Ar'''A + + +D'''M !&&((  t # # # # W:2qtt,, , , , , V91Q44A & & & & "6 " " "b1r!tnn 4 4 4 4 Y<< ' ' ' ' [>\ ) ) ) )~!!(+++++~!!(+++++-.228<<<<<;<@@JJJJJ;<@@JJJJJ""8,,,,, " # ' ' 1 1111 ! " & &x 0 0000 @ A E Eh O OOOO @ A E Eh O OOOO ! " & &x 0 0000 !&&((  t # # # # W:2qtt,, , , , , V91Q44A & & & & "6E 2 2 2b1r!tnn D D D D Y<< ' ' ' ' [>\ ) ) ) ) [>Q{^ + + + + 7 8 45 6 6 6 6 7 8 45 6 6 6 6~!!(+++++~!!(+++++-.228<<<<<;<@@JJJJJ;<@@JJJJJ""8,,,,, " # ' ' 1 1111 ! " & &x 0 0000 @ A E Eh O OOOO @ A E Eh O OOOO ! " & &x 0 0000 !&&((  } , , , , :>>>>??? :;;;;<<< !A$$))A,,1Q447* 7 7 7 ?qttyy||a!ai'f555 6 6 6 6 !A$$11%A A6:AaDD, . . . 1Q44"x#a&&++a..1!4HQNN##!$%&(((),hq!nn_)=)=> ? ? ? ? ? ? ? ?r7c ttttdttdksJtttttdksJtt tttdttdksJttttdzttdksJtt tttttzttt}ttttttzttttzztttzz tt}ttttttz ttttz ztttz z tt}|dksJ|dksJ|dksJtdtzttztttzd tzttdzzzt tz ttd ksJd tdtttzttdzztttz z d td zzzttvsJt tdtd}tttt|ttz|tzttzzz dttz|ttz|tzttzzz zttttz|ttz|tzttzzz ztzttdksJd}tttttztttttz|ksJtttttztttttzttd|ksJtt dtd zztttzdttdksJtt dttd zztttzdttdksJtt tttttttzz dttdksJttt ttdttdtsJttt tttdttd gdksJtdtzttztttzd tzttdzzzt tz ttdd}t| gdksJttttt tttz z tt}|dksJdS)Nru)rDrFrErt2nd_power_series_ordinaryrILiouville_Integral)rDrIr)rDrGrLrJrCrBrHrOrQrNrFrErIr>rMrKrArRr@)rmrDrGrLrJrCrBrHrOrQrNrErIr>rMrKrArRr@) rLrJrCrkrQrN3nth_linear_constant_coeff_undetermined_coefficients1nth_linear_constant_coeff_variation_of_parametersrMrKrArl:nth_linear_constant_coeff_variation_of_parameters_Integral) rmrJrCrQrNrrrKrArr?)rmrLrCrkrNrMrArlRiccati_special_minus2ctttttttztttSr4)rrWr0yrr6r7r8rrz#test_classify_ode..sI|A!Q Q0D0D 1 1 - aGG  r7k) rmrGrLrJrCrQrNr>rMrKrA)rmrDrGrLrJrCrBrHrOrQrNrrMrKrArRr@rCnth_linear_euler_eq_nonhomogeneous_variation_of_parameters_IntegralT)prep)rmrDrGrLrJrCrQrNrErIr>rMrKrA) rmrDrGrLrJrCrQrNrIr>rMrKrA)rQrNdict)rUrFrVrCrr=)rCrArUrV ordered_hints)rGrLrQrNr>rM)rr0rWrr rr*r5r z isinstancerrvrw)ryrzr{ranss r8test_classify_odersq ! !Q1 . .   !add # #'R R R R R 1Q4499Q<<++QqTT 2 27= = = = =& ! ! a1 . .3> > > > >( R! ! qtt+Q//166AQqTTYYq\\!A$$&1add2QqttV;QqTTBBAQqTTYYq\\!A$$&1add2QqttV;QqTTBBA F F F F F G G G G G G G G G G  !AaDD11Q!a/#a&&8!A$$   N N N N N $Qqttyy||^addAg-!Q61r7BAaDDII J J J J : s A ! ! a!fqs1Q44x&781QqTT61QqTT6 !AaDDD<addiillNAaddFQqS1X$5679:;<=aDD B B 5 5 5 5 5  OC$ QqTTAXq))D1qMM9 : :c A A A A QqTTAXq))D1qMM91Q44! # # #&) * * * * 1QT6!A$$))A,,.22AaDD 9 9 6 6 6 6 6 1QqTT1W9QqTTYYq\\1155qtt < <      4!a==1Q4472A66! = =) * * * * l2addA;;!4@@@$ G GGGG ,r!A$$))A,,22AaDDtDDDIIKK L LEEE F F F FQqS1Xaddiill*a!eQqTT1W_.N|B!#>>>r7rcLtttSrrrsr8rrz'test_classify_ode_ics..Xrr7cLtttSrrrsr8rrz'test_classify_ode_ics..\rr7cLtttSrrrsr8rrz'test_classify_ode_ics..`rr7cLtttSrrrsr8rrz'test_classify_ode_ics..drr7rdcLtttSrrrsr8rrz'test_classify_ode_ics..urr7cLtttSrrrsr8rrz'test_classify_ode_ics..yrr7cLtttSrrrsr8rrz'test_classify_ode_ics..}rr7cLtttSrrrsr8rrz'test_classify_ode_ics..rr7cLtttSrrrsr8rrz'test_classify_ode_ics..rr7cLtttSrrrsr8rrz'test_classify_ode_ics..rr7cLtttSrrrsr8rrz'test_classify_ode_ics..rr7cLtttSrrrsr8rrz'test_classify_ode_ics..rr7cLtttSrrrsr8rrz'test_classify_ode_ics..rr7cLtttSrrrsr8rrz'test_classify_ode_ics..rr7cLtttSrrrsr8rrz'test_classify_ode_ics..rr7) r0rWrr*r5r1rsubsrrrrs@@r8test_classify_ode_icsrHs 11a1Q44 B a&C :>>>>>??? Q44)C :>>>>>??? Q44)C :>>>>>??? Q77A,C :>>>>>??? Q441,C :>>>>>??? Q441,CQqTTs#### Q44)CQqTTs#### Q4499Q<<  Q " "A &C :>>>>>??? Q4499Q<<  Q " "A &C :>>>>>??? Q4499Q<<  Q " "A &C :>>>>>??? Q77<<??  1 % %q )C :>>>>>??? addAq ! ! & &q! , ,a 0C :>>>>>??? Q4499Q<<  Q " "AaDD )C :>>>>>??? Q4499Q<<  Q " "AaDDIIaLL$5$5a$;$; >>>>??? Q4499Q<<  Q " "A &C :>>>>>??? Q77<<??  1 % %q )C :>>>>>??? addAq ! ! & &q! , ,a 0C :>>>>>??? Q4499Q<<  Q " "AaDD )C :>>>>>??? Q4499Q<<  Q " "AaDD )CQqTTs#### Q4499Q<<  Q " "A &CQqTTs######r7chtdd\}}}}tdd\ }}}}}} } } } tdt\} }}}}}tdt\}}}}}}}tdt\}}}td }t|||}t|||}t|t |||z| ||||zt||||| ||||zf} d d ||d fd d ||d fd d ||d fd d ||d fd d ||d fd d ||d fd d ||d fd d ||d fd id ||||gd| |||| t |||zzt |||z| |||| |||zt |||zg||d ||d id}!t | |!ksJt|||d z||||z zt|||||z f}"d d ||d fd d ||d fd d ||d fd d ||d fd d ||d fd||z d ||d fd d ||d fd||z d ||d fd id||||gd||d z t |||z||||z z || ||z t |||zg||d ||d id}#t |"|#ksJt|||||||z||||||zt|||||||z||||||zf}$||||gdd||| |||z||||||zt |||z||| |||z||||||zt |||zgd ||d fd d ||d fd d ||d fd d ||d fd d ||d fd d ||d fd d ||d fd d ||d fd i||d ||d id d}%t |$|%ksJt|||||dzzt|||dzf}&d d ||d f||dz d ||d fd d ||d fd d ||d fd d ||d fd d ||d fd d ||d fd d ||d fd id||||gd|| ||dzzt |||z||dz t |||zg||d ||d id}'t |&|'ksJt|||||zt|d zzt|||d zt|d zzf}(d d ||d f|| t|d zzd ||d fd d ||d fd d ||d fd d ||d fd d ||d fd d ||d f|| t|d zzd ||d fd id||||gd|| ||zt|d zzt |||z||d z t|d zzt |||zg||d ||d id})t |(|)ksJdS)Nz k, l, m, nT)Integerz"k1, k2, k3, l1, l2, l3, m1, m2, m3zP, Q, R, p, q, rclszP1, P2, P3, Q1, Q2, R1, R2zx, y, ztrurr?type2F)no_of_equation func_coefftype_of_equationfunc is_linearr\rVrPtype3type4)rrrr\rrVrrrtype1)rrrr rrr!r)*rlmnk1k2k3l1l2l3m1m2m3PQRpqrP1P2P3Q1Q2R1R2rWrrrx1y1eq6sol6eq7sol7eq8sol8eq11sol11eq13sol13s* r8test_classify_sysoders t444JAq!Q)01U_c)d)d)d&BBBBB1x@@@Aq!Q1!()E8!T!T!TBBBBiX...GAq! A aadd1BT!!A$$q\\ b#a!f++aa!QQqTTll* + +R11QQqTT7711QQqTT!!A$$<<3G-H-H ICAAaDD! a!QQqTT1qSTVWVWXYVZVZ\]R^`adeghghijgkgkmncoqr!aL!a1q\1q!!A$$lA11Q44|Q0H]dQqTT11Q44Luaa!aaddmm^C!!A$$KK-G*UVUVWXUYUY[\J]J]-]`a`abcbcdebfbfAaDDa a ` !AAaDD''`qqttQ''`(-)45AaDD!QQqTT13E G GD 3  4 ' ' ' ' b!!A$$'!!A$$qqtt)# $ $bQQqTT!!A$$Y&7&7 8CAAaDD! a!QQqTT1qSTVWVWXYVZVZ\]R^`adeghghijgkgkmncoqr!aL"QQqTT'AqqttQ<Q!aL"QQqTT'AqqttQ+>>AAaDD ""QQqTT((@RSTSTUVUVWXUYUY[\[\]^[_[_abScSc@cqqttQA*+QQqTT1q1aaddA,AqqQRttUVaQqTT1WHz!!A$$""">$/0qttQ!a.@ B BE 4 E ) ) ) ) r!!A$$qqtt)CFFAI% & &2aaddAgc!ffai.?(@(@ AD !QQqTT1!uSVVQYQRTUTUVWTXTXZ[P\^_bcefefgheieiklamop!aL!a1q\1q!!A$$lA11Q44|aaddU3q66ST9_WXZ[Z[\]Z^Z^`aVbde1g!!A$$!ESTSTUVSWSWRWXYXYZ[X\X\R\]`ab]c]cef]fRfqqttQS!A$$'#a&&!)+j1q.A.AARCNOaPQddTUWXWXYZW[W[]^M_ a aE 4 E ) ) ) ) ) )r7c d,ttttdttz z tttddit ttt dtzdzksJttttttz tttddit ttt tksJttttttz ttttttddit ttt tksJtttttttztttddttttddit ttttttzksJttttttz ttztttttz ttz gttttgtddtddit ttt tttzt ttt tttzgksJtdzttdzztz t}t|tttddit ttt tdz  tz t ttt tdz tz gksJt|ttttttddit ttt ttj z  tz t ttt ttj z tz gksJtttttttzttdzz tttzz }t|tttddiddt ttttzttdzdz ztddksJt|tttddidd t ttttzttdzdz ztddksJtt tttt tzgttgtgtdditdiksJtt tttttzt ttzzgttgtt gtddtt"dz ditdt diksJtt tttttzt ttzzgttgtt gtddttttdditdt diksJtt tttttzt ttzzgttgtt gtddit diksJt%tttttttdtz ztddttttddit ttdt tttdzd z tdz z hksJt%tttttttdtz ztddit ttdt ttt tztdzd z zhksJt'd \}}}t tt||zt |tzz| |z|t |tzz| |zz dz zz }tt ttt|ttzdtt|z z ztttd|i|ksJt%tttttttdtz zttttddtddit ttdt tttdzd z hksJt%tttttttdtz ztddttttddit ttdt tttdzd z hksJt)t*d t)t*d t'd\}}}t tttt tzzt,tdzzzt.tdzzz|tdzzd|zz zg}ttg} tt t,t.g} tddttttddt||ddt||ddi} tddttttddt1tttdt|dt1tttdt|di} t|| | | } t|| | | }| |cxkr8tdt dt,|dz|zd|zz t.| |zd |zz iksnJtdtdi}ttttttz tt|t tttdt tzksJttttdttttdtdtdi}tttttttztt|t tttdttzttttdttzzksJdS)Nr?rurdrrLF)rer=rirTzK r f0c Ztttttt tzt ttzzgttgtt gtddttdiS)Nrr?) r%r r0rWr[rC2rr r6r7r8rrz test_solve_ics..sny"QqTT2c!ff9r#a&&y3H*I*I)JQqTTFUWY[T\_`ab_c_cefhijlhmhmop^qrrr7c ^ttttt t tdzz gttgt gttttddiS)Nrur)r%r r0rWrr[rrr6r7r8rrz test_solve_ics..spy"QqTT4QT ??*C*C)DqttfrdUVWXUYUYU^U^_`UaUaUfUfghjkUlUlnoTpqqr7zEI q L)rr0rWrr rrrrrr1r Halfr r%r[rr setrr*r5C3C4r)r\Krf0solEIrLsolsfuncs constantsics1ics2solved_constants1solved_constants2res r8test_solve_icsrs !A$$))A,,1Q44'1AaDD!9 = = = 1Q44a!eai!! " " " " !A$$))A,,1%qtt!A$$ ; ; ;r!A$$A?O?O O O O O !A$$))A,,1%qtt!A$$))A,,2C2CAq2I2I11M N N NRTUVWXUYUY[^_`[a[aRbRb b b b b !A$$))Aq//AaDD(!A$$QqTT1 ! ! !Q5$ % % %(*1Q44Q#a&&(A(A B B B B 1Q4499Q<UW_`acdWeWe;f;f g g g g "add qTT!H; 7 7 7:>!1Q44QS* U U UZ] ] ] ] ] vaddiillAaDDIIaOOA$56QqTTYYq\\=N=NqQR=S=STUWXYZW[W[\]<^___ ` ` AaDD!b1qAvz**+ , , , , vaddiillAaDDIIaOOA$56QqTT!QqTTYYq\\EVEVWXYZE[E[\]<^___ ` ` AaDD!b1qAvz**+ , , , , :rrsss :qqrrrx  HB1 qttR"Q$YAqD(2ad72Qq!tVRU^C D D ED qTTFERR I aDD! ! ! !Q ! !Q ! !Q D aDD! ! ! !Q QqTTYYq!__a##Q QqTTYYq!__a##Q (D "$y$??!$y$??  1     A A AqDFAbDM QBqD!B$K 6       Q441,C !A$$))A,,1%qtt 5 5 5AaDD!A$$s1vv+9N9N N N N N Q4499Q<<  Q " "AaDDIIaLL$5$5a$;$;QqTT1Q44 HC !A$$))Aq//AaDD(!A$$C 8 8 8 1Q441c!ffqttyy||00A66s1vv==>> ? ? ? ? ? ?r7c  td}td}td}td|zt||z||dksJt|t |||zd|z||zzt ||z z ||dksJt|dz||||z|t |||zz||z ||dksJtt |t||z||||dksJtt |t |t||z||z|||dksJtt ||||||dksJtt ||||t |||z||dksJtt ||||t |||z||dksJtt |t |t||z||z|||dksJtt |||z|||dksJt|t t |||dzz||z||dksJtt |t |t||z||z|||dksJtt ||||||dksJtt |t||z||||dksJtt ||||t |||z||dksJtt |t ||||z|||dksJt|t t |t |||dzz||z||dksJdS)Nr0r1rWrrr?ru)rr r(rrrr)r0r1rWs r8test_ode_orderr=s( A As A QqSQQqTT]AAaDD ) )Q . . . . QtAAaDD!}}_qs11Q44x/#a&&(:AAaDD A AQ F F F F QT!!A$$))Aq//)Ad11Q44mmO;aaddBAAaDD I IQ N N N N T!C!II+q!,,aadd 3 3q 8 8 8 8 T!D3qqtt99a333Q771 > >! C C C C T!!A$$1%%qqtt , , 1 1 1 1 T!!A$$1%%d11Q44mm3QQqTT : :a ? ? ? ? T!!A$$1%%d11Q44mm3QQqTT : :a ? ? ? ? T!D3qqtt99a333Q771 > >! C C C C Z!!A$$**AAaDD 1 1Q 6 6 6 6 Qs:a!aiA667771 > >! C C C C Z*Qs11Q44yy[!Q"?"? ?CCQQqTT J Ja O O O O Z!a++QQqTT 2 2a 7 7 7 7 Z#aadd)) Q22AAaDD 9 9Q > > > > Z!a++JqqttQ,?,??1 F F! K K K K Z*QQqTT1a"8"8 8!<.us0155r7)rrrrWrrrrrr r r0rr rrr*r5)rys r8test_homogeneous_orderrWs S1XXAaC0!Q 7 71 < < < < QTCFF3q66M11a 8 8 @ @ @ QUQs1Q3xxZ/A 6 6! ; ; ; ; acDA1$5$551c!ffs1vvo8NN Ax1~~  d1ggqj (*+,a 1 14 > > >s A Qs1Q3xxZ!CFF(2Aq 9 9Q > > > > QqTTYYq\\1a 0 0 8 8 8 addiill]Q.1 5 5 = = = QqTT1a ( ( 0 0 0 Q1a4[!Q / / 7 7 7 QUA & &" , , , , QT1 % % - - - :5566666r7ctdtz}t|t|}}|t r|t rJdSdS)Nr?)rrWrrrxr)rirs r8)test_noncircularized_real_imaginary_partsrxsY QqS A a55"Q%%qAe .u.......r7cdttdz z}|t|ttdz ksJdS)Nr?rur)rrWr)rs r8$test_collect_respecting_exponentialsrsB c!A#hh,C '3AaC)) ) ) ) ) ) )r7c+ttttddiksJttdtzt dztdt dtzt dztdtzt dzhdksJtttt tztddiksJt ttt tztdzt tztdzttztttztth}ttttdztzdzztd|dksJttttdzztttzzttztd|dksJtt dtzttztdztzdzztdt dtzttztdzt dtzzttzt tt dtzztdzt tzt dtzztt tzt dtzztt dtzzttzhdksJtdttz tddiksJttttddiksJtdtztdztzdzztddtztdtzztdzdtzzhdksJtttztddiksJtt tt dtzdzztdt ddtzzhdksJttttdztzdzztdtt tztttztdzt tztdzttzt ttthdksJttttttzztddiksJtttttdtzzztddiksJtttttztddiksJttdzttzt tztttzztztdtdzt tzt tztt ttj t tttztttt tzttztt tztt tzt tztttzt tt tztdzt tzttzh dksJtd tzttdz zttt tdz ztttdz zt tdz ttdz hdd ksJtdtztztddtztdtzzhdksJtdtzt dtzztddtzt dtzzhdksJtt t tz tddiksJttd tdtj hdksJtd t tztdt thdksJtt ttztdt ttzhdksJttttdt ttthdksJtt ttdt ttthdksJtd d t tzzdttzztdtj t tttt thdksJttdztdtj ttdzhdksJtdtzt tzt t ztdtt tzt tt t hdksJtdt dtzzttztdt dtzttzt tt dtzzhdksJttttz tdtj tt ttthdksJttdzdtzztdtj ttdzhdksJtd tzttztdtt tztttzt ttthdksJtttdtzztdtt dtzzttdtzzt dtztdtzhdksJttdzt t ztdtt t ztdzt t zt t hdksJtdt t ztdzt t zz tdtt t ztdzt t zt t hdksJtt dtztdzztdtj ttdzt dtzhdksJttt t ztdtt t zt t hdksJttt dtzztdtj tt dtzhdksJtttt t ztdt tttt t hdksJtt ttdt thdksJttj t dtzdz z tdtj t dtztdtzhdksJtt dtztj t dtzdz zztdt dtztdtzzt dtzt dtzzt dtzhdksJtdtzttzt tztdtj tt ttthdksJtt tdz t dtzdz z tdt tt dtztttdtzhdksJtt tdztddiksJtdtdzztddiksJdS)NtestFrurT)rtrialsetr?rr)rr r r)r&r1rWrrrrrrrr Onerr)ss r8$test_undetermined_coefficients_matchr s +AaDD! 4 4 G G G G +C!d1gg ,>,> B B 1tAww  QqS477]!3!3 4 6 6 6 6 6 6 ,CFF3q66M1 = =      Q3q661a4A;1SVV Qs1vvXs1vvFA +CFFAqD1HqL,A1 E E1%% & & & & + Aq!t c!ffQhQ' , ,8 > > > +AqD! 4 4 G G G G +CFF3qsQw<<,? C CCAaCLL>22 3 3 3 3 +CFFAqD1HqL,A1 E EAc!ffHaAh1SVV  1SVV SVVSVV$% & & & & & & ,CFFAAJ,? C C      +CFFAAaCL,A1 E E      +CFF3q66M1 = =      + 1SVV CFFQs1vvX%)1  1a4A;s1vv#5q#a&&!% As1vv s1vvqQxA#a&&!CFF(3q66/ #a&&#a&&Q-Ac!ffSVV!3#5 6 6 6 6 6 6 ,AaCAE NA > >s1q5zz\1SQZZ<QUSQZZHCC     ,AqDFA 6 6AqD!AqD&>22 3 3 3 3 +AqDQqSM1 = =AqDQqSM?33 4 4 4 4 +CGGAIq 9 9      ,AaDD! 4 4AE7++ , , , , +Bs1vvIq 9 9CFF8,, - - - - +C!HHa 8 8C!HH:.. / / / / +CFFA 6 6CFFCFF#344 5 5 5 5 +CFFA 6 6CFFCFF#344 5 5 5 5 +A#a&&L1SVV8,CQ G GAE3q663q663q66#BCC D D D D +AqD! 4 4AE1ad#344 5 5 5 5 +AaCAJaR,@! D DAc!ffHc!ffc1"gg#>?? @ @ @ @ +Ac!A#hhJs1vv,=q A AC!HHSVVOSVVC!HH_#EFF G G G G +AAJ : :AE1c!ffc!ff#=>> ? ? ? ? +AqD1Q3J : :AE1ad#344 5 5 5 5 +AaCAJ : :Ac!ffHaAhAA#GHH I I I I +Ac!A#hhJ : : s1Q3xxZ3qs88S1XXs1Q3xx 8 : : : : : : ,AqDaRL! < <Ac1"ggIq!tCGG|S!WW#EFF G G G G +Ac1"ggI1S!WW ,Da H HAc1"ggIq!tCGG|S!WW#EFF G G G G +C1II1,?? @ @ @ @ +Ac1"ggIq 9 9Ac1"ggIsA2ww#788 9 9 9 9 +AAaCL! < <AE1c!A#hh#788 9 9 9 9 +CFFS!WW, > > > +CFFA 6 6CFF8,, - - - - ,AFS1XXaZ,? C CAE3qs88S1XX#>?? @ @ @ @ , AaC!&3qs88A:%&  3qs88C!HH#4c!A#hhs1Q3xx6G AaC#       ,AaC#a&&L3q66,A1 E EAE1c!ffc!ff#=>> ? ? ? ? ,CFF1Hs1Q3xxz,A1 E ECFFC!HHc!ffc!A#hh#GHH I I I I +C1IIq 9 9fe_ L L L L +A1Iq 9 9fe_ L L L L L Lr7cddlm}t|ttt ttzt tzzztt ttzt tzdz}t |t tdksJtdzt tdzzt ttzdtzt tzz }t |t tdksJdS)Nr)Aru)rmrLrJrCrkrQrNrrrMrKrArlr) rmrLrBrHrOrQrNrMrRr@) sympy.abcrrWrr0r)rr\s r8test_issue_4785_22462rs  QD1qMM!AaDD( ) )D1qMM 9AaDD @1 DB AaDD ! !&F F F F F Q$1q.!A$$))A,, &1QqTT 1B AaDD ! !&= = = = = = =r7cttdtttt tt ttt zz ttdddddksJttdtttttt dddddksJdS)Nctttttttttzz ttSr4)rr0rWrrr6r7r8rrz!test_issue_4825..s9va1ggll1oo!Aq'' &A1Q44HHr7Trrr6)rVrUrcttttttSr4)rr0rWrrr6r7r8rrz!test_issue_4825..s&vaddiillAaDD99r7)r*r5rr0rWrrr6r7r8test_issue_4825rs :HHIII !Q Q!AaGG)3QqTT E E Er:: ; ; ; ; :99::: ! ! add 6 6 6r:: ; ; ; ; ; ;r7cddlm}tttzt t zzttt zt tzzcxkr#ttzt t zzksnJtt tzztt zz}|tt tgdD]1\}}}t||tzz|t zz|ksJ2dS)Nr) variationsr)sympy.utilities.iterablesrr"r[rWrrr)reryrzr{s r8'test_constant_renumber_order_issue_5308rs444444 RTBqD[ ) )"Q$A+&&     1r!t        BF R!VA:r2rlA..991a AEAE!233q8888899r7c$td\}}}}||z||z||zzg}t|d||zksJt|d|gt|zksJt|d|gtgt|zksJt|||gt|zt|zt|zzgksJt|||gtdt|zt|zt |zzgksJdS)Nze1:3 x yr) variables)r newconstantszC3:5)rr"r[rrr)e1e2rWrexprss r8test_constant_renumberr #s,:&&LBAq T2a4"Q$; E U1X & &"Q$ . . . . U1X! 5 5 5A = = = = U1X!B4 H H HBqD P P P P Uq!f 5 5 5"Q$1r!t 9L L L L L Uq!f76?? S S SXZ[\X\^`ab^beghiei^iWj j j j j j jr7c tdd}td}td}t|||d|dz||zz ||}t d|jDdksJt tttztttzzttztthtttzttzksJt tttztttzztttzzttthtttztttzzksJt tttzthtttzksJt ttztztthttzksJt ttztttdd fzthttzksJdS) NrTr.rwrcFg|]}|jd|S)C)name startswith).0r s r8 z#test_issue_5770..3s-FFFaqv/@/@/E/EFFFFr7r?ru)r rrrlen free_symbolsr#r[rrWrrrrrr)rrr"rs r8test_issue_5770r+.ssAs A A 11a1a4!9,aadd 3 3C FF3+FFF G G1 L L L L CFFRAY.A6R A A 3q66 #a&&     3q66 Bs1vvI-3q66 9BB< H H 3q66 Bs1vvI     BF bT * *bQi 7 7 7 7 B RG , ,Q 6 6 6 6 B!aAY!7!77" > >"q& H H H H H Hr7cttt ttztJttttz ttJdSr4)rrrWrrr6r7r8test_issue_5112_5430r-=sS c!ffWuQxx/ 3 3 ; ; ; QQZA . . 6 6 6 6 6r7cVtdttfddS)Nr0ctttdztdS)Nrufdsjf)rrWrr0sr8rrz!test_issue_5095..Ds1vaaddiillAoqqttWEEr7)rr*r5r1s@r8test_issue_5095r2Bs. A :EEEEFFFFFr7c td}tt|tz}tdtzd|tzz }t t|tzt ddz}t dtzd|tzzd|tzdtzzz }tdtdzzd|tdzzz }t dtzd|tzzd|tzdtzzz tdtdzzd|tdzzz z}tdtzd|tztdzzz }t|t|tdksJt|t|tdksJt|t|tdksJt|t|tdksJt|t|tdksJt|t|tdksJt|t|tdksJdS)Nr0rrrrur) rrrWrrr rrr)r0eq1eq2eq3eq4eq5rrs r8test_homogeneous_functionr:Gsf A a!!A$$h--C qsQqqttVn  C a!fXa^^# $ $C qsQqqttV|a!fqsl+ , ,C qAv!!A$$' " # #C qsQqqttV|a!fqsl+c1QT6AaaddAgI2F.G.GG H HC qsQqqttVad]# $ $C S!QQqTT * *d 2 2 2 2 S!QQqTT * *a / / / / S!QQqTT * *d 2 2 2 2 S!QQqTT * *a / / / / S!QQqTT * *a / / / / S!QQqTT * *a / / / / S!QQqTT * *d 2 2 2 2 2 2r7cjtdtzdttzzdzztdtzdttzzdzz}}||z }t|t |z}t |}|}t |}tt|}t |}dtddf} | |tt| ksJ| |tt| ksJ| |tt| ksJdtzttz } t | } dtzdztz } t | } dtzdttzzd zdtzdttzzdzz }t |}dtzdttzztdzdtzdttzzdzz }t |}| | ttJ| | ttJ| |ttJ| |ttJdS) Nrurrr4r  rir?) rrWr0rrexpandr'rr _linear_coeff_matchr)rdratr5obj1r6obj2r7obj3rr8obj4r9obj5robj6robj7s r8test_linear_coeff_matchrHYs acAaddFlQ AaC!AaDD&L2$5!6qA A#C c((S&& &C c " "D C c " "D c#hh--C c " "D hsA C  # #C1 . .# 5 5 5 5  # #C1 . .# 5 5 5 5  # #C1 . .# 5 5 5 5 Q3!*C c " "D Q37A+C c " "D Q31Q44AAE!Q$O P PC "-AaDD!9 B B Bc I I I I I Ir7c l td\}}tttdtttzz }t |dksJt tt|tdzdz dzz|tztdzdz dzzzt tdzz}t|d |ksJt||d ksJt tt|tdzd zdz tdzdzdz ztdzdzz dzz|ttdzd zdz ztdzdzdz z tdzzzt tdzz}t|dd |ksJt tt|tz|zt tdzz}t|dd|ksJt||d ksJdtdzztttdzdtztttzzdttzz }t |dksJt tt|td z dz tdzzdzz|tzzt tdzz}t|d |ksJt||d ksJtttdttttzzttz}t |dksJt tt|td zdz tdzdz z dzz|tztdz dz dzzzt tdzz}t||ksJtttdtttztttzz }t |dksJt tt|td z dz tdzdz zdzz|tztdzdz tdzdz ztdz z dzzzt tdzz}t||ksJtttdtttzz}t |dksJt tt|tdzdz tdzdz z dzz|tztdz dz dzzzt tdzz}t|dd|ksJt||d ksJdS)NzC1 C2ru)2nd_linear_airyrrrr?rrr<Trrr)r=x0)r=r)rm2nd_hypergeometric2nd_hypergeometric_Integralr)rmrr)rrr4) rr0rWrrr rrr r )r[rr\rs r8test_2nd_power_series_ordinaryrVs W  FB 11a1QqTT6 !B   O O O O O QqTT2q!tAvz?RT1a47Q;%77!AqD''A B BC "6 7 7 73 > > > > r3  9 , , , , QqTT2Az!|q1uqjl2a!eaZ?!CD a1q51*R-1q51*Q,.15 67 AqD''  C "62 > > ># E E E E QqTT2a4"9qAww& ' 'C "6! < < < C C C C r3  9 , , , , ad(QqTTYYq!__ %!QqTTYYq\\(: :AaddF BB    ! ! ! ! ! QqTT21uQwA~)*RT1AadGG; < > > > r3  9 , , , , 11a1addiill+ +add 2B   K K K K K QqTT2q!tAv1Q*+bdQTE!GaK.@@1QT77J K KC "::     11a1Q4499Q<< '!AaDD& 0B   = = = = = QqTT21uRx!Q$q&(1,-dAqDGad1f$qs*Q./023AqD'': ; ;C "::     11a1QqTT6 !B   O O O O O QqTT2q!tCx!Q$q&(1,-1q!teBhl0CCa1ggM N NC "6! < < < C C C C r3  9 , , , , , ,r7ctd\}}}tdztttdzdtztttzz dtzdzttzz}t tt|tdzzdtdzzdz dtdzzzdtzz dzzt tdzz}t |d |ksJt||d ksJdtdzztttdzd tdzztttzz dtdzzdzttzz}t tt|ttztdzd z tdzdz ztdzdz ztzdzzt tdzz}t |d |ksJt||d ksJtdztttdztdztttzz tdzdz ttzz}t tt|tdz dz dtdzzdz z tdzd z z tdzdz ztdz zdzztz |tdzztdz dz tdzdz ztdz zdzzzt tdzz}t ||ksJt||d ksJtdztttdzttttzztdztddz ttzz}t tt|tdzd z tdzdz z dzzttz |ttztdzdz tdzdz z dzzzt tdzz}t |d |ksJt||d ksJttttdztttz|tzttzz }t tt||dztdzzdz |tdzzdz zdzzt tdzz}t |ttd |ksJt||d ksJtttddtz tz tttzz|tz ttzz}t tt|| tdzz|dz z|dz z|dz z|dz zdz |tdzz|dz z|dz z|dz zdz z|tdzz|dz z|dz zdz z |tdzz|dz zdz z|tzz dzzt tdzz}t |ttd |ksJt||d ksJdS)NzC1 C2 arurrir r?r2nd_power_series_regularr<rQrrirP<rM@i@8i@$) rrWr0rr rrr rr )r[rryr\rs r8test_2nd_power_series_regularr^s ""IBA AqttyyA !A#qttyy||"4 4!a1~ EB QqTT2ad7C1HQJ1a4/!A#59:Qq!tWWD E EC "5 6 6 6# = = = = r3  9 , , , , 1a411a !1QT61Q4499Q<<#8 8AadF = Q44< B QqTT2d1gg:q!tBwAa/!Q$q&811DqHI!K q!tGadU2X1R'!A#-1 2356q!tWW= > >C "::     r3  9 , , , , AqttyyA !QqTTYYq\\"2 2adXa^^6KQqTT5Q QB QqTT2q!tBwAa'!+,T!WW4 477 AqDHq!tAv%)*+-.q!tWW5 6 6C "5 6 6 6# = = = = r3  9 , , , , 1Q4499Q?? QqTTYYq\\ )AaC!H 4B QqTT2q!tAqDy|a1fQh.23a1gg= > >C "add!; < < < C C C C r3  9 , , , , 11aQUAIqttyy||3 3qsAaDDj @B QqTT2r!Q$wAA.A6A>uD !Q$AAA&s*+-.q!tVQU^QU-CB-FG !Q$AqQ3!"#$&'1gg. / /C "add!; < < < C C C C r3  9 , , , , , ,r7ctd}td}ttdtd|z|ksJdS)Nrrr[r)r r$r)rrs r8test_issue_15056r`sRs A B !&,,$1B1B"B C Cr I I I I I Ir7c t tz dtzttzz ttz ttttz}tt tdztzzt tdztzt tt tdz tz ztz tzz}t||dksJttt t tzzz}ttttzttt tdz}t||ttdksJdS)NrurQ) r[rWr0rrrrr rrr\rs r8test_issue_15913rcs" Q1QqTT AaDD :addA#6#6 6B SA]] SA]]8BsAqD519~~4Ea4G+K+KK KC r3  9 , , , , r#qbd))| C AaddFA  11a 0B r3! % % 2 2 2 2 2 2r7c^ttdttddS)Nc@ttttt ttgttt tt tgSr4)rr0rWrr1r2r6r7r8rrz"test_issue_16146..sRvqttyy||QqTTYYq\\&BQqTT1Q44QRSTQUQUDVWWr7cttttt ttgttgSr4)rr0rWrr1r6r7r8rrz"test_issue_16146..s=vqttyy||QqTTYYq\\&BQqTTFKKr7)r*r5r6r7r8test_issue_16146rgs0 :WWXXX :KKLLLLLr7cLttdz tttz}tttt}t ||ksJttt tz tttdz}ttttttzztttt th}tt ||ksJttdzdttzz dztttdz}ttttttzzttdzzz}t ||ksJdS)Nrur?r) r0rWrr r[rrrrrrbs r8&test_dsolve_remove_redundant_solutionsrisQ A$$q&!A$$))A,, B QqTT2,,C "::     A$$s1vv+! !Q (B addBAI  1Q44Q 0 0 1C vbzz??c ! ! ! ! A$$'!AaDD&. AaDDIIaOO +B QqTT219r!Q$w& ' 'C "::      r7ctdt\}}td}tt |||||||ztt |||||||zg}t |}t ||dddgfksJdS)NzA BrrTr)rrr r rrr )rBrr\rs r8test_issue_13060rls 5h ' ' 'DAqs A Z!a !!A$$qqtt) , ,bAAaDD!1D1Daadd11Q44i.P.P QB **C r3  D1a&> 1 1 1 1 1 1r7ctd\}}td}d|zt|dz z||zd|dzz|dz z|d|z|dz zz dzzt|||dfzz}t |dd }|d d d ksJtd \}}t |d }t |t tdksJdS)NzN srhog@r?rruTr;rrtermsrzC1,C2r<a Eq(rho(s), C2*(1 - 4.0*s**4*sqrt(N - 1.0)/N + 0.666666666666667*s**4/N - 2.66666666666667*s**3*sqrt(N - 1.0)/N - 2.0*s**2*sqrt(N - 1.0)/N + 9.33333333333333*s**4*sqrt(N - 1.0)/N**2 - 0.666666666666667*s**4/N**2 + 2.66666666666667*s**3*sqrt(N - 1.0)/N**2 - 5.33333333333333*s**4*sqrt(N - 1.0)/N**3) + C1*s*(1.0 - 1.33333333333333*s**3*sqrt(N - 1.0)/N - 0.666666666666667*s**2*sqrt(N - 1.0)/N + 1.33333333333333*s**3*sqrt(N - 1.0)/N**2) + O(s**6)))rrrrrrr-str)Nr rneqnmatchr[rrs r8test_issue_22523rts5 5>>DAq 5//C a%QU CCFF "a1fa!enAaCQK!7K&K SSVVaV $ $&% %C 4e 4 4 4E , -g 6! ; ; ; ; W  FB 6 7 7 7C c??j.K*L*LMM M M M M M Mr7ctdt\}}tdd\}}}}}d\}}}}t|t|||dz|||zz|||||z zz d}t|t|||dz|||||z zzd}||g} t | ||||g|dd|||dd|dd |||ddi \} } | t||t d t d z t d d t dzzztt d|zt d t d z zdz zdz dt dzt d zt t d d zztt d|zt t d d zzdz zdz zksJ| t||t d d ztt d|zt d t d z zdz zd z d t d z tt d|zt t d d zzdz zd z zksJdS)Nzx1, x2rz t k1 k2 m1 m2Tr.)r?r?r?r?rurr?rdrr r[)rrr rrrrr) rx2rrrrrr5r6eqsx1solx2sols r8test_issue_22604r|-s XX . . .FB===Ar2r2NBB RRRUUAq!! !Brr!uuH ,r22a5522a55=/A A1 E EC RRRUUAq!! !B11 $6 6 : :C *CC""Q%%A1aAARARSTUVAWAWXY8:1aAARARSTUVAWAWXY8[\\\NUE Brr!uud1tAww;//bAd1ggI1EFs4PQ77ST9UYZ[^bcd^e^eZeUfUfKfghKhGiGiijll477 T"XX-tDGGaK/@/@@T!WWQYtTXYZT[T[^_T_O`O`E`abEbAcAccdffghh h h h h Brr!uutAww{CQ $q477{:K:K0KA0M,N,NNrQUVY]^_Y`Y`U`befjklfmfmnofoptuyz{u|u|@AvAqBqBgBCDgDcEcEUEFHUHHII I I I I I Ir7c  ttttdttdzzdttzdz z ttttdttdzzdttzdz z fD](}dt |ttvsJ)dS) Niruii rrr4rC)r r0rWrr)des r8test_issue_22462r;s qttyy||S1q[3qtt8D=8 9 9 qttyy||R!aZ!AaDD&(2 3 3555l2qtt444444455r7c:td}td}tt|z ||z||zd}t |dksJdS)NrWrr)rtnth_order_reduciblerr)rrr rrrrWrr\s r8test_issue_23425rBs A A QTE!!A$$))++""$$ $qqttyy{{ 2A 6 6B    ; ; ; ; ; ; ;r7ztoo slow for @slowctd}td}||dz|||dzdz}t|||JdS)NrWrrru1)r rrrrs r8test_issue_25820rKsks A A 1q11a 2 %B "aadd   ' ' ' ' 'r7N)wsympy.core.functionrrrrsympy.core.numbersrrr r sympy.core.relationalr sympy.core.singletonr sympy.core.symbolr r$sympy.functions.elementary.complexesrr&sympy.functions.elementary.exponentialrr%sympy.functions.elementary.hyperbolicr(sympy.functions.elementary.miscellaneousr(sympy.functions.elementary.trigonometricrrrrsympy.integrals.integralsrsympy.polys.polytoolsrsympy.series.orderrsympy.simplify.radsimprsympy.solvers.oderrrsympy.solvers.ode.subscheckr sympy.solvers.ode.oder!r"r#r$r% sympy.solvers.ode.nonhomogeneousr&sympy.solvers.ode.singler'sympy.solvers.deutilsr(sympy.testing.pytestr)r*r+r,sympy.utilities.miscr-C0r[rrrC5C6C7C8C9C10urWrrr0r1r2r9rbrgr}rrrrrrrrr rrrr r+r-r2r:rHrKrNrVr^r`rcrgrirlrtr|rrrr6r7r8rsZBBBBBBBBBBBB333333333333$$$$$$""""""////////99999999========777777999999KKKKKKKKKKKK......&&&&&& ******433333HHHHHHHHHHHHHHNMMMMM777777++++++::::::::::::++++++/6gg.>.>+BBBBB WW4 ( ( ( 1a HSMM HSMM HSMM%%% #3#3#3L888 D?D?D?N\\\~i$i$i$V,*,*,*^Y?Y?Y?vFFF4777B///***lMlMlM^===,;;;999kkk I I I777 GGG 333$777@FFF JJJ(/-/-/-d$-$-$-NJJJ 333MMM   222MMM, I I I555;;;(((((r7