wg| 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>dAdBZvyC)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:zTrealfghcbtt5tddddy#1swYyxYwN)r* ValueErrorr$e/home/mcse/projects/flask/flask-venv/lib/python3.12/site-packages/sympy/solvers/ode/tests/test_ode.pytest_get_numbered_constantsr90s'  %t$%%%s %.c 2ttjt}t|d}|d}|jj dd}|d}|j j dj dd}idtttttdtzdttttd ttttd ttttd ttttd ttttd ttttdttttdttttdttttdtttttdtzdttttdtttdttd|z|tttdttttdttttdtttttd|z |tttz fzdttttttd|ttfttdtzttttttttd d dd}||k(sJt|dttttk(sJy)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_hintrb5s 11B BU #F & 'D q ! B E FD ((--    "1 %C$b1rHQN/B&C$b1rlR!b\ AaD" *2adB<  0AaD"  R!b\  3BqtRL !"QqT2, b1rl r!A$Xa^(;< R!b\ bhq!nx2&FAaD!QSUV R!b\ 3BqtRL ) rMr@rRrKrArlrSrTrUrmrErVr>rVrSrUrmrTrLrGrBrOrHrJrKrMr@rRr>ctdS)N notarealhintr<rr\sr8z%test_dsolve_options..svb~>r7ctdS)N Liouviller<rprqsr8rrz%test_dsolve_options..svb{;r7r?)rWr0rrsortedkeysr(r r[hasrr*r5rr)abcrw Integral_keysr\s @r8test_dsolve_optionsr}asU 1Q499Q<!A$ BrAr.Ar'A +D'M !&&( t ## # W:2qt, ,, , V91Q4A && & "6 "b1r!tn 44 4 Y<< '' ' [>\ )) )~!!(++ +~!!(++ +-.228<< <;<@@JJ J;<@@JJ J""8,, , " # ' ' 11 1 ! " & &x 00 0 @ A E Eh OO O @ A E Eh OO O ! " & &x 00 0 !&&( t ## # W:2qt, ,, , V91Q4A && & "6E 2b1r!tn DD D Y<< '' ' [>\ )) ) [>Q{^ ++ + 7 8 45 66 6 7 8 45 66 6~!!(++ +~!!(++ +-.228<< <;<@@JJ J;<@@JJ J""8,, , " # ' ' 11 1 ! " & &x 00 0 @ A E Eh OO O @ A E Eh OO O ! " & &x 00 0 !&&( } ,, , :>? :;< !A$))A,1Q47* 7 ?qtyy|a!ai'f5 66 6 !A$11%A A6:AaD, . 1Q4"x#a&++a.1!4HQN#!$%&((),hq!n_)=> ? ?? ?r7c tttjtdttdk(sJtttttdk(sJtt ttjtdttdk(sJtttjtdzttdk(sJtt ttjtttzttt}tttjtttzttttzztttzz tt}tttjtttz ttttz ztttz z tt}|dk(sJ|dk(sJ|dk(sJtdtzttzttjtzd tzttdzzzt tz ttd k(sJd tdttjtzttdzztttz z d td zzzttvsJt tdtd}tttjt|ttz|tzttzzz dttz|ttz|tzttzzz ztttjtz|ttz|tzttzzz ztzttdk(sJd}tttttztttttz|k(sJtttttztttttzttd|k(sJtt dtd zzttjtzdttdk(sJtt dttd zzttjtzdttdk(sJtt tttttttzz dttdk(sJttt ttdttdtsJttt ttjtdttdjgdk(sJtdtzttzttjtzd tzttdzzzt tz ttdd}t|jgdk(sJtttjtt tttz z tt}|dk(sJy)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_minus2ctttttj tj tztttSr4)rrWr0yrr6r7r8rrz#test_classify_ode..s9|A!Q Q0D0D 1 - aG 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_oders ! !Q1 .    !ad #'R RR R 1Q499Q<+QqT 27= == =& ! ! a1 .3> >> >( R! ! qt+Q/16AQqTYYq\!A$&1ad2QqtV;QqTBAQqTYYq\!A$&1ad2QqtV;QqTBA F F F F G G G G G G G G  !AaD11Q!a/#a&8!A$  N NN N $Qqtyy|^adAg-!Q61r7BAaDI JJ J : s A ! ! a!fqs1Q4x&781QqT61QqT6 !AaDD<adiilNAadFQqS1X$5679:;<=aD B 5 55 5  OC$ QqTAXq)D1qM9 :c AA A QqTAXq)D1qM91Q4! #&) ** * 1QT6!A$))A,.2AaD 9 6 66 6 1QqT1W9QqTYYq\115qt <    4!a=1Q472A6! =) ** * l2adA;!4@$ GG G ,r!A$))A,2AaDtDIIK LE FF FQqS1Xadiil*a!eQqT1W_.N|B!#>r7rc:tttSrrrsr8rrz'test_classify_ode_ics..Xrr7c:tttSrrrsr8rrz'test_classify_ode_ics..\rr7c:tttSrrrsr8rrz'test_classify_ode_ics..`rr7c:tttSrrrsr8rrz'test_classify_ode_ics..drr7rdc:tttSrrrsr8rrz'test_classify_ode_ics..urr7c:tttSrrrsr8rrz'test_classify_ode_ics..yrr7c:tttSrrrsr8rrz'test_classify_ode_ics..}rr7c:tttSrrrsr8rrz'test_classify_ode_ics..rr7c:tttSrrrsr8rrz'test_classify_ode_ics..rr7c:tttSrrrsr8rrz'test_classify_ode_ics..rr7c:tttSrrrsr8rrz'test_classify_ode_ics..rr7c:tttSrrrsr8rrz'test_classify_ode_ics..rr7c:tttSrrrsr8rrz'test_classify_ode_ics..rr7c:tttSrrrsr8rrz'test_classify_ode_ics..rr7c:tttSrrrsr8rrz'test_classify_ode_ics..rr7) r0rWrr*r5r1rsubsrrrrs@@r8test_classify_ode_icsrHs 11a1Q4 B a&C :>? Q4)C :>? Q4)C :>? Q7A,C :>? Q41,C :>? Q41,CQqTs# Q4)CQqTs# Q499Q<  Q "A &C :>? Q499Q<  Q "A &C :>? Q499Q<  Q "A &C :>? Q7<<?  1 %q )C :>? adAq ! & &q! ,a 0C :>? Q499Q<  Q "AaD )C :>? Q499Q<  Q "AaDIIaL$5$5a$; ? Q499Q<  Q "A &C :>? Q7<<?  1 %q )C :>? adAq ! & &q! ,a 0C :>? Q499Q<  Q "AaD )C :>? Q499Q<  Q "AaD )CQqTs# Q499Q<  Q "A &CQqTs#r7ctdd\}}}}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 | |!k(sJt|||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 |"|#k(sJt|||||||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 |$|%k(sJt|||||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 |&|'k(sJt|||||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 |(|)k(sJy)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_sysoderst4JAq!Q)01U_c)d&BBBBB1x@Aq!Q1!()E8!TBBBBiX.GAq! A ad1BT!A$q\ b#a!f+a!QqTl* +R1QqT71QqT!A$<3G-H ICAaD! a!QqT1qSTVWXYVZ\]R^`adeghijgkmncoqr!aL!a1q\1q!A$lA1Q4|Q0H]dqT1Q4Lua!adm^C!A$K-G*UVWXUY[\J]-]`abcdebfaDa ` AaD'`qtQ'`(-)45aD!QqT13E GD 3 4 '' ' b!A$'!A$qt)# $bQqT!A$Y&7 8CAaD! a!QqT1qSTVWXYVZ\]R^`adeghijgkmncoqr!aL"QqT'AqtQ<Q!aL"QqT'AqtQ>AaD "QqT(@RSTUVWXUY[\]^[_abSc@cqtQA*+QqT1q1adA,AqQRtUVaqT1WHz!A$"">$/0tQ!a.@ BE 4 E )) ) r!A$qt)CFAI% &2adAgc!fai.?(@ AD !QqT1!uSVQYQRTUVWTXZ[P\^_bcefgheiklamop!aL!a1q\1q!A$lA1Q4|adU3q6ST9_WXZ[\]Z^`aVbde1g!A$!ESTUVSWRWXYZ[X\R\]`ab]cef]fRfqtQSA$'#a&!)+j1q.AARCNOPQdTUWXYZW[]^M_ aE 4 E )) )r7c F"tttjtdttz z tttddit ttt dtzdzk(sJtttjtttz tttddit ttt tk(sJtttjtttz ttttjtjtddit ttt tk(sJtttjttttztttddttjtjtddit ttttttzk(sJtttjtttz ttzttjtttz ttz gttttgtddtddit ttt tttzt ttt tttzgk(sJtdzttdzztz jt}t|tttddit ttt tdz  tz t ttt tdz tz gk(sJt|ttttjtjtddit ttt ttjz  tz t ttt ttjz tz gk(sJtttttttzttdzz ttjtzz }t|tttddiddt ttttzttdzdz ztddk(sJt|tttddidd t ttttzttdzdz ztddk(sJtt tttt tzgttgtgtdditdik(sJtt tttttzt ttzzgttgtt gtddtt"dz ditdt dik(sJtt tttttzt ttzzgttgtt gtddttjtjtdditdt dik(sJtt tttttzt ttzzgttgtt gtddit dik(sJt%tttjtttjtdtz ztddttjtjtddit ttdt tttdzd z tdz z hk(sJt%tttjtttjtdtz ztddit ttdt ttt tztdzd z zhk(sJt'd \}}}t tt||zt |tzz| |z|t |tzz| |zz dz zz }tt ttjt|ttzdtt|z z ztttd|i|k(sJt%tttjtttjtdtz zttjtjtddtddit ttdt tttdzd z hk(sJt%tttjtttjtdtz ztddttjtjtddit ttdt tttdzd z hk(sJt)t*d t)t*d t'd\}}}t tttt tzzt,tdzzzt.tdzzz|tdzzd|zz zg}ttg} tt t,t.g} tddttjtjtddt|j|ddt|j|ddi} tddttjtjtddt1ttjtdt|dt1ttjtdt|di} t|| | | } t|| | | }| |cxk(r4tdt dt,|dz|zd|zz t.| |zd |zz ik(sJJtdtdi}tttjtttz tt|t tttdt tzk(sJttjtjtdttjtjtdtdtdi}tttjttttztt|t tttdttzttjtjtdttzzk(sJy)Nr?rurdrrLF)rer=rirTzK r f0c tttttt tzt ttzzgttgtt gtddttdiS)Nrr?) r%r r0rWr[rC2rr r6r7r8rrz test_solve_ics..sZy"QqT2c!f9r#a&y3H*I)JQqTFUWY[T\_`ab_cefhijlhmop^qrr7c ttttt t tdzz gttgt gttj tjtddiS)Nrur)r%r r0rWrr[rrr6r7r8rrz test_solve_ics..s[y"QqT4QT ?*C)DqtfrdUVWXUYU^U^_`UaUfUfghjkUlnoTpqr7zEI q L)rr0rWrr rrrrrr1r Halfr r%r[rr setrr*r5C3C4r)r\Krf0solEIrLsolsfuncs constantsics1ics2solved_constants1solved_constants2res r8test_solve_icsrs !A$))A,1Q4'1AaD!9 = 1Q4a!eai! "" " !A$))A,1%qt!A$ ;r!A$A?O OO O !A$))A,1%qt!A$))A,2C2CAq2I11M NRTUVWXUY[^_`[aRb bb b !A$))Aq/AaD(!A$QqT1 ! ! !Q5$ %(*1Q4Q#a&(A BB B 1Q499QUW_`acdWe;f gg g "ad qT!H; 7:!1Q4QS* UZ] ]] ] vadiilAaDIIaOA$56QqTYYq\=N=NqQR=STUWXYZW[\]<^_ ` AaD!b1qAvz*+ ,, , vadiilAaDIIaOA$56QqT!QqTYYq\EVEVWXYZE[\]<^_ ` AaD!b1qAvz*+ ,, , :rs :qrx HB1 qtR"Q$YAqD(2ad72Qq!tVRU^C D ED qTFERR I aD! ! ! !Q ! !Q ! !Q D aD! ! ! !Q QqTYYq!_a#Q QqTYYq!_a#Q (D "$y$?!$y$?  1  A A AqDFAbDM QBqD!B$K 6    Q41,C !A$))A,1%qt 5AaD!A$s1v+9N NN N Q499Q<  Q "AaDIIaL$5$5a$;QqT1Q4 HC !A$))Aq/AaD(!A$C 8 1Q41c!fqtyy|00A6s1v==> ?? ?r7c td}td}td}td|zt||z||dk(sJt|t |||zd|z||zzt ||z z ||dk(sJt|dz||j ||z|t |||zz||z ||dk(sJtt |t||z||||dk(sJtt |t |t||z||z|||dk(sJtt ||||||dk(sJtt ||||t |||z||dk(sJtt ||||t |||z||dk(sJtt |t |t||z||z|||dk(sJtt |||z|||dk(sJt|t t |||dzz||z||dk(sJtt |t |t||z||z|||dk(sJtt ||||||dk(sJtt |t||z||||dk(sJtt ||||t |||z||dk(sJtt |t ||||z|||dk(sJt|t t |t |||dzz||z||dk(sJy)Nr0r1rWrrr?ru)rr r(rrrr)r0r1rWs r8test_ode_orderr=s A As A QqSQqT]AaD )Q .. . QtAaD!}_qs1Q4x/#a&(:AaD AQ FF F QT!A$))Aq/)Ad1Q4mO;adBAaD IQ NN N T!C!I+q!,ad 3q 88 8 T!D3qt9a33Q71 >! CC C T!A$1%qt , 11 1 T!A$1%d1Q4m3QqT :a ?? ? T!A$1%d1Q4m3QqT :a ?? ? T!D3qt9a33Q71 >! CC C Z!A$*AaD 1Q 66 6 Qs:a!aiA6771 >! CC C Z*Qs1Q4y[!Q"? ?CQqT Ja OO O Z!a+QqT 2a 77 7 Z#ad) Q2AaD 9Q >> > Z!a+JqtQ,??1 F! KK K Z*QqT1a"8 8!.us015r7)rrrrWrrrrrr r r0rr rrr*r5)rys r8test_homogeneous_orderrWs S1XAaC0!Q 71 << < QTCF3q6M11a 8 @@ @ QUQs1Q3xZ/A 6! ;; ; acDA1$551c!fs1vo8NN Ax1~  d1gqj (*+,a 14> >s A Qs1Q3xZ!CF(2Aq 9Q >> > QqTYYq\1a 0 88 8 adiil]Q.1 5 == = QqT1a ( 00 0 Q1a4[!Q / 77 7 QUA &" ,, , QT1 % -- - :56r7ctdtz}t|t|}}|j t r|j t rJyy)Nr?)rrWrrrxr)rirs r8)test_noncircularized_real_imaginary_partsrxsA QqS A a5"Q%qAe u. . r7cpdttdz z}|t|ttdz k(sJy)Nr?rur)rrWr)rs r8$test_collect_respecting_exponentialsrs2 c!A#h,C '3AaC) )) )r7c"ttttddik(sJttdtzt dztdt dtzt dztdtzt dzhdk(sJtttt tztddik(sJt ttt tztdzt tztdzttztttztth}ttttdztzdzztd|dk(sJttttdzztttzzttztd|dk(sJtt dtzttztdztzdzztdt dtzttztdzt dtzzttzt tt dtzztdzt tzt dtzztt tzt dtzztt dtzzttzhdk(sJtdttz tddik(sJttttddik(sJtdtztdztzdzztddtztdtzztdzdtzzhdk(sJtttztddik(sJtt tt dtzdzztdt ddtzzhdk(sJttttdztzdzztdtt tztttztdzt tztdzttzt ttthdk(sJttttttzztddik(sJtttttdtzzztddik(sJtttttztddik(sJttdzttzt tztttzztztdtdzt tzt tztt ttjt tttztttt tzttztt tztt tzt tztttzt tt tztdzt tzttzh dk(sJtd tzttdz zttt tdz ztttdz zt tdz ttdz hdd k(sJtdtztztddtztdtzzhdk(sJtdtzt dtzztddtzt dtzzhdk(sJtt t tz tddik(sJttd tdtjhdk(sJtd t tztdt thdk(sJtt ttztdt ttzhdk(sJttttdt ttthdk(sJtt ttdt ttthdk(sJtd d t tzzdttzztdtjt tttt thdk(sJttdztdtjttdzhdk(sJtdtzt tzt t ztdtt tzt tt t hdk(sJtdt dtzzttztdt dtzttzt tt dtzzhdk(sJttttz tdtjtt ttthdk(sJttdzdtzztdtjttdzhdk(sJtd tzttztdtt tztttzt ttthdk(sJtttdtzztdtt dtzzttdtzzt dtztdtzhdk(sJttdzt t ztdtt t ztdzt t zt t hdk(sJtdt t ztdzt t zz tdtt t ztdzt t zt t hdk(sJtt dtztdzztdtjttdzt dtzhdk(sJttt t ztdtt t zt t hdk(sJttt dtzztdtjtt dtzhdk(sJtttt t ztdt tttt t hdk(sJtt ttdt thdk(sJttjt dtzdz z tdtjt dtztdtzhdk(sJtt dtztjt dtzdz zztdt dtztdtzzt dtzt dtzzt dtzhdk(sJtdtzttzt tztdtjtt ttthdk(sJtt tdz t dtzdz z tdt tt dtztttdtzhdk(sJtt tdztddik(sJtdtdzztddik(sJy)NtestFrurT)rtrialsetr?rr)rr r r)r&r1rWrrrrrrrr Onerr)ss r8$test_undetermined_coefficients_matchr sa +AaD! 4 GG G +C!d1g ,> B 1tAw QqS47]!3 4 6 66 6 ,CF3q6M1 =    Q3q61a4A;1SV Qs1vXs1vFA +CFAqD1HqL,A1 E1% && & + Aq!t c!fQhQ' ,8> > +AqD! 4 GG G +CF3qsQw<,? CCAaCL>2 33 3 +CFAqD1HqL,A1 EAc!fHaAh1SV  1SV SVSV$% & && & ,CFAAJ,? C    +CFAAaCL,A1 E    +CF3q6M1 =    + 1SV CFQs1vX%)1 1a4A;s1v#5q#a&!%% As1v s1vqQxA#a&!CF(3q6/ #a&#a&Q-Ac!fSV!3#5 6 66 6 ,AaCAE NA >s1q5z\1SQZ<QUSQZHC   ,AqDFA 6AqD!AqD&>2 33 3 +AqDQqSM1 =AqDQqSM?3 44 4 +CGAIq 9    ,AaD! 4AEE7+ ,, , +Bs1vIq 9CF8, -- - +C!Ha 8C!H:. // / +CFA 6CFCF#34 55 5 +CFA 6CFCF#34 55 5 +A#a&L1SV8,CQ GAEE3q63q63q6#BC DD D +AqD! 4AEE1ad#34 55 5 +AaCAJaR,@! DAc!fHc!fc1"g#>? @@ @ +Ac!A#hJs1v,=q AC!HSVOSVC!H_#EF GG G +AAJ :AEE1c!fc!f#=> ?? ? +AqD1Q3J :AEE1ad#34 55 5 +AaCAJ :Ac!fHaAhAA#GH II I +Ac!A#hJ : s1Q3xZ3qs8S1Xs1Q3x 8 : :: : ,AqDaRL! <Ac1"gIq!tCG|S!W#EF GG G +Ac1"gI1S!W ,Da HAc1"gIq!tCG|S!W#EF GG G +C1I1,? @@ @ +Ac1"gIq 9Ac1"gIsA2w#78 99 9 +AAaCL! <AEE1c!A#h#78 99 9 +CFS!W,> > +CFA 6CF8, -- - ,AFFS1XaZ,? CAEE3qs8S1X#>? @@ @ , AaC!&&3qs8A:%& 3qs8C!H#4c!A#hs1Q3x6G AaC#    ,AaC#a&L3q6,A1 EAEE1c!fc!f#=> ?? ? ,CF1Hs1Q3xz,A1 ECFC!Hc!fc!A#h#GH II I +C1Iq 9fe_ LL L +A1Iq 9fe_ LL Lr7cddlm}t|ttt ttzt tzzztt ttzt tzdz}t |t tdk(sJtdzt tdzzt tjtzdtzt tzz }t |t tdk(sJy)Nr)Aru)rmrLrJrCrkrQrNrrrMrKrArlr) rmrLrBrHrOrQrNrMrRr@) sympy.abcrrWrr0r)rr\s r8test_issue_4785_22462rs QD1qM!AaD( ) )D1qM 9AaD @1 DB AaD !&F FF F Q$1q.!A$))A, &1QqT 1B AaD !&= == =r7cttdtttt j tt ttt zz ttdddddk(sJttdtttj ttt dddddk(sJy)Ncttttj tttttzz ttSr4)rr0rWrrr6r7r8rrz!test_issue_4825..s.va1gll1o!Aq' &A1Q4Hr7Trrr6)rVrUrcttttjtttSr4)rr0rWrrr6r7r8rrz!test_issue_4825..svadiilAaD9r7)r*r5rr0rWrrr6r7r8test_issue_4825rs :HI !Q Q!AaG)3QqT Er: ;; ; :9: ! ! ad 6r: ;; ;r7cddlm}tttzt t zzttt zt tzzcxk(rttzt t zzk(sJJtt tzztt zz}|tt tgdD]*\}}}t||tzz|t zz|k(r*Jy)Nr) variationsr)sympy.utilities.iterablesrr"r[rWrrr)reryrzr{s r8'test_constant_renumber_order_issue_5308rs4 RTBqD[ )"Q$A+&  1r!t     BF R!VAr2rlA.91a AEAE!23q8889r7ctd\}}}}||z||z||zzg}t|d||zk(sJt|d|gt|zk(sJt|d|gtgt|zk(sJt|||gt|zt|zt|zzgk(sJt|||gtdt|zt|zt |zzgk(sJy)Nze1:3 x yr) variables)r newconstantszC3:5)rr"r[rrr)e1e2rWrexprss r8test_constant_renumberr #s:&LBAq T2a4"Q$; E U1X &"Q$ .. . U1X! 5A == = U1X!B4 HBqD PP P Uq!f 5"Q$1r!t 9L LL L Uq!f76? SXZ[\X\^`ab^beghiei^iWj jj jr7c tdd}td}td}t||j|d|dz||zz ||}t |j Dcgc] }|j jds|"c}dk(sJttttztttzzttztthtttzttzk(sJttttztttzztttzzttthtttztttzzk(sJttttzthtttzk(sJtttzt ztt httzk(sJtttzt#ttdd fzthttzk(sJycc}w) NrTr.rwrCr?ru)r rrrlen free_symbolsname startswithr#r[rrWrrrrrr)rrr"rr s r8test_issue_5770r(.ssAs A A 11a1a4!9,ad 3C 3++Faqvv/@/@/EF G1 LL L CFRAY.A6R A 3q6 #a&   3q6 Bs1vI-3q6 9BB< H 3q6 Bs1vI   BF bT *bQi 77 7 B RG ,Q 66 6 B!aAY!77" >"q& HH HGs . I!I!cttt ttztJttttz ttJyr4)rrrWrrr6r7r8test_issue_5112_5430r*=sA c!fWuQx/ 3 ;; ; QQZA . 66 6r7cBtdttfdy)Nr0crttjtdztdS)Nrufdsjf)rrWrr0sr8rrz!test_issue_5095..Ds$vadiilAoqtWEr7)rr*r5r.s@r8test_issue_5095r/Bs A :EFr7c 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|tdk(sJt|t|tdk(sJt|t|tdk(sJt|t|tdk(sJt|t|tdk(sJt|t|tdk(sJt|t|tdk(sJy)Nr0rrrrur) rrrWrrr rrr)r0eq1eq2eq3eq4eq5rrs r8test_homogeneous_functionr7Gs A a!A$h-C qsQqtVn C a!fXa^# $C qsQqtV|a!fqsl+ ,C qAv!A$' " #C qsQqtV|a!fqsl+c1QT6AadAgI2F.GG HC qsQqtVad]# $C S!QqT *d 22 2 S!QqT *a // / S!QqT *d 22 2 S!QqT *a // / S!QqT *a // / S!QqT *a // / S!QqT *d 22 2r7ctdtzdttzzdzztdtzdttzzdzz}}||z }t|t |j z}t |}|}t |}tt|}t |}dtddf} |j|tt| k(sJ|j|tt| k(sJ|j|tt| k(sJdtzttz } t | } dtzdztz } t | } dtzdttzzd zdtzdttzzdzz }t |}dtzdttzztdzdtzdttzzdzz }t |}| j| ttJ| j| ttJ|j|ttJ|j|ttJy) Nrurrr1r  rir?) rrWr0rrexpandr'rr _linear_coeff_matchr)rdratr2obj1r3obj2r4obj3rr5obj4r6obj5robj6robj7s r8test_linear_coeff_matchrEYs+ acAadFlQ AaC!AaD&L2$5!6qA A#C c(S& &C c "D C c "D c#h-C c "D hsA C  # #C1 .# 55 5  # #C1 .# 55 5  # #C1 .# 55 5 Q3!*C c "D Q37A+C c "D Q31Q4AAE!Q$O PC "-AaD!9 Bc II Ir7c  td\}}ttjtdtttzz }t |dk(sJt tt|tdzdz dzz|tztdzdz dzzzt tdzz}t|d |k(sJt||d k(sJt tt|tdzd zdz tdzdzdz ztdzdzz dzz|ttdzd zdz ztdzdzdz z tdzzzt tdzz}t|dd |k(sJt tt|tz|zt tdzz}t|dd|k(sJt||d k(sJdtdzzttjtdzdtzttjtzzdttzz }t |dk(sJt tt|td z dz tdzzdzz|tzzt tdzz}t|d |k(sJt||d k(sJttjtdtttjtzzttz}t |dk(sJt tt|td zdz tdzdz z dzz|tztdz dz dzzzt tdzz}t||k(sJttjtdttjtztttzz }t |dk(sJt tt|td z dz tdzdz zdzz|tztdzdz tdzdz ztdz z dzzzt tdzz}t||k(sJttjtdtttzz}t |dk(sJt tt|tdzdz tdzdz z dzz|tztdz dz dzzzt tdzz}t|dd|k(sJt||d k(sJy)NzC1 C2ru)2nd_linear_airyrrrr?rrr<Trrr)r=x0)r=r)rm2nd_hypergeometric2nd_hypergeometric_Integralr)rmrr)rrr1) rr0rWrrr rrr r )r[rr\rs r8test_2nd_power_series_ordinaryrSs W FB 11a1QqT6 !B  O OO O QqT2q!tAvz?RT1a47Q;%77!AqD'A BC "6 73 >> > r3 9 ,, , QqT2Az!|q1uqjl2a!eaZ?!CD a1q51*R-1q51*Q,.15 67 AqD' C "62 ># EE E QqT2a4"9qAw& 'C "6! < CC C r3 9 ,, , ad(QqTYYq!_ %!QqTYYq\(: :AadF BB   ! !! ! QqT21uQwA~)*RT1AadG; > > r3 9 ,, , 11a1adiil+ +ad 2B  K KK K QqT2q!tAv1Q*+bdQTE!GaK.@@1QT7J KC ":   11a1Q499Q< '!AaD& 0B  = == = QqT21uRx!Q$q&(1,-dAqDGad1f$qs*Q./023AqD': ;C ":   11a1QqT6 !B  O OO O QqT2q!tCx!Q$q&(1,-1q!teBhl0CCa1gM NC "6! < CC C r3 9 ,, ,r7c td\}}}tdzttjtdzdtzttjtzz dtzdzttzz}t tt|tdzzdtdzzdz dtdzzzdtzz dzzt tdzz}t |d |k(sJt||d k(sJdtdzzttjtdzd tdzzttjtzz dtdzzdzttzz}t tt|ttztdzd z tdzdz ztdzdz ztzdzzt tdzz}t |d |k(sJt||d k(sJtdzttjtdztdzttjtzz tdzdz ttzz}t tt|tdz dz dtdzzdz z tdzd z z tdzdz ztdz zdzztz |tdzztdz dz tdzdz ztdz zdzzzt tdzz}t ||k(sJt||d k(sJtdzttjtdztttjtzztdztddz ttzz}t tt|tdzd z tdzdz z dzzttz |ttztdzdz tdzdz z dzzzt tdzz}t |d |k(sJt||d k(sJtttjtdzttjtz|tzttzz }t tt||dztdzzdz |tdzzdz zdzzt tdzz}t |ttd |k(sJt||d k(sJttjtddtz tz ttjtzz|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 |k(sJt||d k(sJy)NzC1 C2 arurrir r?r2nd_power_series_regularr<rNrrirP<rJ@i@8i@$) rrWr0rr rrr rr )r[rryr\rs r8test_2nd_power_series_regularr[s "IBA AqtyyA !A#qtyy|"4 4!a1~ EB QqT2ad7C1HQJ1a4/!A#59:Qq!tWD EC "5 6# == = r3 9 ,, , 1a411a !1QT61Q499Q<#8 8AadF = Q4< B QqT2d1g:q!tBwAa/!Q$q&811DqHI!K q!tGadU2X1R'!A#-1 2356q!tW= >C ":   r3 9 ,, , AqtyyA !QqTYYq\"2 2adXa^6KQqT5Q QB QqT2q!tBwAa'!+,T!W4 47 AqDHq!tAv%)*+-.q!tW5 6C "5 6# == = r3 9 ,, , 1Q499Q? QqTYYq\ )AaC!H 4B QqT2q!tAqDy|a1fQh.23a1g= >C "ad!; < CC C r3 9 ,, , 11aQUAIqtyy|3 3qsAaDj @B QqT2r!Q$wAA.A6A>uD !Q$AAA&s*+-.q!tVQU^QU-CB-FG !Q$AqQ3!"#$&'1g. /C "ad!; < CC C r3 9 ,, ,r7ctd}td}ttdtd|z|k(sJy)Nrrr[r)r r$r)rrs r8test_issue_15056r]s;s A B !&,$1B"B Cr II Ir7c t tz dtzttzz ttz ttttz}tt tdztzzt tdztzt tt tdz tz ztz tzz}t||dk(sJttt t tzzz}ttttztttjtdz}t||ttdk(sJy)NrurN) r[rWr0rrrrr rrr\rs r8test_issue_15913r`s Q1QqT AaD :adA#6 6B SA] SA]8BsAqD519~4Ea4G+KK KC r3 9 ,, , r#qbd)| C AadFA 11a 0B r3! % 22 2r7cHttdttdy)Nctttjtt tjtgttt tt tgSr4)rr0rWrr1r2r6r7r8rrz"test_issue_16146..s@vqtyy|QqTYYq\&BQqT1Q4QRSTQUDVWr7ctttjtt tjtgttgSr4)rr0rWrr1r6r7r8rrz"test_issue_16146..s0vqtyy|QqTYYq\&BQqTFKr7)r*r5r6r7r8test_issue_16146rds :WX :KLr7cRttdz ttjtz}tttt}t ||k(sJttt tz ttjtdz}ttttttzztttt th}tt ||k(sJttdzdttzz dzttjtdz}ttttttzzttdzzz}t ||k(sJy)Nrur?r) r0rWrr r[rrrrrr_s r8&test_dsolve_remove_redundant_solutionsrfs A$q&!A$))A, B QqT2,C ":   A$s1v+! !Q (B adBAI 1Q4Q 0 1C vbz?c !! ! A$'!AaD&. AaDIIaO +B QqT219r!Q$w& 'C ":  r7c.tdt\}}td}tt |||||||ztt |||||||zg}t |}t ||dddgfk(sJy)NzA BrrTr)rrr r rrr )rBrr\rs r8test_issue_13060ris 5h 'DAqs A Z!a !A$qt) ,bAaD!1Dad1Q4i.P QB *C r3 D1a&> 11 1r7ctd\}}td}d|zt|dz z||zd|dzz|dz z|d|z|dz zz dzzt|||dfzz}t |dd }|d d d k(sJtd \}}t |d }t |t tdk(sJy)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 rkeqnmatchr[rrs r8test_issue_22523rqs 5>DAq 5/C a%QU CF "a1fa!enAaCQK!7K&K SVaV $&% %C 4e 4E , -g 6! ;; ; W FB 6 7C c?j.K*LM MM Mr7cRtdt\}}tdd\}}}}}d\}}}}t|t|||dz|||zz|||||z zz d}t|t|||dz|||||z zzd}||g} t | ||||g|dd||jj |dd|dd ||jj |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 zk(sJ| 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 zk(sJy)Nzx1, x2rz t k1 k2 m1 m2Tr.)r?r?r?r?rurr?rdrr rX)rrr rrrrr) rx2rrrrrr2r3eqsx1solx2sols r8test_issue_22604ry-s XX .FB=Ar2r2NBB RRUAq! !Br!uH ,r2a52a5=/A A1 EC RRUAq! !B11 $6 6 :C *CC"Q%A1aAARARSTUVAWXY8:1aAARARSTUVAWXY8[\NUE Br!ud1tAw;/bAd1gI1EFs4PQ7ST9UYZ[^bcd^eZeUfKfghKhGiijll47 T"X-tDGaK/@@T!WQYtTXYZT[^_T_O`E`abEbAccdffgh hh h Br!utAw{CQ $q47{:K0KA0M,NNrQUVY]^_Y`U`befjklfmnofoptuyz{u|@AvAqBgBCDgDcEUEFHUHHI II Ir7c tttjtdttdzzdttzdz z tttjtdttdzzdttzdz z fD]}dt |ttvrJy) Niruii rrr1rC)r r0rWrr)des r8test_issue_22462r|;s qtyy|S1q[3qt8D=8 9 qtyy|R!aZ!AaD&(2 355l2qt44445r7ctd}td}tt|z ||j j z||j zd}t |dk(sJy)NrWrr)rtnth_order_reduciblerr)rrr rrrrWrr\s r8test_issue_23425rBsi A A QTE!A$))+""$ $qtyy{ 2A 6B   ; ;; ;r7ztoo slow for @slowctd}td}||dz||j|dzdz}t|||Jy)NrWrrru1)r rrrrs r8test_issue_25820rKsPs A A 1q11a 2 %B "ad  '' '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*r/r7rErHrKrSr[r]r`rdrfrirqryr|rrr6r7r8rsBB33$"/9=79KK.& *4HHN7+::+/6g.>+BBBBB W4 ( 1a SM SM SM% #3L8 D?D?N\~i$V,*^Y?vF47B//*lM^=,;9k I7 G 3$7@F J(/-/-d$-NJ 3M 2M, I5;((r7