wgg *dZddlmZmZddlmZddlmZmZm Z m Z ddl m Z m Z ddlmZddlmZmZddlmZmZdd lmZmZmZdd lmZmZmZmZdd lm Z m!Z!dd l"m#Z#dd l$m%Z%m&Z&m'Z'm(Z(m)Z)m*Z*m+Z+ddl,m-Z-m.Z.ddl/m0Z0ddl1m2Z2m3Z3ddl4m5Z5ddl6m7Z7m8Z8ddl9m:Z:m;Z;mZ>ddl?m@Z@ddlAmBZBmCZCddlDmEZEmFZFddlGmHZHmIZImJZJmKZKmLZLddlMmNZNddlOmPZPmQZQddlRZRe8dZSe8dZTedZUe8dZVe7dZWe7dZXed \ ZYZZZ[Z\Z]Z^Z_Z`ZaZbed!\ZcZdZed"Zfd#Zgd$Zhd%Zid&Zjd'Zkd(Zldqd)Zmdqd*Znd+Zod,Zpdqd-Zqgdfd.Zrd/Zsd0ZteQd1ZueQd2ZveQd3ZweQd4Zxd5ZyeQd6ZzeQd7Z{eQd8Z|eQd9Z}d:Z~eQd;Zd<ZeQd=ZeQd>ZeQd?Zd@ZeQdAZeQdBZeQdCZdDZeQdEZeQdFZeQdGZeQdHZdIZeQdJZeQdKZeQdLZdMZeQdNZdOZdPZdQZeQdRZdSZdTZeQdUZeldVZeldWZeldXZeldYZeldZZeld[Zeld\Zeld]Zeld^Zeld_Zeld`ZeldaZeldbZeldcZelddZeldeZeldfZeldgZeldhZeldiZeldjZeldkZeldlZeldmZeldnZeldoZdpZy)ra This File contains test functions for the individual hints used for solving ODEs. Examples of each solver will be returned by _get_examples_ode_sol_name_of_solver. Examples should have a key 'XFAIL' which stores the list of hints if they are expected to fail for that hint. Functions that are for internal use: 1) _ode_solver_test(ode_examples) - It takes a dictionary of examples returned by _get_examples method and tests them with their respective hints. 2) _test_particular_example(our_hint, example_name) - It tests the ODE example corresponding to the hint provided. 3) _test_all_hints(runxfail=False) - It is used to test all the examples with all the hints currently implemented. It calls _test_all_examples_for_one_hint() which outputs whether the given hint functions properly if it classifies the ODE example. If runxfail flag is set to True then it will only test the examples which are expected to fail. Everytime the ODE of a particular solver is added, _test_all_hints() is to be executed to find the possible failures of different solver hints. 4) _test_all_examples_for_one_hint(our_hint, all_examples) - It takes hint as argument and checks this hint against all the ODE examples and gives output as the number of ODEs matched, number of ODEs which were solved correctly, list of ODEs which gives incorrect solution and list of ODEs which raises exception. ) Derivativediff)Mul)EIRationalpi)EqNe)S)Dummysymbols)imre)LambertWexplog)asinhcoshsinhtanh)cbrtsqrt) Piecewise)acosasinatancossecsintan)Eierfi)hyper)Integral integrate)rootof)FunctionSymbol)airyaiairybibesseljbessely lowergamma)NonElementaryIntegral) classify_odedsolve)allhints_remove_redundant_solutions) FirstLinear ODEMatchErrorSingleODEProblemSingleODESolverNthOrderReducible) checkodesol)raisesslowNxuyfgzC1:11a b cz^Hint did not match the example {example}. The ODE is: {eq}. The expected hint was {our_hint}zDifferent solution found from dsolve for example {example}. The ODE is: {eq} The expected solution was {sol} What dsolve returned is: {dsolve_sol}zFsolution found is not correct for example {example}. The ODE is: {eq}asolution returned by dsolve is incorrect when using {hint}. The ODE is: {eq} The expected solution was {sol} what dsolve returned is: {dsolve_sol} You can test this with: eq = {eq} sol = dsolve(eq, hint='{hint}') print(sol) print(checkodesol(eq, sol)) zdsolve raised exception : {e} when using {hint} for the example {example} You can test this with: from sympy.solvers.ode.tests.test_single import _test_an_example _test_an_example('{hint}', example_name = '{example}') The ODE is: {eq} zTested hint was : {hint} Total of {matched} examples matched with this hint. Out of which {solve} gave correct results. Examples which gave incorrect results are {unsolve}. Examples which raised exceptions are {exceptions} cfd}|S)Nc}g}|dD]}|d|d|d|d|d|jdg|d|jd|d||d|jdd|d|jdd |d|jd d|d|jd d|d|jd d|d d }|j||S)NexampleseqsolXFAILfuncr;F simplify_flagTcheckodesol_XFAILdsolve_too_slowcheckodesol_too_slowhint) rErFrGrH example_namer;rIrJrKrLrM)getappend)solverrDexampletemprHs h/home/mcse/projects/flask/flask-venv/lib/python3.12/site-packages/sympy/solvers/ode/tests/test_single.pyinnerz _add_example_keys..inners.vj) "GZ(1$7j)'259 +G488"Ez*7377vf~N 'z*7377F &z 27 ; ? ?PT U%+J%7%@%D%DEXZ_%`"("4W"="A"ABSTY"Z'-j'9''B'F'FG]^c'dv D OOD ! ")rHrUs` rT_add_example_keysrXs& 7NrVc|D]9}|s|ds|r|dst|d|d}|ddk7s,t|d;y)Nr;rMT) solver_flag xpass_msg)_test_particular_exampleprint) ode_examples run_slow_testrRresults rT_ode_solver_testrbsU'GFOPWX^P_ )'&/7PTU + " $ &% & 'rVcttdgz}t}|D]%}|jdsd|vrt |||'y)Ndefault _Integralseries)listr2_get_all_examplesendswith_test_all_examples_for_one_hint)runxfail all_hints all_examplesour_hints rT_test_all_hintsrosOX {*I$&LJ   [ )X-A ',IJrVcrt|tk(rtfd|DSj|S)Nc3@K|]}j|ywN)dummy_eq).0sub_dsol expected_sols rT z"_test_dummy_sol..sNx<((2Ns)typerganyrs)rv dsolve_sols` rT_test_dummy_solr{s2 JN:NNN$$Z00rVcPt}|D]}|d|k(s t||y)NrN)rhr])rnrNrmrRs rT_test_an_exampler}s0$&L8 > "l 2 $Xw 78rVc|d}|d}|d}||dv}|d}ddd}|d} |d } |d } |d } d } |r1|t||vr#tj||| }t||t||vrt||d< | st ||| |}nt |dk(r|d}n|}|r4|gk7r.d}t|tk(r3|D]-}|jtrt!|| }n||v}|s-n0n.||v}|D]%}|jtst!|| }'|r$t"j||||}t|t%t |Dcgc]}d}}t |dk(rd}| s`| s^t'|||d|k7rL||d<d} t(j||||}|r"t*j||}t|d|z|d<| r |r |dz|z|d<|S#t $rb}g}||d<|stjtjt|||||d<|r|st|dd} Yd}~d}~wwxYwcc}w)NrErFrNrGrHr\)msgr[rIrJrKrLT)rRrErn match_listsimplifyrMrexception_list)erMrRrErF)rRrErFrzTr)solve_for_func unsolve_list)rMrErFrz)rRrEzAssertionError: zis now passing for the hintr[)r0 hint_messageformatAssertionErrorr1len Exception traceback print_exc exception_msgstrr^rxrghasr r{expected_sol_messageranger9dsol_incorrect_msgcheckodesol_msg)rn ode_examplerZrErvrRxfailrHrarIrJrKrLxpassmessagerzrexpect_sol_checksub_soliexpected_checkodesols rTr]r]s T Bu%L.)G  G, ,E v Db )Fo.M#$78!"34O&'=> E <D1 1"))'b8)TG ) )<D))&| ##B}(S |$a'!-aJ!-J :r>$ J%+G{{5)/>w /S+S(+2*+D(' $.\#A +TG{{5)/>w /S+S(T .55g"R^ku5v$W--7,HF5M U"),I"IH"TF;  Mc J'.F# $##%)003q6RY^`0aF5M5fUm$E >$Qs?(G/! I/ I8AIIc|gk(r t}ggg}}}|D]x}||dv}|r|s|rt||}||jdgz }||jdgz }||jdgz }|`|d} | dk7skt|dz|Tt |} t |t |z t |z } t j || | ||} t| yy)NrGrrrrr\)rMmatchedsolveunsolve exceptions)rhr]rOr^rcheck_hint_msgr) rnrmrkrrrrrrar match_countsolveds rTrjrj'sr(* /12rn J# % K00 E   )(K@fjjb11  >"55 &**%5b99  -CBwfUm$ %*o Z\!223~3FF##;f^jxF#G c rVc~tttjtttt}t |t t fdt t fdt t fdt t fdtttjtttttzzttt}t|t tfdt|jdusJtttjtttzttt}|jdk(sJtttjtdttjtd zttjtd z ttt}|jdk(sJtttjtd ttjtd zttd zz ttt}|jd k(sJtttjtd tttjtd zzttd zz ttt}|jdk(sJy) Nc$jSrr)matchesrQsrTz&test_SingleODESolver..Fs(8rVc$jSrrget_general_solutionrsrTrz&test_SingleODESolver..Gs(C(C(ErVc$jSrr)_matchesrsrTrz&test_SingleODESolver..Hs(9rVc$jSrr)_get_general_solutionrsrTrz&test_SingleODESolver..Is(D(D(FrVc$jSrrrrsrTrz&test_SingleODESolver..Rs&"="="?rVFrT) r6r?r<rr7r:NotImplementedErrorr4r5rorder is_autonomous)problemrQs @rTtest_SingleODESolverrAs qtyy|QqT15G W %F  89  EF  9:  FG qtyy|ad1Q4i71qAG  !F =?@  !F >> u $$ $qtyy|ad2AaD!11QOQRSTQUWXYG ==A  qtyyA11a@1Q47JAaDRSTG  D (( (qtyyA1Q499Q?1BBQqT1WLaPQdTUVG  E )) )rVc"ttyrr)rb)_get_examples_ode_sol_linear_coefficientsrWrVrTtest_linear_coefficientsrfs >?rVcttjttttz z }t|tt}tttjtzttz }t|tt}d}d}||vr||vsJ||vr||vsJt t t t y)N(1st_homogeneous_coeff_subs_dep_div_indep(1st_homogeneous_coeff_subs_indep_div_dep)r?r<rr0rb>_get_examples_ode_sol_1st_homogeneous_coeff_subs_dep_div_indep0_get_examples_ode_sol_1st_homogeneous_coeff_best)eq1c1eq2c2sdisids rTtest_1st_homogeneous_coeff_oderjs A$))A,1a C c1Q4 B AaDIIaL.1Q4 C c1Q4 B 4C 4C b=S]* * b=S]* *STEFrVcHttdttdyNT)r`)rbrrrWrVrT,test_slow_examples_1st_homogeneous_coeff_oderysScghEUYZrVc"ttyrrrb;_get_examples_ode_sol_nth_linear_constant_coeff_homogeneousrWrVrT*test_nth_linear_constant_coeff_homogeneousrs PQrVc&ttdyrrrWrVrT8test_slow_examples_nth_linear_constant_coeff_homogeneousrsP`derVc"ttyrr)rb%_get_examples_ode_sol_2nd_linear_airyrWrVrTtest_Airy_equationrs :;rVc"ttyrr)rb_get_examples_ode_sol_lie_grouprWrVrTtest_lie_groupr 45rVc8ttjt}tttz |zttdzttztdzttzdz z z}t |dk(sJt t y)Nrr) factorableseparable_reduced lie_groupseparable_reduced_Integral)r?r<rr!r0rb'_get_examples_ode_sol_separable_reduced)dfrEs rTtest_separable_reducedrst 11B ad(B#ad1Q4i1a4!9q=9: :B   & && &<=rVc&ttdyr)rbrrWrVrT$test_slow_examples_separable_reducedrs  >1 EB QqTBHS1qs7QAY2G)K%La A#ON&8Aq1uIq1!34558(1a G D ;::61Q462,a01 2C &"?@ @@ @ r3 9 ,, ,rVc"ttyrr)rb8_get_examples_ode_sol_2nd_nonlinear_autonomous_conservedrWrVrT'test_2nd_nonlinear_autonomous_conservedrs MNrVcttjtdtttz}t t dt tdt tttzz z tttfttzt t dt tdt tttzz z tttfttz g}t|dd}t||D]\}}|j|rJt||Dcgc]}|jc}ddgk(sJycc}w)Nrr+2nd_nonlinear_autonomous_conserved_IntegralFrMrr)r?r<rrr r%rr_urr1ziprsr9doit)rEactualrass rT0test_2nd_nonlinear_autonomous_conserved_integralrs  11a4!: %B!Dab2(>&>!>??"adLbSTfUx$rAhtBx44455AaDzBBFKMF BJUZ [F66"!zz!}} rf5AFFH5 69i:P PP P5sE, c"ttyrr)rb'_get_examples_ode_sol_2nd_linear_besselrWrVrTtest_2nd_linear_bessel_equationrs <=rVc ttttttjtzz}tttt ttttt t t tzzg}t||dt}|tttt t t tzzgk(sJttyNr) r?r<rr rrrr3rb#_get_examples_ode_sol_nth_algebraic)eqnsolns solns_finals rTtest_nth_algebraicrs A$1adiil" "C !c!f  !bRTl # %E.sE1a@K 2adBs2a4yL12 22 289rVc&ttdyrrb2_get_examples_ode_sol_nth_linear_var_of_parametersrWrVrT>test_slow_examples_nth_linear_constant_coeff_var_of_parametersr sGW[\rVc"ttyrrrrWrVrT0test_nth_linear_constant_coeff_var_of_parametersr s GHrVcd}ttjtddttjtdzzttjtzdtzz tttzz }t |tt|d}t |tt|d}||k7sJt ||ddd k(sJt ||ddd k(sJy) N:nth_linear_constant_coeff_variation_of_parameters_IntegralrrTrFrrr)r?r<rrrr1r9)rnrEsol_simp sol_nsimps rT@test_nth_linear_constant_coeff_variation_of_parameters__integralrs LH 11a1QqTYYq!_, ,qtyy| ;ac AC!H LBb!A$X=Hr1Q4h?I y  r81U Cy PP P r81U Cy PP PrVc&ttdyrrb_get_examples_ode_sol_1st_exactrWrVrTtest_slow_examples_1st_exactr 4DIrVc"ttyrrrrWrVrTtest_1st_exactrrrVc:tttttttzttdzz ttj tzz }t |ttdd}t ||ddsJy)NrF1st_exact_Integralrrr)rr?r<r rr1r9)rEsol_1s rTtest_1st_exact_integralrsg QqTaAaD kAaD!G+QqTYYq\9 9B 2qte2F GE r5% @@ @rVc&ttdyr)rb)_get_examples_ode_sol_nth_order_reduciblerWrVrT&test_slow_examples_nth_order_reducibler!s >dSrVc&ttdyr)rb:_get_examples_ode_sol_nth_linear_undetermined_coefficientsrWrVrTFtest_slow_examples_nth_linear_constant_coeff_undetermined_coefficientsr$sO_cdrVc&ttdyrrb_get_examples_ode_sol_separablerWrVrTtest_slow_examples_separabler(rrVcttttttttzt j ztz }d}|t|vsJtty)N3nth_linear_constant_coeff_undetermined_coefficients) r rr?r<rr Halfr0rbr#)rErns rT8test_nth_linear_constant_coeff_undetermined_coefficientsr,sP D1qM1QqT6AFF?Q. /BDH |B' '' 'OPrVc d}t}||tttztt|tttzdk(sJ||tttztt|tttzdk(sJ|tt|tttz|tttdzdk(sJ||tttztd|ttzttztdzdk(sJ||tttd|tttdz|tttdzdk(sJ||tttd|tttdzdk(sJt t y)Ncltt|tttj Srr)r8r6r?r<r)rEs rTrz*test_nth_order_reducible..s#$%5b!A$%BCLLNrVFrrrT)rr>r?r<r=rbr )FDs rTtest_nth_order_reducibler1s]NAA Qq1vq!_q1qz ) *e 33 3 Qq1vq!_q1qz ) *e 33 3 QqT!AaD!*_q1q!} , - 55 5 Qq1vq!_q1QqT1a0 0 1U :: : QqtQ]QqtQ] *QqtQ] : ;u DD D QqtQ]QqtQ] * +t 33 3>?rVc"ttyrrr&rWrVrTtest_separabler3rrVc tttdtztz tt tttdtzzt k(sJt tyr)r&rr?r<r r%rb _get_examples_ode_sol_factorablerWrVrTtest_factorabler6 sN d1QqS6"9o%q )htB1Q3K7H!.L-L LL L56rVc&ttdyr)rbr5rWrVrTtest_slow_examples_factorabler8&s 5TJrVc"ttyrr)rb_get_examples_ode_sol_riccatirWrVrTtest_Riccati_special_minus2r;+s 23rVc"ttyrr)rb*_get_examples_ode_sol_1st_rational_riccatirWrVrTtest_1st_rational_riccatir>/s ?@rVc"ttyrr)rb_get_examples_ode_sol_bernoullirWrVrTtest_BernoullirA4s 45rVc"ttyrr)rb _get_examples_ode_sol_1st_linearrWrVrTtest_1st_linearrD8s 56rVc"ttyrr)rb#_get_examples_ode_sol_almost_linearrWrVrTtest_almost_linearrG<s 78rVc d}ttttttz tttttttzdz zttttdzdz z tt}ttttttz tttttdz ztttttdzzdz z tt}||vsJ||vsJ|dz|vsJ|dz|vsJt t y)N Liouvillerre)r0rr?r<rb_get_examples_ode_sol_liouville)rMnot_Liouville1not_Liouville2s rTtest_Liouville_ODErM@s D!$qtQ-/AaDadAq9I4I!4K"K QqT1 q#aD"N!$qtQ-/D1q!4DQ4F"F $qtQ- 1#d$N ~ %% % ~ %% % + ^ 33 3 + ^ 33 345rVctd\}}}}}td}d}tt||d|dzzdtt||dz|dzzz dt|zz}|t |vsJ|||z||zt|||zz||dzzt|||dzz}|t |vsJt t y)Nz x t a b cr>nth_linear_euler_eq_homogeneousr r$)rr(rr?r0rb'_get_examples_ode_sol_euler_homogeneous)r<trbcr>rnrEs rT*test_nth_order_linear_euler_eq_homogeneousrVOsK(MAq!Q A0H adAq !Q$ D1q!$4!4QT!9 9BqtG CB |B' '' ' 1Q4!A#d1Q4m# #a1fT!A$1-=&= =B |B' '' '<=rVctd\}}tdd\}}}}d}|dztt||dzd|dzztt||dzz d t|zz|z}|t|t|vsJ||dzztt||dz||ztt||zz|t|zz|t |zz}|t|t|vsJt t y) Nzx tza b c dTintegerDAqD1JAq!QMH Ad1Q4A AqDadAq)9!9 9BqtG Ca GB |B!- -- - 1a4QqT1a 1Q3tAaD!}#4 4q1v =#a& HB |B!- -- -CDrVc rtd\}}tdd\}}}}d}t|dztt||dzd|ztt||zz dt|zz|dz}|t |t|vsJt||d zztt||d z||dzztt||dzz||ztt||zz|t|zz|t |z}|t |t|vsJt ty) Nzx, tz a, b, c, dTrX:nth_linear_euler_eq_nonhomogeneous_variation_of_parametersr r)rr rr?r0rrb$_get_examples_ode_sol_euler_var_parar\s rTEtest_nth_order_linear_euler_eq_nonhomogeneous_variation_of_parametersrdks% 6?DAqt4JAq!QKH AqDad1Q !A#d1Q4l"2 2R!W !AadF4!Qq>$9 9AaCQqT! ? d4!a##aQqT1 o5!A$>1Q4ArEBq!tG+,- !fQTk!A$))Aq/)CF1HQqTYYq\,AACF1Q4KO1Q4CAK"SQ[.89:,-kldQqTYYq!_$q11~5!<1Q4A1b52ae83bBh>ArEIBqRTuHTUV 1Q4Aq)1Q4bdRT)*+[15 5hsA\6c dttttdztttttztttttzzdttttt t tzzt tdzdz zgdttdztttttzdtzttttzz dttzztdztttttt tzztddtdzzzzgdttdztttttztttttzz dttzz t ttz ttttt tdzzzt tdzdz z t tdz z tz gdttdztttttzdtzttttzzdttzz t tdzt tz tttttdzz t tdzzztddt tdzzz tdd t tdzzz tdd t tzz td d z gdttdzttttttzdtdzztttttzz d tzttttzzd ttzz t ttttttzt tdzzzttdzzztdd t tzz tddz gddtdzzttjtdzttztdtztt dtzdz zztttttttt tdz zt tt tdz zztdt tztt dtzdz zdz zzgddtdzzttjtdzttztt dtzdz zttttttztt tdz zt ttztt tdz zzdtt dtzdz zdz z dtt dtzdz zdz z gdddS)NrZrrrrrra @rrk rQr)euler_undet_01euler_undet_02euler_undet_03euler_undet_04euler_undet_05euler_undet_06euler_undet_07rz) r?r<r rrrrrr{rr rrWrVrTr[r[s SaDAd1Q4A&&4!a=8!<1Q4bQi#a&!)A+567 Ad1Q4A&&1T!A$])::QqtVCQTJ1Q4BAIAq!t(;;<=> Ad1Q4A&&4!a=81QqT6A3q6!8L1Q4"r!Q$w,Q14s1vby@!CDE Ad1Q4A&&1T!A$])::QqtVCSVQYQTUVQWEWX1Q4AqD2ad7*Xa]3q619-DDxPQRT~VYZ[V\^_V_G__bjklmobpqtuvqwbww{CDEGJ{KKLM Ad1Q4Aq))AadF4!a3C,CCac$qQRtUV-FWWZ[\]^_\`Z``befgbhi1Q4A1a4"QT'1HQN3q64IIHUWY[L\\]^1fQqTYYq!_$qt+d1Q3iC!HQJ.GG1Q4a"SQ]"2RCF1H 5E"EQPSTUPVWZ[^_`ab_b[cde[eWfHfghHh"hijk 1fQqTYYq!_$qt+c#ac(1*o=1Q4DGCAqM1BtAwJs3q6!8}4LLqQTUXYZ[\Y\U]^_U_Q`O`abObbefgjknopqrorkstukugvevwxexxyzA$( (rVcx dttttdztttttzdtzttttzz dttzztdztttttt tzztdzdz zzgdtdtdzzt ttttzdtzt tttzzdttzz tdzttztttttdzz t tzztttzdz zdttzdz z dttzdtzz zdttzdtdzzz z gdttdztttttzdtzttttzz dttzztdzttztttttt tzztttzzdttzz zgdtdztttttzdtzttttzz dttzzttz tttttzt tdzzzttdz ztddzgdtt tttttdfzttttztz zttttt ttzzttzttz gdtdzttj tdztttj tzzdttzzdtz z tttttdttzzt tdttzzzdd tzz zgdd d S) Nr`rrrrkrrarr) euler_var_01 euler_var_02 euler_var_03 euler_var_04 euler_var_05 euler_var_06rz) r?r<r rrrrrrrr"r rrWrVrTrcrcsQaDAj1q!,,qs:adA3F/FF1Q4OQRTUQUV1Q4BAI1Q./01 1a4QqT1a((1Q3tAaD!}+<? d:adAq))AaC 1Q40C,CCa!fLsSTvU1Q4A1a4#a&(2Xa^CDE 1vg:adQF33j1q6II1LL1Q4bQi#a&02a589: dQqTYYq!_$q1Q499Q<'77!ad(BQqSH1Q4C#a&M)Bs1SV8},<rrrWrVrTr@r@ s aD1Q499Q<!A$&1q0!41Q4BqD1H&'%&diilQqtV#1Q4C!H%& d1Q499Q<!#1Q4$rAaCx.)2adDacN+CD rVcdttddttjtzttdzztttz z dtdzzztttt d t t t dttzdz zzdzdtzz gdidS) NRiccati_special_minus2 riccati_01rrrrrz)r?r<rr rr!rrrWrVrTr:r:#s-aD! ! nqtQw&1a/!AG);1Q447(3rDGCFN1,<'<#==AAaCHIJ  rVcBdtttttjtttdzzdz dtttt dt t dt dztzz zt t dt dztzz z gdttdzttttzdttztz zdtdzz ztttdt ztz tt tzzz gddtdzzttttztdttzttttzdz zz ttdz ttzztttt dtdzzzt tzz gdttjtdz d ttzz ttdzz tttdt zd t tzzt t tz z gdtdzdtzdtz zttzz ttdzzttttzttttt tdzzdzzt tdzzdz z gdtdzttttztdzztdttdzzttttzzz ttzttttt tztz dzzt tdzzdz z gdttdz ttttzd tdzzd tzz d ztdz dzdtzdz dzzz ztttdt ztzdt zz d td zzz dtdzzzdtd zzz dtdzzzdtzz dzdt ztdzzdt ztzz d t zzdtdzzzdtd zzz dtdzzzdtd zzz dtdzzzd tzz dz z gdttttd tdzzdzttdzztz zdtdzztz d zttzttdz zz zd tdzzdtzz dzttdz dzzz ztttt td zz tdzzdtzz dzt tzt z tdzztd zz tdzzdtzz dzz gdttttdtzttdzdzztz z tttt t zttdzzztzt tdzzdz z gdttttttttztd dz dtzz z tdz tdd z z ttdzzdtzd z tddz z z ztd dz z dtdzzdtzz dzdtd zzdtdzzz z ztttdt ztdzzdt ztdzzzd t ztdzzzd!t ztd"zzzd#t ztd$zzzd%t ztdzzzd&t ztdzzzd't ztd zzzd(t ztdzzzd)t ztd zzzd*t ztdzzzd+t ztd zzzd,t ztdzzzd-t ztzzd.t zztt dtzzzdt dtzzz td/t ztdzzd0t ztdzzzd1t ztd"zzzd2t ztd$zzzd3t ztdzzzd%t ztdzzzd4t ztd zzzd5t ztdzzzd6t ztd zzzd7t ztdzzzd8t ztd zzzd9t ztdzzzd9t ztzzd:t zzt dtzz zz gdd; d<S)=N1st_rational_riccatirrrrrrrkrrrrm<^H!ra9:iiirrb(rfrPii| ri0H iiFi4"i<ii2ҍi"0 iimi^iHi iq!tVac\A%1q51* 561Q42#1*q!t+ac1A51r AqD8H1a48ORSUVRV8V aC991q!AaC1q1$5a$771Q41"R%!AqD&.1,rAqDy1}=>?AaD!$a!fad1fqsl&;qsQqT!V|QqTSTW>T qSUQqT!V^?' tAv'&),QTDF):S)@2ad7RPQSTPTWCT(U'VW1Q4"R%2+2ae 3c"fQUlBT"WQPRU]R "HQUN#BYq!t^,.4Ri1n=?Frz!Q$O BJq!tO&rk!Q$./19"QT1ABDMbLQRTUQUDUVl1a4 #,B,q.14=R<@CDC!H*MPQQTUVWXUXQYzZ"R%2+Bq"u ,s2vae|;d2gaemKeTVhWXZ[W[m[Ri1n%byA~.07 1a4@BI"*QPQT/R RK1 '{1a4/02:2+ad2BCEMb[QR]S RKac(#$ %&' kAE ErVcdttdtttjttttzztdztttttt tdzdz zzt dt tztt dtzdz zdz z t tdz dz zgdidS)N 1st_linear linear_01rrrz) r?r<r rrrrr r#rWrVrTrCrCs!aD11!A$&1-1Q4"qQT!V},d1gd2h.>tDGAIaK?P.PQR.RRTWYZ\]Y]X]^_X_T``ab  rVc%td}td\}}}}}dttidttttttj tzzt ttdt ttt tz ggdddttttj ttttzzd zzt ttt td ttzttd tzd z zz ttd z d z zt ttdggd dd ttj ttttd zzzttj td tttzzzt ttt tt zttt zzt ttt ttd z d z zgd dttj ttttd zzzttj td ttzzt ttt ttd z d z zt ttt ttztttzzgd dttj td zdz ttj td zdz zt ttt tz t ttt tzt ttt d tzzt ttt d tzz gd dttj td tttz ttj tzt ttt!t ttt ttzzdzz t ttt!t ttt ttz zdzz t ttt gdddttj td zdz ttttj tzdz zt ttt tz t tttt d tzz t tttt d tzzt ttt tzgd dt#tttdzd t#tttd zzz dzt ttt tz t ttt tzgd dttd zt#tttdzzd ttd zzt#tttdzzz ttd zt#tttd zzzd ttzt#tttdzzz dttzt#tttd zzzd ttzt#tttzz t#tttdzzd t#tttd zzz dzt ttt tz t tttt d tzz t tttt d tzzt ttt tzgd dtdzttd zzd tdzzttzt#tttd fzztdzt#tttd fd zzzd td zzttzt#tttzzd td zzt#tttzt#tttd fzzdtd zzttd zzz dtd zzttzt#tttd fzz td zt#tttd zzzdtzttzt#tttzz dttd zzzt ttt t%d tztt'd tzzt ttt t%td tztt'td tzzgdddttj td tttz ttj td tttzzt ttt!t tt tdt z zttzzdz z t ttt!t tt tdt z zttz zdz z t ttt!t dtt tdt z zttzzz z t ttt!t dtt tdt z zttz zz z gddd tttj tttd zz t t)dt!|d zz |ttft tzgd!gdd"ttj td zttd zz t ttdt d zd t ztzz td zzz gd!gdd#ttj td zttz t ttt d zdz t tzd z z td zdz zgd d$ttj td zttd zz t ttt ttzt ttt tt zgd d%ttj td zttd zz t ttdt d zd t ztzz td zzz gd d&ttj t|tdzz|td zzz|td zzz|tzz|zd'zz t ttt t+dt||tzz|td zzz|td zzz|tdzzzz tzgddt td tztt#ttt ttt t+tt-td tzz tzt ttt t+t-td tztzgdd(t#tttd ztd zz t ttt d tt/dd zzdz z t ttt d tt/dd zzdz zgd tttj td ztttzz t ttt tt ztttzzgd d)d*S)+z some hints are marked as xfail for examples because they missed additional algebraic solution which could be found by Factorable hint. Fact_01 raise exception for nth_linear_constant_coeff_undetermined_coefficientsr>za0, a1, a2, a3, a4rfact_01r) separable 1st_exactrr1st_homogeneous_coeff_bestrrrrZ1nth_linear_constant_coeff_variation_of_parametersr`r*rofact_02r)rrrfact_03rrfact_04fact_05rrfact_06TrErFr;fact_07fact_08fact_09rkrfact_10rrbfact_11rrErFrKfact_12rfact_13fact_14fact_15fact_16fact_17grp)fact_18fact_19fact_20rz)r rr?r<rr rrr r#rr*rr+r rrrr,r-r/r%rr)r>a0a1a2a3a4s rTr5r5s c A12NBr"R aDK dQqT!A$))A,&&1Q4 R!b1f-.? KdAaDIIaL1a')*1Q4"tAwtBx/T!WQYq[0AAA31uQw<OPRTUVWXUY[\R]^9K$tyy|AaDAI%! !Q!AaD&(@AAaD"VQBZ-"VQBZ-78AaD"S!Q$q\/9RS%K.tyy|AaDAI%! !Q!A$(>?1Q4CAaL)2adBs1vI3q6 4I+JK/K8tyy|Qq 1Q499Q<?1#451Q4a "QqT26"2BqtR!A#X4F1Q4QSVWXYVYQYHZ[9KBtyyAs1Q4y(!A$))A,6 qtS"c$s)R!V"459:; < qtS"c$s)R!V"459:; < qtRL CKVtyy|Qq 1Q4! ! #4Q#671Q4a "QqTDacN?";BqtT"qQRs(^? @ qtSSDBKa!89A=>? @ qtSQRRU ^R!V%? @ qtSQRRU ^R!V%? @    YKp!A$))A,!a'(3q!t9q!A$i@"q&IJqK|diilAo!a'1Q4BEAbDFNQT1234}KHdiilAo!$1Q4Qq2a46)AqDF234IKRdiilAo!a'1Q4CF#R!baRj%9:SK\diilAo!a'1Q4BEAbDFNQT1234]KhdiilBq!tGbAg-1a47"Q$>CtLL1Q4hqb2a4i"QT'.ABq!tG.KbQRTUQUg.U)V'VXYZZ[\iKx1q53z!A$/011Q4hrD1Q3L'8!<<=r!A$XVZ[\]^_`]`[aVbdeMfHf?gh 1q!1$q!t+1Q4a8Aa= 00223R!b1QQRST EUCUVWCW>W5XY ! !Q!AaD&(1Q4CGbQi/01OKO OrVcddlm}tdd}td}|tj t}d|ttdz|tdzz|z|td zzd zt |tttd tz zd z td d zt |td td tzz ttd tz zd z td d zzdz t |td td tzzttd tz zd z td d zzdz gd t|tz|zdtz|tdzzzd zt |tttd|d tzzz tdtzz t |tttd|d tzzz tdtzzgd t|zt|tzzd zt |tt|tz tt zgd tt|tz|zt|tzd tzzt |ttttz ttd dzz gd t|tt |ttz|tzzzt |ttz|tzdzt |ttttt |d zdf| tz|ztz d z ttz|d zz dfztt zgd ddS)Nr)r"ATpositiver? almost_linearrrrrrr) almost_lin_01 almost_lin_02 almost_lin_03 almost_lin_04 almost_lin_05rz)'sympy.functions.special.error_functionsr"r)r(r<rr rrrrrrr)r"rr?r]s rTrFrF)s:sT"A A ! ! A$aDd1Q47l1nqtQw&*1Q4"S1X+/HQN:; 1Q4"tAwqy.2c!A#h;?Xa^"CCAEF 1Q4"tAwqy.2c!A#h;?Xa^"CCAEFH!fQh1QqT1W$q(1Q4$Qr!A#wYBqD 9::;R!dBQRSUVWXYVYSZQZN\_`bcd`d\eKeFf=gh sQqtV|a1Q4"r!u*c1"g-./ AaD k!mc!A$i'!A#-1Q4RTAhq!n$44567 !QadA&1-..adA>1EI1Q4"y Bq1uaLaRTAX\A-s1v5q1u=tD FFGJA2wOPQ1! !rVc td}td}dttt ttttz t ttttdz zt tttdzdz z t ttt ttttzzz gdt tttt zttt ttt ttttzzz gdt ttttz t ttttdz zt tttdzdz z ttt ztttz jt ttt ttttzzz gdt ttttdttz t tttdzzzdtz t tttzzt tttttt tzz t tttttt tzzgdtt ttttztttz t tttdzzztt tttzzt ttttttt zz t ttttttt zzgdt ttttzj ttdt ttt tttz zgdt ttttz t ttttdz zt tttdzdz z ttt ztttz t ttt ttttzzz gdtdzt tttz|ttzttdzzt tttdzzzt tttdfzt ttttddtdzdz zztt|dzdz t|dz|zdz z|ttfzdgdd d S) Nrr>rIrrrrr) liouville_01 liouville_02 liouville_03 liouville_04 liouville_05 liouville_06 liouville_07 liouville_08rz)r)r r?r<rr rrrrexpandrr.rr/)r_ys rTrJrJTsSs A sBaD1Q4mAoQqT1a 0 22T!A$]A5Ea5GG1Q4QRT ]+,-1S!A$Z<A&1Q4QRT ]+,- adAq4!a#3A#55QqT1 q8H8JJCQRSTQUPUJVWZ[\]^[_W``hhj1Q4QRT ]+,- 1Q4A1Q4adA(:!::QqSadA=NN1Q4$rBs1vI~../AaD$rBs1vI~:N1OP QqT1a  1QqT6$qtQ-*:#::QtAaD!}_L1Q4$rBsA2wJ//0"QqT4RQBZ;P2QR !C!I+##Aq)1-1Q4R"Q$Y() AaD!}QadAq!1!!33d1Q4mQ6Fq6HH#qQRte*TUXYZ[\Y]U^^1Q4QRT ]+,- d4!Q<1QqT6AaD!G#3T!A$q\1_"DDtAaDSTVWRXGYY2:hqmQT!V<<1Q4bQA%6q%8%:!;;;bQQRTUEVWXEXYZEZA[>[[\]dii1o!A$))A,.1Q4bRTl*+, dii1o! !Q/!AaDIIaL.@1Q4bRTl*RAaC[89:dii1o! !Q/!AaDIIaO2CC !A$))A,1Q4bRTl*RAY6C!HDEF dii1o!A$))Aq/ 111Q4bdR1W,r#bd)|;<=dii1o!A$))Aq/ 111Q4bdRT!WHQJ%77"Sa^:KKLMdii1o!A$))Aq/ 111Q4bdRAaC[02c!A#h;>?@dii1o!A$))Aq/ 11AaDIIaL@1Q4bd3q6k)Bs1vI51SV CbQiORTUXYZU[R[[^`adefag^gghidii1o! ! a/1Q4d1gd2rAv;&77a@@A 1Q4d1gd2rAv;//a889; 1Q4Aq)JqtQ1,EE1Q4bQA!a5!III 23q1a4!A#:>U;; ;<=>? MLP PrVc>tt }ttjtd}dttjtdzd|zzttjtzttz tz }t d}t dt }t dd \}}}}} t d } d ttid |t|zz ttttttdz ztz ttttdzd z z dtzdz z zztt zzdz gddd||z ttttttdz ztz ttttdz z zztt zzdz gddd|dttjtzzdttzzdz tttttdtzzttt zzdzgddd|dttjtzzdttzzdttzz tttttdtzzttt zzdttzzgddd|dttjtzzdttzztttzz ttttddz tdz zttdtzzttt zzttttzzz zgddd|dttjtzzdttzzttz tttttdtzzttt zzttdz zdttzdz z gddd|dttjtzzdttzzttz tttttdtzzttt zzdttzdz zttdz zgddd|dttjtzzdttzzddttzzdttzzz tttttdtzzttt zzttzttdz zdttzdz z dzgddd|ttjtzttztdzz tttdtztdzzttttdzdz zttttdzdz zztt dz zzgddd|dttjtzz dttzz d tzttzz dtt zz tttt ttzdtt zz ttdtzzzttdtzzzgddd!|dttjtzz dtdtzzttzz ttttttdtzzzdtdtzzttzdz z tdtzttzdz z gddd"ttjtdd|zz ttztz ttzttttttdz z tttzztt zztt tzzttzzgddd#|ttjtztdzz dtzz tttttdzdz zttt zzgddd$|ttjtztz tdtzz tttttz tdtzdz z tdtzdz z tdzdz zttt zzgddd%|ttzdtzttzz tttttdzz ttzttzttzzgddd&|dttzzttdtzzz tttttdzdz z tdtzzttd'z ztdtzzzgddd(|dttjtzzttztdztt zz tttttttdzdz zzztt zgddid)ttjtdd|zzdttjtzzttzdtt zz tdztt zztttttttttdzd*z z tdz zzzzztt zgddd+|dttjtzzdttzztdtzz tdzz tttttt ztdzdz ztt#ddzz ttz tdtzzzt#d,dzgddd-|dttjtzz dttzzttt zz ttttttzttdtzzzdtzdztt zd.z zgddd/|ttjtzdttzz tz tdtzz tttt#d0d.tdz z ttd1tzzzttdz ztdtzzzgddd2|ttzttz tt z ttttttzttdz z ttzztt dz zgddd3ttjtdd|zz dttjtzzttz ttz tttttttttdz zzzzzttzgddd4|ttztj$z tdtzdz z ttttj$tdtzdz z tttzztttzzgddd5ttjtdttjtz tdtztj$tdtzdz z zz ttttttt zztttzzd6tdtzzd7tdtzzzd8ztdtzzd9z zgddd:ttjtddttjtdzzttjtzdtzz ttz ttz tttttdzzttttdz z zzttzzt tt&tdz zzzttzzgddd;|ttzttdz z tdtzdz zttttdtzd'z tttzzttdz zttzzgdddttjtdttjtz t)tz ttttttt zztttzztt)tzdz zgd?d@ttjtddAttzz t)dtzz tttttdBtzzttd,tzzzt)dtzdCz z gd?dDttjtdttjtz t)tz ttz ttttttt zztt)tzdz zttdz zttzzgd?dEttjttttztttdz zz ttttttztttzztdzttdz zdz z tttdz zdz zgd?dFttjtdttztdzttdz zz ttttttztttzztdzttdz zdz z tdzttdz zdz ztdzttdz zdz zdtdzzttdz zdz z dtzttdz zdz z gd?ttjtdttz ttdz z tttttt ztttd0zdz zttzzgd?ttjtdttzttdz dzz ttttttztttzztttdz zdz z dzgd?ttjtddttjtdzzttjtzdtzz tttzz tttttdzttttzzdz dzzztttzzttzzt t&tzzttzzgd?t||j||||z ||j|zzd||zz ||zz| t||||z| ztt| |t||dzzd|zz t|z zzd|zz zztt|| t||dzzd|zz t|z zzd|zz zzg||dGt||z||j||z||z||j|zz||z|tt| z|zzt||tt| |t||dzzd|zz t|z zzd|zz ztt|| t||dzzd|zz t|z zzd|zz zz|tt| z|zz||z| dzzt|z|z| zz dz z z g||dGtt+tttddt+tttzz ttzttttttttzzttzttdz zgd?dHdIS)JNrrrrSr=)clszR L C E_0 alphaTromegar*undet_01rrrrundet_02raundet_03rrundet_04rbundet_05rundet_06undet_07undet_08rkundet_09undet_10rmundet_11undet_12undet_13undet_14undet_15undet_16rundet_17undet_18rundet_19rundet_20rQundet_21rrjundet_22undet_23undet_24undet_25iiundet_26undet_27undet_28undet_29undet_30rundet_311irundet_32undet_33undet_34r)undet_35undet_36undet_37undet_38undet_39undet_40rz)rr<r?rrr(r)r r{rrrr r rrr|rr+r}rr) r@f2rUrSr=RLCE_0alphar's rTr#r#vs QBA 11aB !A$))Aq/AbD 1Q499Q</!A$6:A A!A!"3TBAq!S% 7OEIaDv!A#g1Q4C!HqBBAbL1Q3r64I1J,JCQRPRG+SSVWWXYv!e1Q4C!HqBB1H,=sA2w+FFJKLv1QqTYYq\>!AadF*Q.1Q4C1I3r7 2Q678v&1QqTYYq\>!AadF*RAY61Q4C1I3r7 2Qs1vX=>?'v21QqTYYq\>!AadF*S1X51Q4!A$r'AbD.2c"Q$i<"S!W*+DqQqSz+QRST3v>1QqTYYq\>!AadF*SV31Q4C1I3r7 2SVBY>3q6"LMN?vJ1QqTYYq\>!AadF*SV31Q4C1I3r7 2Qs1vXb[@3q6"9LMNKvV1QqTYYq\>!AadF*a!CF(lQs1vX.EF1Q4C1I3r7 2SV;c!fQhF3q6RSSVWWXYWvb1Q499Q<!A$&A-1Q4A13qay{+;(;bQtAwYq[AQ>Q(QSVXYWYZ[W[S\'\\]^cvn1QqTYYq\>!AadF*QqSQZ7"S!W*D1Q4!CFQsA2wY.C1I=3qs8 KLMovz1QqTYYq\>!Ac!A#hJs1v$551Q4bQqSk)Ac!A#hJs1v,=a,??#ac(3q6/RSBSSTU{vFdii1o"$qt+a/#a&81Q4SVAXbdCG(;;rBqDy#a&>PPQRGvR1Q499Q<!Q$&1,1Q4ad1fr#qb'z123Sv^1Q499Q<!#c!A#h.1Q4a#ac(1*,s1Q3x{:QT!VCbaRjPQR_vj1Q4i!A#c!f*$1Q4"q!t)SV+rAvs1vo=>?kvv1QqT6kAc!A#hJ&1Q4"q!tAv+s1Q3x/2"9c!A#h2FFGHwvB1QqTYYq\>!AaD(1a4QB<71Q4"q"q!tBw,//aR89:CvNdii1o"$q11~5!!AadF*QsA2wY61Q4CFRAaC[0AaC!GS!W3DR3GGHIivt1Q499Q<!AaD&(1,s1Q3x71Q4"b)AaC/"SAY,>"qs(CPQRSPSHATTUVuv@1Q4i#a& 3r7*1Q4CFb1Q3hA%66aRBCDAvLdii1o"$q11~5!MvX1Q4i!&& 3qs8A:-1Q4#ac(1*,r#a&y82c!f9DEFYvddii1o! ! ,s1Q3x#ac(1*9L/MM1Q4baRj2c!f94C!H r#aPQc({8RUX8XZ]^_`a^aZb7bcg7gghievrtyyA1Q499Q?!22QqTYYq\AAaCG AQ1Q4adb1b1Q3h<&7Q%??22PQRSPS8 CTVYZ[V\B\\]^ sv@1Q4i#a&("S1XaZ/1Q4QqS"r#a&y0B1Hc!f3DDEFAvLdiilQ1Q4a !MvZ1Q499Q<#ad)+1Q4baRj3sQw<>9:;[vhdii1o! ! ,tAw61Q4baRj2c!f94qay{BCDivrdii1o1Q4'$qs)31Q4C1I3qs8 3d1Q3ilBCDsv|dii1o! ! ,tAw6Q?1Q4baRj1T!W9Q;6"qs(CF9JJKL}vHdii1o!$qQU|31Q4CFRAY.Ac!a%j1BBQs1q5z\RS^STUIvRdii1o!$q!tC!H}4AaD"SV)bQi/!Q$s1q5z/"2DDq!tCPQTUPUJWXGXX[\^_[_`cdehidi`j[jkl[llopqrtuquouvyz{~zwApABCpCCFGHIFIJMNORSNSJTFTUVFVVWXSv^dii1o!$s1q5z11Q4CGrAc"gIaK'7Q&??@A dii1oad"C!HQJ/1Q4CFRAY.3q1u:a?!CDEdii1o!A$))Aq/ 11AaDIIaL@1Q3FQqSQ1Q4adAc!A#hJqL1$455bdCF8JJbSUVWSWiY\]^Y_M__`a!A$))Aa.1a4! ! #44q!A#wqt|CUK1Q41URQBD1a4!A#4FtAw4N0N,OQRSTQT,U(V%VV SQBa1fqsl+DG334ac: ; ;<=>! !A#adii!n$qs1Q499Q<'77!A$>C%PQ N@RT1Q4CAQq!tVac\(:47(B$B CQqS IJJ SQBa1fqsl+DG334ac: ; ;< c!E'!)n ac%(lQqSU5[81< =>?@! AaD!Q'!JqtQ,?*??!A$FAO1Q4"r!t)SV+c!fQh678evz zrVc*#td\}}d}d}d}||z}td}dttidttj tttz t ttt ttzgdd tttj tzttz t ttt tzgdd ttj tttzt ttt ttzgdd ttd zd ztd zd zttj tzz t tttt ttzgddttj tttz ttz d z t ttt ttz d z gddttj td tttz zd z t t ttztttzt gddtttd zzttj tzttdzz d td zzttj tzz t ttt ttdt ztztzdz zz t tzd z z d z t ttt ttdt ztztzdz zzt tzd z z d z gdddttd zd z d ttztttzzttj tzz t tttt td zzdt ztzzdt zzd z t tttt td zzdt ztzzdt zzd zgdddtttzttj tztd ttd zzzt tttt tttz gddddttd ztttztttttj tzzt tttzttttd z d z zttttd zd z z tttzt gdddttttztd ztttzttj tzz|d ztttzttj tzz t ttt!t t|d z td zzz d t"zzt ttt!t t|d z td zzzgdddttj tttttzz t ttt ttz gddtd z tttzttj tzd tztttzz t ttt"t%t td zd tzz d zztd tzzz t ttt%t td zd tzz d zztd tzzgddttj ttttttzttz z t tttt ttzgddtttj tzd ttd zztttzzt tttt tz gddddttj ttttd zzzt ttdt ttd z d z zzgdd tttd ztd zd tzzd zztttzttzttj tzzt ttttd t td zzd tzzz  t ttttd t td zzd tzzz gdttj tttzt ttt tt zgdtttd ttzztttd ttzzttj tzzt ttt"t!t ttd zz d z z t ttt!t ttd zz d z gdd tz ttj tztttd zzz t ttt tt ztz d ztd z z gdttt tttztzdtzttd zzz t tttd tt tdtd zzzd zzdz t tttdtt tdtd zzzd zzdz gdttj tttttzz t tttdt ttzz gd!gd"tttj tzd zttd zz t ttt td zz t td zzz gdt|j |d d#t&zt|zz z t t|d$tt&t d#|zz zzd$zt&z gt|d%ttj tt ttzz t ttt(tt tzzgd||||d zzz |t+||zz t ||d&t-t d'|zz z g||dd(ttj tttttz z t ttttt t ttzd z z gdd) d*S)+Nza,t`g#@g?vr separable_01r separable_02 separable_03 separable_04rr separable_05 separable_06 separable_07rrTr separable_08 separable_09)rErFr;rJ separable_10 separable_11 separable_12 separable_13 separable_14 separable_15 separable_16r separable_17rrogfffffI@gXvH ?rg Q|%Qg+R4J?)rErFrHrJ) separable_18 separable_19 separable_20 separable_21 separable_22 separable_23 separable_24 separable_25 separable_26 separable_27rz)rr(r?r<rr rrr rr!rrrrrrr rr>rrr)rSrrr@rf1r[s rTr'r'~sr  %.CAa A A A QB AaDcdiilQqT!1Q4CF#$c ! ! nqt#1Q4A cdiilSV#1Q4c!f%&c dAgkQTAXqtyy|331Q4R$q'\*+,!c*diil3q6!AaD(1,1Q4CFQ'(+c4diila#ad)m,q0A2!9s1Q4y("-.5c>d1a4i! ! $qtQw.1a4! ! 1DD1Q41"tAqtAvzA~$677"Q$(CAEF1Q41"tAqtAvzA~$677"Q$(CAEFH ?cLdAgkQqtVa!f_adiil::1Q4$r!Q$w2a/!B$6:;;< 1Q4bAg"Q&2-1235 McZAhqtyy|#d1qtQw;&771Q4b3s1v;./01! [ch!a%jQqT"S1Yqtyy|%;;1SV8c#ad)a-022SQqTQ5G5IICPQRSPTIUWYZ[ictQqT{QT#ad)^AaDIIaL881a4AaD >!A$))TU,;VV qtd2dAqD51a4<0011AbD8 9 qtT"T1a4%!Q$,//0 1  ucFdiilQqT#a&[(1Q4CF#$GcP1uc!A$i! ! ,qs3qt9}< qtR$r1a4!A#:>23qs8;<< = qtT"adQqSj1n-c!A#h67 8 Qc`diilQqT#ad)^CF221Q4RAY()acj! ! nAaD!G T!A$Z771Q4RT#$! kcxdiilQ!q\)1Q4badU1Wo-./ycB!A$'lAqD1Q3JN+q1v!}adiil.JJ qtd3q"q!t)ac/2344 5 qtT#aadQqS123 4 CcTdiilQqT!1Q4CG$% !fS1Q4[ 3q6#a!f+#5adiil#BB1Q4d2c!fai<0223R!d2c!fai<>PQR>R5ST 1uadiil"Q!q\11Q4"S!W*q.1,q1u567 d4!a= 1$qs1Q47{21Q4$q'$r#a1f+~'9"::1<= 1Q4abQq!tVnq011!346diilSQqT]*1Q4Rc!f-./011"QqT1W,1Q42#1*sQTz234diilAeAgadlN+1Q4.s1b;Ma;O6O3P/QQThhjkklm!diilR!W$1Q4C1I&' 1! ?"QAaD)9%991Q4,T"7J17L2L-MMNO!! diilS1]*1Q4c!fWbQi!m4567cg grVcdtttttttztttttzttj tzzt ttt tttz  dtzzt ttt tttz gdddtzttzdzttz tttz ttdzz ttj tzzt tttttdzz tt tt tdzzzzgdgddddtzttttzzdttzttztttz ttj tzzt ttttztttztdzzttdzztgdgddtttttttzttdzz ttj tzz t ttttzttd zd z ztgdddtzttztdzttdzzttj tzzt tdzttzttd zd z ztgdd d tttttttzttdzz ttj tzz t ttttzttd zd z ztgd d tttdzttdzzztdzttzttttdzttdzzz z ttj tzz t tttd tdttdztdzz zzttttz zdttztz dtdttdztdzz zzzz z d tdttdztdzz zztdtdttdztdzz zttztz z dttdzztdzz zzdttztz dtdttdztdzz zzzz z d ttttz zttztdttztz dtdttdztdzz zzzzz zd ttztdtdttdztdzz zttztz z dttdzztdzz zztdttztz dtdttdztdzz zzzzz zgdddt tdzttj tzd tzttzztttz z dt tttttz td zz gdtttdzz tttzdz tz ttj tzzdDcgc]B}t tt|tttdzzdzzdztz Dc}dtttzdz ttj ttdztttzz zzt ttttttdzzdttzz z t ttttttdzzdttzz zgdtdztttzttj ttztttzzt tttttt ztdzz  tzt tttttt ztdzz gdd dScc}w)a Example 7 is an exact equation that fails under the exact engine. It is caught by first order homogeneous albeit with a much contorted solution. The exact engine fails because of a poorly simplified integral of q(0,y)dy, where q is the function multiplying f'. The solutions should be Eq(sqrt(x**2+f(x)**2)**3+y**3, C1). The equation below is equivalent, but it is so complex that checkodesol fails, and takes a long time to do so. rrTrrr)rErFrGr;rJ)rErFrGr;rF)rErFr;rIrErFrIrmirA)rErFr;rKrr)rr) 1st_exact_01 1st_exact_02 1st_exact_03 1st_exact_04 1st_exact_05 1st_exact_06 1st_exact_07 1st_exact_08 1st_exact_09 1st_exact_10 1st_exact_11rz)r?r<r rrr rrr rrrrrr)rs rTrr2s aD!fS1YQAaD !1!A$))A,!>>1Q4$r#a&y/)AbD012adDCFO3LMs1Q4x!|QqT!QqTAXqtQw$6qtyy|$CC1Q4R!Q$Y1"S"q!t_2D)EEFGH cAaDQK1QqT6CF?S1Y#>! ! "LL1Q4A;QqT*QT1AaD!G;R@A !A$i1S1Y;1q0!A$))A,>>1S1Y;1q*B/0c!A$h!Q$1q.!A$))A,661a4!9qtQwqy("-. !A$i1S1Y;1q0!A$))A,>>1S1Y;1q*B/0QTAaD!G^$$1QqT 1Q4$q!tadAg~:N3N(OQRSTQUQZQZ[\Q]']]3q6 QtA!a1 $% %eAaDFm 3S1XaZ 4AaD!GAqDL !!6"# #%&tA!a1 ,<'=%= AQ1qA%&qt+A- -!A$' !Q$ >?&@ QqT!ba!A$'!Q$,.// /&1 1 %!Q-!aQqT!ba!A$'!Q$,6F1G.G!GHI  J !A$s1tA!a1 ,-ad2144q1qyA~EFF C!HQJDQqT1WQT\!122 2 3 5  567  Aadiil"QqS1X-Aq8!<1Q4"s1v+q!t+,-d1a4iAaDFQJ>1Q499Q<77?FG!1Q4!DAqD1--1145G 1vzQqTYYq\1a4!AaD&=991Q4T"q!t)aAh"6778"QqT1tBQRTUQUIXYZ]^_Z`X`L`GaCa:bc 1uc!A$i!A$))A,q.QqT"::1Q4$r#qb'z!Q$//"45r!A$RQBZPQSTPT_@U7VWSMQ QHHs>Am.cxtt }ttjtd}dttjtdzd|zzttjtzttz tz }dtt|t|zz t ttt ttdz ztz t tttdzdz z dtzdz z zztt zzdz gdd ||z t ttt ttdz ztz t tttd z z zztt zzdz gdd ttjtdz t ttt tzgdd |dttjtzzdttzzd z t ttt td tzzttt zzdzgdd |dttjtzzdttzzd ttzz t ttt td tzzttt zzdttzzgdd |dttjtzz d ttzz dtzttzz dtt zz t ttt ttzdtt zz t td tzzzttd tzzzgdd |dttjtzzttztdztt zz t ttt tttdzd z zzztt zgdd |dttjtzz dttzzttt zz t ttt ttzttdtzzzdtzdztt zdz zgdd ttjtdd|zz dttjtzzttz ttz t ttt tttt tdz zzzzzttzgdd |dttjtzzttztt tz z t ttt ttttzzztt zgdd |ttzdttz dzttz z t ttt tttdz dz ztttdzdz z ttzttttdz dz ztttdzdz z ttzzgdd ttjtd dtz z t ttt ttzzt tdzzztdztttdz zzzgdd ttjtddttjtdzzttjtzdtzz tttzz t ttt tdzzttt td z z dtttzzdz zdtt tzzdz zzzdtdtztzzdz zdtztttzzzdtztt tzzz ttzztttttzd z zdtztttzzdz zdtztt tzzdz z zzdtztdtztzzdz zdtttzzz dtt tzzz ttzzttttzzz gdttjtddttjtdzzttjtztttzz t ttt ttt td z z zzdtdtztzzdz zttzztttttzd z zzzdtztdtztzzdz zttzzttttzzz gdtttttdttzzttz t ttt tdz z tdtzdz ztdtzzttttzdtttdzzdz z dttdzzdz ztdtzzzgdd ddS)NrrrrrrrTrrarrrbrmrrkrQrr)var_of_parameters_01var_of_parameters_02var_of_parameters_03var_of_parameters_04var_of_parameters_05var_of_parameters_06var_of_parameters_07var_of_parameters_08var_of_parameters_09var_of_parameters_10var_of_parameters_11var_of_parameters_12var_of_parameters_13var_of_parameters_14var_of_parameters_15rz)rr<r?rr r{rrrr rr|rr}rr)r@rSrUs rTr r s  QBA 11aB !A$))Aq/AbD 1Q499Q</!A$6:AGaD!A#g1Q4C!HqBBAbL1Q3r64I1J,JCQRPRG+SSVWWXY!e1Q4C!HqBB1H,=sA2w+FFJKLdiilQ1Q4a !1QqTYYq\>!AadF*Q.1Q4C1I3r7 2Q6781QqTYYq\>!AadF*RAY61Q4C1I3r7 2Qs1vX=>?1QqTYYq\>!AadF*QqSQZ7"S!W*D1Q4!CFQsA2wY.C1I=3qs8 KLM1QqTYYq\>!AaD(1a4QB<71Q4"q"q!tBw,//aR89:1QqTYYq\>!AadF*QsA2wY61Q4CFRAaC[0AaC!GS!W3DR3GGHIdii1o"$q11~5!1QqTYYq\>!AaD(3r71941Q4"q"s1v+.QB7891Q4i!CF(1*SV++1Q4"s3q6A:q003s1vz?13DD a&SVaZ**SQ!_Q->>AFGHI dii1o!#1Q4bdR1W,q!tR#a&(]/CCDEdii1o!A$))Aq/ 11AaDIIaL@1Q3FQqSQ1Q4adb1b1Q3h3qs8A.E#qbQRd) TU .U+V&VYZ[^_`ab_bcd_d[eYefhYh&hklmnknorstuvsvowkw&wz{|}z}BDECEFGCGH{H'HJMNOJP&PPTVYZ\^abcdadefaf\fijklilmpqrstqtmuiuvwiw\wz{|}z}~ACDBDEFBF~GzGHIzI\IYJTJ A#c!A#a%j. TAaCjT)+,S!AY;T78;AS??AB3qs8LMNdii1o!A$))Aq/ 11AaDIIaL@3qs8K1Q4rArAaCxL01S1QZ<?BCFJJbSTVX[\]^[^_`[`V`SaNadefgdghklmnolopqlqhrdrsuduNuwz{|w}M}}ABCFGHIJGJCKAKKLM1q!$q1v-A61Q4"qs(S1XaZ/QqS9R#c!f+=M CA N 1 > Q{1}>-.1!A#h=7789 o]a arVcdtttdzttjtdztttjtzztdzdz ttzzttttt dtzt tdtzzgdtdzttjtdztttjtzztdzdzttzzttttt dtztzt tdtztzzgdtdzttjtdztttjtzztdzttzzttttt dtzt tdtzzgdtdzttjtdztttjtzzdtdzztd d z z ttzzttttt td d z d tzzt ttd d z d tzzzgdtdzttjtdztttjtzztdzdz ttzzttttt d tdzdz zt td tdzdz zzgdtdzttjtdzdtzttjtzztdzdz ttzzttttt td dz tdzdz zt ttd dz tdzdz zzttz gdtdzttjtdztttjtzztdztd dz z ttzzttttt td dz tzt ttd dz tzzgdtdzttjtdzd tzttjtzz dtzdzttzzttttdztt ddttzzt tddttzzzzgdtttjtdzttjtz dtd zzttzztttttt td dz tdzzt ttd dz tdzzzzgdtdz dzttjtdztdz ttjtzz dtdz dzzttzzttttdz tt d dtzdz zt td dtzdz zzzgdttjttdtz ttjtzzttzttttt td dz tzt ttd dz tzzttz gdtdzttjtdztttjtzztdztdzztdzz tdzz ttzzttttt ttdztttdztdzz zzt tttdztttdztdzz zzzgdd dS)N2nd_linear_besselrrrrrQrrmr) 2nd_lin_bessel_012nd_lin_bessel_022nd_lin_bessel_032nd_lin_bessel_042nd_lin_bessel_052nd_lin_bessel_062nd_lin_bessel_072nd_lin_bessel_082nd_lin_bessel_092nd_lin_bessel_102nd_lin_bessel_112nd_lin_bessel_12rz)r?r<rr rr,rr-rr rrrUrTrWrVrTrrs(aDdAaDIIaO$q!A$))A,'771a4!8QqT/I1Q4GAqM)Bwq!},<<=> dAaDIIaO$q!A$))A,'771a48QqT/I1Q4GAaCO+b1a.@@AB dAaDIIaO$q!A$))A,'771a41+E1Q4GAqM)Bwq!},<<=> dAaDIIaO$q!A$))A,'772ad7AaDF?AaD:PP1Q4GAaDFAaC002gad1fac6J3JJKL dAaDIIaO$q!A$))A,'771a4!8QqT/I1Q4GAq!tAv..GAq!tAv4F1FFGH dAaDIIaO$qsAaDIIaL'99QTAXqtOK1Q4"WT"XaZAa882gd2hqjRSUVRVWXRX>Y;YY[_`a[bbcd dAaDIIaO$q!A$))A,'771a4!A$q&=!A$:NN1Q4GAaDFA..GAaDFA4F1FFGH dAaDIIaO$qsAaDIIaL'99QqS1WadNJ1Q4Ar'!QtAwY"77"WQ$q' =R:RRSTU 11a!AaDIIaL01QT6!A$;>1Q4BwqtAvq!t44r'!A$q&!Q$:O7OOPQR sQh! !Q(AaC11+==1Q3( 1Q4O1Q4!a%"WQ!a%8"82ga1q>Q;Q"QRSTdii1o!AaDIIaL 001Q471Q4"WQqT!VQ//"WQqT!VQ5G2GGaPQR d1Q499Q?"Qqtyy|^3q!tAqDyA~17LaPQd6RR1Q4GDAJ$q!tAqDy/0ABBRPTUVXYUYPZ\]^bcdfgcghiklhlcl^m\mHnEnnopq<@ @rVc dttttdz zttjtdztddz dtzz ttjtzzdttzzt ttt ttddz zzt tddz tddz ftddz ftztt dtd dz ftzzgdttdz zttjtdztddz tzttjtzzttzt ttt dtz tddz zzt tddz dftddz fdtz ztt td dz d ftd dz fdtz zztdz tddz zz gdttdz zttjtdztdtddz tzzttjtzzttzt ttt dtz td dz zzt tddz dftd dz fdtz ztt td dz d ftd  dz fdtz zztdz td dz zz gdttddz z dttd dz zzd z dttddz zzdz z tdtzdz zzttzdttddz z tzdzzttddz ztzdzzz ttttdfzt ttttddz zt ttddz zzt tddz td dz ftd dz fttddz zztt tddz td  dz ftddz fttddz zzzzttddz zdz tddz zz gddddS)N2nd_hypergeometricrrrrr)rrrrrrPrmi`iQ rIi*rk-brrfTTrp)2nd_2F1_hyper_012nd_2F1_hyper_022nd_2F1_hyper_032nd_2F1_hyper_04rz) r?r<rr r rr$rrrWrVrTrrAs)aD1gadii1o%1a1adiil(BBQqtVK1Q4A!QKqtAvqtAv.>1a 1(MMPRSXYaefgheidijkdkcmopSqPqqrs 1gadii1o%1a1Q499Q<(??!A$F1Q4"a!eqtAv..uad1fa[1Q46)QQRU/SS UQqTE!GR=AaD57*a!e 4 45781u!Q6GHIJ1gadii1o%1qtAvax11(EE!L1Q4"a!equQw//qtAvqkAbE!G:qSTu0UU UQqTE!GR=AaD57*a!e 4 45781u"a6HIJK1Q46{lDQqT!V,Q.a!A$q&k1A"1DDqvaxPQzQRSTUVSWWY\_`cdefcghici_j^j _ ^ZAaDF a!#Z$%'1!A$A'?@1Q4QrU2X1qtBw<qtAv!uQw>ORSTURVWXRXQZ\]`abc`def`f\g8h(h UAaDGaeVBY'!A$q&A!QK @ @)ABCDqtAv;QR?VWXZV[\^V^B_`ab % rVc\ dttttjtdtttzt ttzt t dttdtzt tzz dtzzdttzz z tttfttzt t dttdtzt tzz dtzzdttzz z tttfttz gddttjtdtttzdttz zt tdt dtdtzdttddzzz dt tzz z tttfzttzt tdt dtdtzdttddzzz dt tzz z tttfzttz gddttjtdtttzt t dttdttzzz tttfttzt t dttdttzzz tttfttz gddttjtdtttzt t dttdt!tzz z tttfttzt t dttdt!tzz z tttfttz gddttjtdt#ttzt t dttdtzt#tzz dtdtdzz zz z tttfttzt t dttdtzt#tzz dtdtdzz zz z tttfttz gddgd d d S) N"2nd_nonlinear_autonomous_conservedrrFryrrrrErFrIrG)%2nd_nonlinear_autonomous_conserved_01%2nd_nonlinear_autonomous_conserved_02%2nd_nonlinear_autonomous_conserved_03%2nd_nonlinear_autonomous_conserved_04%2nd_nonlinear_autonomous_conserved_05rz)r?r<rrrr r%rrrrrrr rrrrrWrVrTrrbs9aDdii1oAaD )C!I5 x$rAbDRL01R47!CG)CDDr1Q4jQSUXYSY Z x$rAbDRL01R47!CG)CDDr1Q4jQSUXYSY Z  .dii1oQqT *QqtV3 tAwx$qtaHQN0B.B'BQs2wY'N"O ORTVWXYVZQ[\\^`cd^d e tAwx$qtaHQN0B.B'BQs2wY'N"O ORTVWXYVZQ[\\^`cd^d e  .dii1oAaD ) x$rAc"gI~..QqT ;R!V D x$rAc"gI~..QqT ;R!V D  .dii1oQqT * x$rAd2hJ//"ad!#c!Q$qt)hq!n*D&Ea&GGcRSfUV diiladad1Q4i!m453q!tAaDy>!#c!Q$qt)hqm*C&DR&GGcRSfUV dQqT1a4!9%%3q!tAaDy>!#c!Q$qt)a-&8&::BQKHI11!Q!q1Q46{1Q4/?!.C*C DDaH2c!ad1f+ad*++aAa!A$q&kM!A$'4IA4M0N.Nq.PPRTWZ[\W]R]^_11!Q!q1vk/ ::A>1Q4$q'$q"s1v+"771=>adDGDBQK$991Q3?@B11AadAg1QqT !91Q4 ?AqD1qLSTVWSWXYZ[X\^_X_S_D_A` aacde1Q4a1a4j)* 111q!A$ 22A6 qtd1gXd1b3q6k?33A5 6 qtT!WT!R#a&[/2214 5 111a1 q!A$q&z ::A>C!qM>!#aAaD k!m3R#a&[AB 111a1 a/31Q4Bq!tGaK)* 111q1qt 3Q 77;SV8c!A$(m+c!A$i7!AaD&@"EFdiiladad1Q4i!m451Q4$q'$qHQqM)9'943r"u:PQTWPXAXZ]^`ac^cZd@d;ehklnoqlqhrstvwswhw;w{CDEFG{H;H(HAx!}SBZqRU|a#h'>BrE &J!KcRTUWRWjYZ\]Y]o!]`hijkl`m mm(npqrsuvrvpv(w#xxxyzq'$r!Xa]""D!C2J,S*A3r"u:)M$NQTUWXZUZQ[\]_`\`Q`$`cklmnocp#ppAx!}SBZqRU|a#h'>BrE &J!KcRTUWRWjYZ\]Y]o!]`hijkl`m mmnpqrsuvrvpvwDG)QT$qHQqM))43r"u:C1H#bQSe*0T+UX[\^_a\aXbcdfgcgXg+gjrstuvjw*wwAx!}SBZqRU|a#h'>BrE &J!KcRTUWRWjYZ\]Y]o!]`hijkl`m mmnpqrsuvrvpvwxxyy zz{||| @BCDEGHDHBHI IJ 1Q4$q'$qHQqM))43r"u:C1H#bQSe*0T+UX[\^_a\aXbcdfgcgXg+gjrstuvjw*wwAx!}SBZqRU|a#h'>BrE &J!KcRTUWRWjYZ\]Y]o!]`hijkl`m mmnpqrsuvrvpvwx xxy zq'$r!Xa]""D!C2J,S*A3r"u:)M$NQTUWXZUZQ[\]_`\`Q`$`cklmnocp#ppAx!}SBZqRU|a#h'>BrE &J!KcRTUWRWjYZ\]Y]o!]`hijkl`m mmnpqrsuvrvpvwDG)QT$qHQqM))43r"u:C1H#bQSe*0T+UX[\^_a\aXbcdfgcgXg+gjrstuvjw*wwAx!}SBZqRU|a#h'>BrE &J!KcRTUWRWjYZ\]Y]o!]`hijkl`m mmnpqrsuvrvpvwxxyy zz{|| | @BCDEGHDHBHI  IJ 1Q4aa8Aa=(($#be* q3x0GRPRU/S*TWZ[]^`[`WabcefbfWf*fiqrstuiv)vvAx!}SBZqRU|a#h'>BrE &J!KcRTUWRWjYZ\]Y]o!]`hijkl`m mmnpqrsuvrvpvwx xxy zq'$r!Xa]""D!C2J,S*A3r"u:)M$NQTUWXZUZQ[\]_`\`Q`$`cklmnocp#ppAx!}SBZqRU|a#h'>BrE &J!KcRTUWRWjYZ\]Y]o!]`hijkl`m mmnAqD& T!W9ad4!Xa]"2 2D!C2J,S:QSVWYZ\W\S]9]4^adeghjejaklmoplpap4ps{|}~tA4A!AAx!}SBZqRU|a#h'>BrE &J!KcRTUWRWjYZ\]Y]o!]`hijkl`m mm!npqrsuvrvpv!wxx yy zz{|| | @BCDEGHDHBHI  IJ 1Q4aa8Aa=(($#be* q3x0GRPRU/S*TWZ[]^`[`WabcefbfWf*fiqrstuiv)vvAx!}SBZqRU|a#h'>BrE &J!KcRTUWRWjYZ\]Y]o!]`hijkl`m mmnpqrsuvrvpvwx xxy zq'$r!Xa]""D!C2J,S*A3r"u:)M$NQTUWXZUZQ[\]_`\`Q`$`cklmnocp#ppstuvyABCDEyFvFtFGJKMNPKPGQtQSWYZ[^_abd_d[eYec(Y 2JXSbe*QT/S*,4QqMR:t::<=adFCEFtAwYPQSTPTUYZ[\]_ghijk_l\lZlnrtuvyz|}zwAuADEHKDLuLNQRTUWRWNXtXoYbe*QT/oQqMn"["$%a!A&6$6s2b5z$A43rRTu:XY\_X`I`befhikfkblHlCmpstvwytypz{|~{pCCKLMNOCPCP%P[PSTUVXYUYSY[YVZQZF[[ \\]^^ ^abdefgijfjdjak  kl%m.!5>d1Q47lQz!A$2221Q42r!Q$w{+,-Cei irVcPtd\}}}dttidtttjtzttdzzttdzzdttzz dtzz gddddtttjtzttdzzttdzzdttzz dtzz gdd d t td zt tttzd td zzttdzzzdtdzzdzttd zzz dgdd dttjttttt tzz dztdzz z t tzgdgddttjtdzt tttgdgddt ttjttdzttzt tttttd ztdd zzgddttjt|ttzz|t|tzzz t tttt ||zz|tt||zzzzt| tzz||zz t|| ft |tzzt| tzzdfgddttjtdtzttzztttdz zz t ttttdzdz zttdz zgddddtzzttjtzdzdttt zz t ttt tdtzdzz dzgddtdzttjtzttz tdzttdtz z zzt tttttz tdtz zgdgddtdzttdzztt tttzzt ttdttdzzz gddttttdtzttzztttdz zz t ttttdz ttdzdz zzgddttttttttzztdtzz t ttttt tttdtztttztzzgddttttttttzztdtzdz z t tttttt zttzdz gddtttttzttztttzz t ttttttzz ttztz gdd tttttzttz tt tz z t ttttt t tzzgdd!ttjtttz ttjtttzzt tttttzt ttttt zgdttjtttjtttz zt tttttzt tttgdttjtttz ttjtttzzt ttttt zt tttttzgdttjtttjtttzzt tttt ttttt zgdd"d#S)$NrAr lie_group_01rrT)rErFrKrL lie_group_02r lie_group_03rrrr lie_group_04ro lie_group_05r lie_group_06rr lie_group_07 lie_group_08 lie_group_09 lie_group_10r lie_group_11 lie_group_12 lie_group_13 lie_group_14 lie_group_15 lie_group_16 lie_group_17) lie_group_18 lie_group_19 lie_group_20rz)rr?r<rr rrrrrrr rr r%)rrTrUs rTrr sUgGAq!aDm! ! nad1f%1q11QqT69!A#= $ m! ! nad1f%1q11QqT69!A#=mAj1q))AadF1Q47N:a1fqj!A$PQ'=QQSTUm*diiladQs1vXo1!Q$66Q?+m6diilAo1Q4 ~7mB11q!tAaDy)1Q4C1Ix1~5567CmLdiilQqtV#aAaCj01Q4bS!a%[1SAE^3C%CS!AY$NPQTUPU$V 1qb $cAaCiaRT*D1345MmXdiilQqS1X%#q!te* 41Q4"q!tAv+sAqD5z123Ymb1Q3w11&*QsAaD5z\91Q4R1q\A-./0cmldAaDIIaL!AaD(1a4A1I+>>1Q4"s1v+s2a4y012mmxd1Q47lQz!A$2221Q4BAI'(ymB1Q4lQqS1X%#q!te* 41Q4adUR!Q$q&[123CmL1Q4lQqT#a&[(3qs831Q4c!fWrHS1Xc#a&k5I1,M'MNOPMmV1Q4lQqT#a&[(3qs8A:51Q4CQL3q61A567Wm`QqT! nqt#aAh.1Q4"qQx-#a&0!345amjQqT! nqt#aAh.1Q4BSV,-./kmttyy|AaD QqTYYq\!A$%671Q4CF#R!baRj%9:um@diiladiilQqT121Q4CF#R!b\2 tyy|AaD QqTYYq\!A$%671Q4CG$b1r#a&y&9: diilAaDIIaL1-.1Q4 b1r#qb'z23Smq qrVcxdttttjtdtttzz ttttt tzt ttzzgdttjtddtzttzzttttt dtddz z tzzt tdtddz z tzzzgdddS)N2nd_linear_airyrrrr)2nd_lin_airy_012nd_lin_airy_02rz) r?r<rr rr*rr+r rWrVrTrr{ s&aDdii1o!A$&1Q4F1I6!9 456 dii1o!AaD(1Q4FA!QK<>22RAaDF |A~8N5NNOP  rVc Utdd}tdd}tdDcgc]"}ttdzdtzzdz |$c}\}}}}}tdDcgc]"}ttdzd tzz d z|$c}\}} } } } tdDcgc]/}ttdzd td zzz d tzzd z|1c}\} }}}}tdDcgc]"}ttdztd zz dz|$c}\}}}}}tdDcgc]}ttdztz d z|!c}\}}}}}t d }dt tidt tj tddt tj tzztt tttt dtzzzgddt tj tdd t tj tzz dt tzztt tttt tzzt tzgddt tj tdt tz tt ttt t ztt tzzgddt tj td t tj tdzdt tj tzz tt tttt dtzzztt dtzzzgddddt tj tdzdt tj tzz d t tzztt ttt tdz ztt ttd d zzzgdddtt tj tddt tj tzzt tz dtt ttt tdtdzzztt t tdd zzzzgdddt t ttd t t ttdzdt t ttzz dt tzz tt ttt d tzztt t dtdzzzztt tdtdzzzzgdddt tj td t tj td z d t tj tdzz d t tj tzztt tttt dtzzztt tzztt dtzzzgdddt tj td d t tj td zzt tj tdzd t tj tzz dt tzz tt ttt t ztt tzztt td tzztt tdtzzzt dtzzzgddd t tj td |dzt tzz tt ttt t| tzztt t|tzzzttt|tzzzttt|tzzzgddd!t tj tdd|zt tj tzz dt tzz tt ttt t|t|dzdzz zztt t|t|dzdzzzzzgddd"t tj tdd |zt tj tzzd#|dzzt tzz tt ttt d$|ztzztt d|ztzzzgddd%t tj td tt ttttzzttdzzzttd zzzgddd&t tj tdd t tj tzzd t tzztt ttttzzt dtzzgddd'd t tj td zdt tj tdzzt tj tzt tz tt ttttzzt t ztt td z zzgddd(t tj td dt tj tdzz d#t tj tzzd)t tzz tt ttttttzzzzt dtzzgddd*t tj tdd|zt tj tzz |dzt tzztt ttttzzt |tzzgddid+t tj td d t tj td zztt ttttzzttdzzztt dtzzzgddd,t tj td dt tj tdzz tt ttttzztt td tzzztt tdtzzzgddd-t tj td dt tj td zzdt tj tdzz d#t tj tzz d.t tzztt ttttzzt dtzztttzzt dtzzzgddd/d.t tj td zd0t tj tdzz d t tj tzzdt tzztt ttt t ztt t d z zztt tdz zztt ttddzzzgddd1t tj td d)t tj tdzz d2t tzztt ttttzzt dtzztttzzt dtzzzgddd3t tj tddt tj tzz dt tzztt tttdtzzttdtzzzt tzgddd4t tj tdt tj tz t tztt tttttd zdz zttttd zdz zzt tdz zgddd5t tj td dt tj tdzzdt tzztt ttttdtzztttd tzzztttdtzzztttd tzzzgddd6t tj tdd t tj tzz d7t tzztt tttd tzzttd tzzzt dtzzgddd8t tj td d t tj tdzzd t tzztt ttttzztttdzztttzztttdzzzgddd9t tj td d)t tzztt tttttd zzttttd zzzt tztt dtzzzgddd:t tj td d t tj tdzztt ttttzzttdtzzzttdtzzzgddd;t tj tddt tj td zzt tj tztt tttttzzttzztt tzzttzzgdddt tj tddt tj tzzdt tzz tt tt t |tzzt t#|tzttt%|tzzttt%|tzzzzzt t#|tzttt%|tzzttt%|tzzzzzgdd?d@t tj tdd t tj tzz t tztt ttt |tzztt | tzzzt t | tzzzt t#| tzttt%| tzzttt%| tzzzzzgdd?t tj tdd t tj td zz d t tj tzzt tztt ttt | tzztt |tzzztt |tzzztt |tzzzt t |tzzzgdd?t tj tddt tj tdzz dt tj td zzdt tj tzzdAt tzz tt tt t dtzzt&t t|zzzt t#|tzttt%|tzzttt%|tzzzzzt t#|tzttt%|tzzttt%|tzzzzzgdd?t tj tddt tj tdzz dt tj tdzzt tj tdzdt tj td zz t tztt ttttzzt t|zztttzzttt%|zz t t&tzzttt%|zzzt tt#|zzzt(t*tzzttt%|zz t,tzt.zttt%|zzzt tt#|zzzgdd?ttdt tj tttzt tj tzdtt ttttdtd d ztzdz zzttdtd d ztzdz zzgdt|t tj tttzt tj tzdtt tttttt|z zztttt|z zzgdtt0t tj tttzt tj tzdtt ttttttt0z zzttttt0z zzgdtt2t tj tttzt tj tzdtt tttt tt2 tzzztt tt2tzzzgdt tj ttt4t tzztt ttt t tt4 zztt ttt4 zzzgdtdBt tj ttzt tzdtt ttttd z ztttd z zzgdtdBt tj ttzt ttt ttt t d z ztt td z zzgdtt tj ttd t t ttzz dt tzzdtt tttt tzzt tzgdtt tj ttd t t ttzz d t tzzdtt ttttzzt dtzzgdtt tj ttdt t ttzzd t tzzdtt tttttdzzttttdzzzt t zgdt tj ttd t tzztt tttdtzzttdtzzzgddCdDScc}wcc}wcc}wcc}wcc}w)ENrTrr)realrrrrrdirr%nth_linear_constant_coeff_homogeneouslin_const_coeff_hom_01rrlin_const_coeff_hom_02lin_const_coeff_hom_03lin_const_coeff_hom_04rkrjrlin_const_coeff_hom_05lin_const_coeff_hom_06rrlin_const_coeff_hom_07lin_const_coeff_hom_08lin_const_coeff_hom_09lin_const_coeff_hom_10lin_const_coeff_hom_11lin_const_coeff_hom_12rbilin_const_coeff_hom_13lin_const_coeff_hom_14lin_const_coeff_hom_15lin_const_coeff_hom_16ralin_const_coeff_hom_17lin_const_coeff_hom_18lin_const_coeff_hom_19lin_const_coeff_hom_20rQlin_const_coeff_hom_21%lin_const_coeff_hom_22rlin_const_coeff_hom_23lin_const_coeff_hom_24lin_const_coeff_hom_25lin_const_coeff_hom_26rlin_const_coeff_hom_27lin_const_coeff_hom_28lin_const_coeff_hom_29lin_const_coeff_hom_30lin_const_coeff_hom_31lin_const_coeff_hom_32lin_const_coeff_hom_33rplin_const_coeff_hom_342rm)lin_const_coeff_hom_35lin_const_coeff_hom_36lin_const_coeff_hom_37lin_const_coeff_hom_38lin_const_coeff_hom_39lin_const_coeff_hom_40lin_const_coeff_hom_41lin_const_coeff_hom_42lin_const_coeff_hom_43lin_const_coeff_hom_44lin_const_coeff_hom_45lin_const_coeff_hom_46lin_const_coeff_hom_47lin_const_coeff_hom_48rz)r)rr'r<rr?rr rrr{rrr|r rr}rrr~C7C8C10C9r rr>)rrrrrrrrr6r7r8r9r10r11r12r13r14r15r16r17r18r19r20r21r22r23r24r25rs rTrr s{ sT"AsA>CAhG&A1q!4GBBB>CAhG6!Q$1*q.!4GBBCPUVWPXY1vadSAXoQ&>&BAFYCc3DI!HMqvadQTkB&6:MCc3@EaI1vadQhlA6ICc3 AA;aDjdii1o!A$))A,.1Q4bRTl*+,j dii1o!A$))A,.1Q471Q4"r#a&y.#a&012 jdii1o!$1Q4CGbQi/01j dii1o! !Q/!AaDIIaL.@1Q4bRTl*RAaC[89:!j,! !Q"QqTYYq\/1AadF:1Q4C!Hr#aA.>*?'??@A-j811a1QqTYYq\>1AaD8!<1Q4C2Q< 011BsA2tAw{;K7L4LLMN9jD1Q4AadAq!11BtAaD!}4DDq1vM1Q4C!Hr#qb!d1g+.>*?'??"SBQUVWQXLIYEZBZZ[\EjPdii1o! !Q/!AaDIIaO2CC !A$))A,1Q4bRTl*RAY6C!HDEF Qj^dii1o!A$))Aq/ 11AaDIIaOC !A$))A,1Q4 1Q4CGbQi/2c47(1*o3E3tTUwWXy>HY3Y[^_abc_c[d2ddef _jldii1o1QqT )1Q4Ca O+bT!WQY.??"SaQR^BSSVXY\]abc]def]fYgVgghimjxdii1o!AaDIIaL 001QqT691Q4C1tAqD1H~#5 677"SAQPQTTUXDVAW=X:XXYZyjDdii1o!AaDIIaL 002ad71Q4<?1Q4C1QK"S1QZ-789EjPdii1o1Q4bdR1W,r!Q$w678Qj\dii1o!A$))A,.1Q471Q4"r!t)SAY./0]jh! !Q!AaDIIaO"33adiilBQqTI1Q4"r!t)S!W,r#ac({:;<ijtdii1o!A$))Aq/ 11Bqtyy|OCa!fL1Q4"q"r!t)},c!A#h678uj@dii1o!AaDIIaL 001a4!9<1Q4"r!t)S1X-./AjLdii1o!A$))Aq/ 111Q4bdR1W,r#bd)|;<=MjXdii1o!A$))Aq/ 111Q4bdRT!WHQJ%77"Sa^:KKLMYjddii1o!A$))Aq/ 11BqtyyA4FF 1Q499Q<QqT'"1Q4"r!t)SAY."r!t)S1X1EEFG ejr11a 2adii1o#55!A$))A,F1Q4O1Q4CGbaRTl2RAaC[@2c!HUVXYNJZF[C[[\]sj~dii1o!A$))Aq/ 11BqtG;1Q4"r!t)SAY."r!t)S1X1EEFGjJdii1o!A$))A,.1Q471Q4"S1X+3qs8 3SV;<=KjVdii1o! ! ,qt31Q4"S471--3qay{3C0CCS1XMNOWjbdii1o!A$))Aq/ 11AadF:1Q4 3tAwqy>Bs4719~--3tAwqy>0AABs4PQ7ST9~DUUWX cjpdii1o!A$))A,.AaD81Q4"S1X+3qs8 3S1X=>?qj|dii1o!A$))Aq/ 11AadF:1Q4"r!t)S47^3rBqDy#aQi.6PPQR}jHdii1o!A$&1Q4"S47^+bQtAwY.??QG"SQSTUQUY,VWXIjTdii1o!A$))Aq/ 111Q4bdRAaC[02c!A#h;>?@Uj`dii1o!A$))Aq/ 11AaDIIaL@1Q4rBqDy#a&00BAIs1v3EEFGajldii1o! !Q/!A$61Q4"Sa1--3tAwqy{3C0CCS!AYN c$q'!)A+ CQ ! $4!4 4c!A#h >?@A mjzdii1o!A$))Aq/ 11AaD81Q4C$Qx!|"4 4553qd1gPQkAR?R;S8SS S4a1 %% & &')+C$tAw{2C0C,D)DEFG {jJdii1o1Q499Q</!AaD&81Q4 3r!t9 s2b6!8}3r"vax=(82c"R&(m;K(KLL bfQh-2c"R&(m+bRVAX.>> ? @AB! Kj\dii!nq11~-!41Q4 3r!t9 r#bd)|#bRTl2 bfQh-2c"R&(m+bRVAX.>> ? @AB! ]jpdii!ns1Q499Qq>11D114EE!L1Q4CAJCAJ6CAJFCPSTUPUJVY[\_`cde`e\fYffgh dii1o!A$))Aq/ 11Aadii1o4EE1Q499UV<WZ\_`ab_cZcc1Q4 3qs8 S3Z-  bgai.Bs2c719~-3r#wqy>0AA B C bgai.Bs2c719~-3r#wqy>0AA B CDE !dii21Q499Q?!22QqtyyA5FF A$))Aq/adii1o-.01!51Q4"r!t)S3Z/b2a4iQr#wY5O3P 19c!BsG)n $4%&)!BsG)n355;=19c!BsG)n:T8U q52:s1RW9~ %9&'*1RW9~86678! a1Q499Qq++adiil:A>1Q4bQA%6q%8%:!;;;bQQRTUEVWXEXYZEZA[>[[\] QqTYYq1%%! ! 4a81Q4bQtAwY//"S47^2CCDE adii!A&&115q91Q4bQtBxZ002c!DH*o3EEFG QqTYYq1%%! ! 4a81Q4bd1gXaZ002c$q'!)n3DDEF dii1o!A$&1Q4C48 ,,r#aaRj//AABC 1Q499Q?"QqT)1-1Q4C!Hr#ac({234 1Q499Q?"AaD)1Q4C1I3qs8 345 11a1T!A$]?2QqtV;Q?1Q4"r#a&y.#a&012 11a1T!A$]?2QqtV;Q?1Q4"r!t)S1X-./11a1T!A$]?2QqtV;Q?1Q4"S47^+bQtAwY.??aRHIJdii1o!A$&1Q4C!Hr#ac({234M jn n HGYMIs"'Bi,%'Bi1"4Bi6,'Bi;)$Bjc dtttttz ttttz ztttz ttttz zttttz zttj tzz t t ttt ttttttz ztz z gddtttj tzttz tttttz zz t t tt tt ttttz dz dz zt ttttz dzdz z gddtttttz zttttttz zz tttttz zttj tzzt t ttttt tz ttttz zdz zttt tz ttttz zdz zgddttj ttttz z dttttz z zt tttttt tz dtzzzt tttttt tzzgddtttttz zttztttj tzz t ttt dtt tz z tzgddd d S) NrTrrrFryrp)dep_div_indep_01dep_div_indep_02dep_div_indep_03dep_div_indep_04dep_div_indep_05rz) r?r<rr rr rrrrr rWrVrTrr s?aDd1fS1a[ AadF3qtAv;$6QqT!V$DadiiPQl#RR3q62AaDQqT!V$4Q$6 7789! ! nqt#aAaDF m33q63r7SQqT!Vq%9!%;;c#ad1f+PQ/>RST>TTUVAaDF mad3qtAv;..3qtAv;qtyy|1KK3q62QqTE!G S1a[ 8 ::S!A$q\#aPQdSTf+=UVW=WWXYdiilQqT!V#aAaDF m31Q4T"s1v+..2567AaD!DcRSfDUBU9VWAaDF mad"Qqtyy|^31Q4ac!fo1234 9!% %rVc 2dttdttjtddttzztdzz ztttttdzdtzzdzz t ddz gdidS) Nlinear_coefficientslinear_coeff_01rrrkrmrrz)r?r<rr rrrWrVrTrr s*aDdiila!AaD&j1q5111Q4QTAaCZ!^,x1~=>?  rVc dtttttttttz zdtzz ttttzzt ttt t  tt t t dzzzgdddttzt tttz zttttjtzzdtzt tttz zttjtzz t tttt dt tttz zz gddtdzzttzttdzztttdzzdtdzzz ttjtzzt ttt dt ztdtd zzt d t zzdz ztz gddttttdztttzz zttjtzttz t tttt dtt ttz dzzz gdd tttztttz ttjtzz t ttt ttdttdztdzz zz ttttz zgdtttjtzttz tttttz zz t ttdtztt tzzgdtdzttdzzdtzttzttjtzz t ttttt tzz t ttttt tzzgdttdzttttdztdzz ztttzz ttjtzzt ttt tt tdt zz z zt ttt tt tdt zz z zgddd d S)NrrrTrprrrrrr)1st_homogeneous_coeff_best_011st_homogeneous_coeff_best_021st_homogeneous_coeff_best_031st_homogeneous_coeff_best_041st_homogeneous_coeff_best_051st_homogeneous_coeff_best_061st_homogeneous_coeff_best_071st_homogeneous_coeff_best_08rz) r?r<rrr rrrrrr rWrVrTrr s1aD daAaDF mac)4!a=881Q4#b'(A2c2#'l?";;<= &!fS1Q4[ 1Q4! ! #44qs3q1v;qtyyQR|7SS3qt9b1S1Q4[=012&1fQqTkAaD!G#q1qy1QT6'91Q499Q<&GG1Q4QrTHR1WSBZ-?$@$BBCAEFG &4!a!AaD&())1Q499Q<7!A$>3qt9b1T1"QqT'A+%6#6678&!A$h!ad(AaDIIaL003q62DQqT1WQT\)9$: ;;d1Q46lJKL& ! ! nqt#aAaDF m31Q41T"Q$Z()& dQqT1Wnqs1Q4x! ! 441Q4$q"q&z**+R!d1b1f:6F-GH& dAg4!a!Q$//!AaD&8!A$))A,FF1Q4"T1"a!B$h-0012adBtQBAbDM?Rr?r@rrr{r|r}r~rrrrrrTrUrrrrrrrXrbror{r}r]rjrrrrrrrrrrrrrrrrr r rrrrr!r$r(r,r1r3r6r8r;r>rArDrGrMrVr^rdrRr[rcr@r:r=rCr5rFrJrr r#r'rr rrrrrrrrrrrhrWrVrTrZs <333*".9GGKKA:[[[>/;*'HH72G::4- 3K 3K 3Z 3K SM SM+27+;'BBBBC ' 1a  * " .'J18M`<>4"*J@ G G[[ RRff<66>>RRCC-OOQ>>::]]I Q QJJ66A TTeeJJQQ @6677 KK4AA679 6 6 > E ; ;77t))X##L2  HHX  VVt''T00fTTnQQjDDNppf]]@eePAAH@//dkk\ssl$yyx &&P  55prV