gg 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(Zldrd*Zmdrd+Znd,Zod-Zpdrd.Zqgdfd/Zrd0Zsd1ZteQd2ZueQd3ZveQd4ZweQd5Zxd6ZyeQd7ZzeQd8Z{eQd9Z|eQd:Z}d;Z~eQd<Zd=ZeQd>ZeQd?ZeQd@ZdAZeQdBZeQdCZeQdDZdEZeQdFZeQdGZeQdHZeQdIZdJZeQdKZeQdLZeQdMZdNZeQdOZdPZdQZdRZeQdSZdTZdUZeQdVZeldWZeldXZeldYZeldZZeld[Zeld\Zeld]Zeld^Zeld_Zeld`ZeldaZeldbZeldcZelddZeldeZeldfZeldgZeldhZeldiZeldjZeldkZeldlZeldmZeldnZeldoZeldpZdqZdS)sa 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} c"fd}|S)Nc}g}|dD]6}|d|d|d|d|d|dg|d|d|d||d|dd|d|dd |d|d d|d|d d|d|d d|d d }||8|S)NexampleseqsolXFAILfuncr;F simplify_flagTcheckodesol_XFAILdsolve_too_slowcheckodesol_too_slowhint) rErFrGrH example_namer;rIrJrKrLrM)getappend)solverrDexampletemprHs _/var/www/html/ai-engine/env/lib/python3.11/site-packages/sympy/solvers/ode/tests/test_single.pyinnerz _add_example_keys..innersStvvj) " "GZ(1$7j)'259 +G488"EEz*7377vf~NN 'z*7377FF &z 27 ; ? ?PT U U%+J%7%@%D%DEXZ_%`%`"("4W"="A"ABSTY"Z"Z'-j'9''B'F'FG]^c'd'dv  D OOD ! ! ! !)rHrUs` rT_add_example_keysrXs)& 577NrVFc|D]P}|s|ds |r |dst|d|d}|ddkrt|dQdS)Nr;rMT) solver_flag xpass_msg)_test_particular_exampleprint) ode_examples run_slow_testrRresults rT_ode_solver_testrbs'' GFO  PWX^P_  )'&/7PTUUU + " $ $ &% & & & ''rVcttdgz}t}|D]-}|dsd|vrt |||.dS)Ndefault _Integralseries)listr2_get_all_examplesendswith_test_all_examples_for_one_hint)runxfail all_hints all_examplesour_hints rT_test_all_hintsrostX {*I$&&LJJ   [ ) ) X-A-A ',IIIIJJrVct|tkrtfd|DS|S)Nc3BK|]}|VdSN)dummy_eq).0sub_dsol expected_sols rT z"_test_dummy_sol..s1NNx<((22NNNNNNrV)typerganyrs)rv dsolve_sols` rT_test_dummy_solr{sM JNNNN:NNNNNN$$Z000rVcdt}|D]}|d|krt||dS)NrN)rhr])rnrNrmrRs rT_test_an_exampler}sG$&&L88 > "l 2 2 $Xw 7 7 788rVc|d}|d}|d}||dv}|d}ddd}|d} |d } |d } |d } d } |r>|t||vr,t||| }t||t||vr8||d< | st ||| |}nt |dkr |d}n|}n}#t $rp}g}||d<|stjtt|||||d<|r|st|dd} Yd}~nd}~wwxYw|rn|gkrgd}t|tkr:|D]6}|trt!|| }n||v}|rn7n4||v}|D]-}|trt!|| }.|r-t"||||}t|dt%t |D}t |dkrd}| ss| sqt'|||d|krZ||d<d} t(||||}|r+t*||}t|d|z|d<| r |r |dz|z|d<|S) NrErFrNrGrHr\)msgr[rIrJrKrLT)rRrErn match_listsimplifyrMrexception_list)erMrRrErF)rRrErFrzcg|]}dS)TrrWrtis rT z,_test_particular_example..s#P#P#P!I#P#P#PrVr)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_solexpected_checkodesols rTr]r]s T Bu%L.)G  G, ,E v Db ) )Fo.M#$78!"34O&'=> E* <D11 1 1"))'b8)TTG )) )<D))))&| # .#B}(SSS |$$a''!-aJJ!-J   J'.F# $ &#%%%)003q66RY^`0aaF5M 5 fUm$$$EEEEEE  " I:r>>$ J%%+G{{5))E/>w /S/S+S((+2*+D('$.\#A +TTG{{5))T/>w /S/S+S( ..55g"R^ku5vv$W---#P#PuS=N=N7O7O#P#P#P <  A%%'0$' I( I"2z4NNNRfff18~. %"4";";bVbnx";"y"y&I&5&<&,HF5M  UU U"),I"IH"TF;  Ms+4C E*A&EEcH|gkrt}ggg}}}|D]}||dv}|r|s|rt||}||dgz }||dgz }||dgz }|#|d} | dkrt|d|nt |} t |t |z t |z } t || | ||} t| dSdS)NrGrrrrr\)rMmatchedsolveunsolve exceptions)rhr]rOr^rcheck_hint_msgr) rnrmrkrrrrrrar match_countsolveds rTrjrj'sZr(** /12rn J# % % K00  E     )(K@@fjjb111  >"555 &**%5b999  -CBwwfUm$$$*oo Z\!2!223~3F3FF##;f^jxF#GG c rVc&ttttttt}t |t t fdt t fdt t fdt t fdttttttttzzttt}t|t tfdt| dusJttttttzttt}|j dksJttttdtttd ztttd z ttt}|j dksJttttd tttd zttd zz ttt}|j d ksJttttd ttttd zzttd zz ttt}|j dksJdS) Nc,Srr)matchesrQsrTz&test_SingleODESolver..Fs(8(8rVc,Srrget_general_solutionrsrTrz&test_SingleODESolver..Gs(C(C(E(ErVc,Srr)_matchesrsrTrz&test_SingleODESolver..Hs(9(9rVc,Srr)_get_general_solutionrsrTrz&test_SingleODESolver..Is(D(D(F(FrVc,SrrrrsrTrz&test_SingleODESolver..Rs&"="="?"?rVFrT) r6r?r<rr7r:NotImplementedErrorr4r5rorder is_autonomous)problemrQs @rTtest_SingleODESolverrAsqttyy||QqTT155G W % %F  8 8 8 8999  E E E EFFF  9 9 9 9:::  F F F FGGG qttyy||add1Q44i71qAAG  ! !F =????@@@  ! !F >>  u $ $ $ $qttyy||add2AaDD!<11QOQRSTQUQUWXYYG =A    qttyyA11a@1Q447JAaDDRSTTG  D ( ( ( (qttyyA1Q4499Q??1BBQqTT1WLaPQddTUVVG  E ) ) ) ) ) )rVc.ttdSrr)rb)_get_examples_ode_sol_linear_coefficientsrWrVrTtest_linear_coefficientsrfs>?????rVc2ttttttz z }t|tt}ttttzttz }t|tt}d}d}||vr||vsJ||vr||vsJt t t t dS)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 c1Q44 B AaDDIIaLL.1Q44 C c1Q44 B 4C 4C b==S]]]] b==S]]]]STTTEFFFFFrVc^ttdttddSNT)r`)rbrrrWrVrT,test_slow_examples_1st_homogeneous_coeff_oderys2ScghhhhEUYZZZZZZrVc.ttdSrrrb;_get_examples_ode_sol_nth_linear_constant_coeff_homogeneousrWrVrT*test_nth_linear_constant_coeff_homogeneousrsPQQQQQrVc2ttddSrrrWrVrT8test_slow_examples_nth_linear_constant_coeff_homogeneousrsP`deeeeeerVc.ttdSrr)rb%_get_examples_ode_sol_2nd_linear_airyrWrVrTtest_Airy_equationrs:;;;;;rVc.ttdSrr)rb_get_examples_ode_sol_lie_grouprWrVrTtest_lie_groupr455555rVcttt}tttz |zttdzttztdzttzdz z z}t |dksJt t dS)Nrr) factorableseparable_reduced lie_groupseparable_reduced_Integral)r?r<rr!r0rb'_get_examples_ode_sol_separable_reduced)dfrEs rTtest_separable_reducedrs 11B add(B#ad1Q44i1a4!9q=9:: :B    & & & & &<=====rVc2ttddSr)rbrrWrVrT$test_slow_examples_separable_reducedrs >1 EB QqTTBHS1qs7QAY2G)K)K%L%La A#ON&8Aq1uIq11!344558(1a G D ;;::661Q4462,a001 2 2C &"?@@@ @ @ @ @ r3  9 , , , , , ,rVc.ttdSrr)rb8_get_examples_ode_sol_2nd_nonlinear_autonomous_conservedrWrVrT'test_2nd_nonlinear_autonomous_conservedrsMNNNNNrVcvtttdtttz}t t dt tdt tttzz z tttfttzt t dt tdt tttzz z tttfttz g}t|dd}t||D]\}}| |sJt|d|DddgksJdS)Nrr+2nd_nonlinear_autonomous_conserved_IntegralFrMrc6g|]}|SrW)doit)rtss rTrzDtest_2nd_nonlinear_autonomous_conserved_integral..s 555AFFHH555rVr)r?r<rrr r%rr_urr1ziprsr9)rEactualrars rT0test_2nd_nonlinear_autonomous_conserved_integralrsL 11a4!:: %B!Dab2(>(>&>!>???"addLLbSTfUUx$rAhtBxx4444555AaDDzBBBFKKMF BJUZ [ [ [F66""!zz!}} r55f555 6 69i:P P P P P P PrVc.ttdSrr)rb'_get_examples_ode_sol_2nd_linear_besselrWrVrTtest_2nd_linear_bessel_equationrs<=====rVc tttttttzz}tttt ttttt t t tzzg}t||dt}|tttt t t tzzgksJttdSNr) r?r<rr rrrr3rb#_get_examples_ode_sol_nth_algebraic)eqnsolns solns_finals rTtest_nth_algebraicr s A$$1addiill" "C !c!ff   !bRTl # # %E.sE1a@@K 2addBs2a4yyL112 2 2 2 2899999rVc2ttddSrrb2_get_examples_ode_sol_nth_linear_var_of_parametersrWrVrT>test_slow_examples_nth_linear_constant_coeff_var_of_parametersrsGW[\\\\\\rVc.ttdSrrr rWrVrT0test_nth_linear_constant_coeff_var_of_parametersrsGHHHHHrVcxd}tttddtttdzztttzdtzz tttzz }t |tt|d}t |tt|d}||ksJt ||ddd ksJt ||ddd ksJdS) 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 11a1QqTTYYq!__, ,qttyy|| ;ac AC!HH LBb!A$$X===Hr1Q44h???I y r81U C C Cy P P P P r81U C C Cy P P P P P PrVc2ttddSrrb_get_examples_ode_sol_1st_exactrWrVrTtest_slow_examples_1st_exactr4DIIIIIIrVc.ttdSrrrrWrVrTtest_1st_exactrrrVctttttttzttdzz tttzz }t |ttdd}t ||ddsJdS)NrF1st_exact_Integralrrr)rr?r<r rr1r9)rEsol_1s rTtest_1st_exact_integralr"s QqTTaAaDD kAaDD!G+QqTTYYq\\9 9B 2qtte2F G G GE r5% @ @ @@@@@@rVc2ttddSr)rb)_get_examples_ode_sol_nth_order_reduciblerWrVrT&test_slow_examples_nth_order_reducibler%s>dSSSSSSrVc2ttddSr)rb:_get_examples_ode_sol_nth_linear_undetermined_coefficientsrWrVrTFtest_slow_examples_nth_linear_constant_coeff_undetermined_coefficientsr(sO_cddddddrVc2ttddSrrb_get_examples_ode_sol_separablerWrVrTtest_slow_examples_separabler,rrVc ttttttttzt jztz }d}|t|vsJttdS)N3nth_linear_constant_coeff_undetermined_coefficients) r rr?r<rr Halfr0rbr')rErns rT8test_nth_linear_constant_coeff_undetermined_coefficientsr0si D1qMM1QqTT6AF?Q. / /BDH |B'' ' ' ' 'OPPPPPrVc 4d}t}||tttztt|tttzdksJ||tttztt|tttzdksJ|tt|tttz|tttdzdksJ||tttztd|ttzttztdzdksJ||tttd|tttdz|tttdzdksJ||tttd|tttdzdksJt t dS)Nctt|tttSrr)r8r6r?r<r)rEs rTrz*test_nth_order_reducible..s.$%5b!A$$%B%BCCLLNNrVFrrrT)rr>r?r<r=rbr$)FDs rTtest_nth_order_reducibler5sNNAA 1QQq1vq!__qq1qzz ) * *e 3 3 3 3 1QQq1vq!__qq1qzz ) * *e 3 3 3 3 1QqTT!!AaDD!**_qq1q!}} , - - 5 5 5 5 1QQq1vq!__qq1QqTT1a00 0 1 1U : : : : 1QQqttQ]]QQqttQ]] *QQqttQ]] : ; ;u D D D D 1QQqttQ]]QQqttQ]] * + +t 3 3 3 3>?????rVc.ttdSrrr*rWrVrTtest_separabler7rrVc 6tttdtztz tt tttdtzzt ksJt tdSr)r&rr?r<r r%rb _get_examples_ode_sol_factorablerWrVrTtest_factorabler: sl d1QqS66"9oo%q ) )htB1Q3K7H7H!.L.L-L L L L L566666rVc2ttddSr)rbr9rWrVrTtest_slow_examples_factorabler<&s5TJJJJJJrVc.ttdSrr)rb_get_examples_ode_sol_riccatirWrVrTtest_Riccati_special_minus2r?+s233333rVc.ttdSrr)rb*_get_examples_ode_sol_1st_rational_riccatirWrVrTtest_1st_rational_riccatirB/s?@@@@@rVc.ttdSrr)rb_get_examples_ode_sol_bernoullirWrVrTtest_BernoullirE4s455555rVc.ttdSrr)rb _get_examples_ode_sol_1st_linearrWrVrTtest_1st_linearrH8s566666rVc.ttdSrr)rb#_get_examples_ode_sol_almost_linearrWrVrTtest_almost_linearrK<s788888rVc `d}ttttttz tttttttzdz zttttdzdz z tt}ttttttz tttttdz ztttttdzzdz z tt}||vsJ||vsJ|dz|vsJ|dz|vsJt t dS)N Liouvillerre)r0rr?r<rb_get_examples_ode_sol_liouville)rMnot_Liouville1not_Liouville2s rTtest_Liouville_ODErQ@sJ D!$qttQ--/AaDDaddAq9I9I4I!4K"K QqTT1 q#aDD""N!$qttQ--/D1q!4D4DQ4F"F $qttQ-- 1#dd$$N ~ % % % % ~ % % % % + ^ 3 3 3 3 + ^ 3 3 3 3455555rVc(td\}}}}}td}d}tt||d|dzzdtt||dz|dzzz dt|zz}|t |vsJ|||z||zt|||zz||dzzt|||dzz}|t |vsJt t dS)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_homogeneousrZOsK((MAq!Q A0H addAq  !Q$ D1q!$4$4!4QT!9 9BqttG CB |B'' ' ' ' ' 11Q44!A#d11Q44mm# #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 dS) Nzx tza b c dTinteger>DAqD111JAq!QMH Ad1Q44A AqDaddAq)9)9!9 9BqttG Ca GB |B!-- - - - - 1a4QqTT1a  1Q3tAaDD!}}#4 4q1v =#a&& HB |B!-- - - - -CDDDDDrVc &td\}}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 tdS) Nzx, tz a, b, c, dTr\: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_parametersrhksk 6??DAqt444JAq!QKH AqDadd1Q !A#d1Q44ll"2 2R!W > !AadF4!Qq>>$9 9AaCQqTT! .|s\gggYZ&A1a4"QT'!9C1H!Ds1u!Lq!PRSTTgggrVrrSrrrrrri=rf 2nd_power_series_regularr^rErFrGTrErFrJ) euler_hom_01 euler_hom_02 euler_hom_03 euler_hom_04 euler_hom_05 euler_hom_06 euler_hom_07 euler_hom_08 euler_hom_09rMrHrD)rr?r<r rrrrrrC3r rC4C5C6r)r1r2r3r4r5s rTrVrVzsgg^cde^f^fgggBBB5aDDD1qMM!!#a1fT!A$$1-=-=&==qAA1Q44bHQNN!222334 1Q44!D1qMM/!++a1fT!A$$15E5E.EEqII1Q44DGGbAg-../ 1Q44!D1qMM/!++ad4!a3C3C.CCQGG1Q44"r#a&&y.!Q$.//0 1Q44!D1qMM/!++a1fT!A$$15E5E.EE1TRSTURVRVXY[\^_M`M`H``bcdd1Q44AqD2a4"QT'1223 addRQqTT1 -a//"QT'$qttQ:J:J2JJQPQTRVWXYZW[W[]^`acdReReMeeghii1Q44ArBs1vvI~3q6619 <=>>? d4!a###aQqTT1 o5!A$$>1Q44ArEBq!tG+,,- !ffQTk!A$$))Aq//)CFF1HQqTTYYq\\,AACFF1Q44KO1Q44CAKK"SQ[[.899:,-kldQqTTYYq!__$q11~5!<1Q44A1b52ae83bBh>ArEIBqRTuHTUUV 1Q44Aq))1Q44bdRT)**+[115 5 5rVcdttttdztttttztttttzzdttttt 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 gddtdzztttdzttztdtztt dtzdz zztttttttt tdz zt tt tdz zztdt tztt dtzdz zdz zzgddtdzztttdzttztt dtzdz zttttttztt tdz zt ttztt tdz zzdtt dtzdz zdz z dtt dtzdz zdz z gdddS)Nr^rrrrrre @rrs rUr)euler_undet_01euler_undet_02euler_undet_03euler_undet_04euler_undet_05euler_undet_06euler_undet_07r) r?r<r rrrrrrrr rrWrVrTr_r_sSaDDAd1Q44A&&&4!a==8!<<1Q44bQi#a&&!)A+5667 Ad1Q44A&&&1T!A$$]])::QqttVCQTJJ1Q44BAIAq!t(;;<==> Ad1Q44A&&&4!a==81QqTT6A3q66!8LL1Q44"r!Q$w,Q14s1vvby@!CDDE Ad1Q44A&&&1T!A$$]])::QqttVCSVVQYQTUVQWQWEWXX1Q44AqD2ad7*Xa]]3q6619-DDxPQRT~~VYZ[V\V\^_V_G__bjklmobpbpqtuvqwqwbww{CDEGJ{K{KKLLM Ad1Q44Aq)))AadF4!a3C3C,CCac$qQRttUV--FWWZ[\]^_\`\`Z``befgbhbhii1Q44A1a4"QT'1HQNN3q664IIHUWY[L\L\\]]^1fQqTTYYq!__$qtt+d1Q3iiC!HHQJ.GG1Q44a"SQ]]"2RCFF1H 5E"EQPSTUPVPVWZ[^_`ab_b[c[cde[eWfWfHfghHh"hijjk 1fQqTTYYq!__$qtt+c#ac((1*oo=1Q44DGGCAqMM1BtAwwJs3q66!8}}4LLqQTUXYZ[\Y\U]U]^_U_Q`Q`O`abObbefgjknopqrorkskstukugvgvevwxexxyyzA$$( ( (rVc< 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 gdtdztttdzttttzzdttzzdtz z tttttdttzzt tdttzzzdd tzz zgdd d S) Nrdrrrrsrrerr) euler_var_01 euler_var_02 euler_var_03 euler_var_04 euler_var_05 euler_var_06r) r?r<r rrrrrrrr"r rrWrVrTrgrgsSQaDDAj1q!,,,qs:addA3F3F/FF1Q44OQRTUQUVV1Q44BAI1Q./001 1a4QqTT1a(((1Q3tAaDD!}}+<>? d:addAq)))AaC 1Q440C0C,CCa!fLsSTvvU1Q44A1a4#a&&(2Xa^^CDDE 1vvg:addQF333j1q6I6II1LL1Q44bQi#a&&02a55899: dQqTTYYq!__$q1Q4499Q<<'77!add(BQqSH1Q44C#a&&MM)Bs1SVV8}},<rrrWrVrTrDrD s_ aDD1Q4499Q<<!A$$&1q0!441Q44BqD1H&&'%&ddiillQqttV#1Q44C!HH%%& dd1Q4499Q<<!#1Q44$rAaCx..))2addDacNN+C+CD  rVcdttddtttzttdzztttz z dtdzzztttt d t t t dttzdz zzdzdtzz gdidS) NRiccati_special_minus2 riccati_01rrrrr)r?r<rr rr!rrrWrVrTr>r>#s-aDD! ! nqttQw&1a/!AG);1Q44477(3rDGGCFFN1,<'<#=#==AAaCHIIJ   rVcddttttttttdzzdz dtttt dt t dt dztzz zt t dt dztzz z gdttdzttttzdttztz zdtdzz ztttdt ztz tt tzzz gddtdzzttttztdttzttttzdz zz ttdz ttzztttt dtdzzzt tzz gdtttdz 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_riccatirrrrrrrsrrrru<^H!re9:iiirrf(rkrTii| ri0H iiFi4"i<ii2ҍi"0 iimi^iHi iq!tVac\A%1q51* 561Q442#1*q!t+ac1A51r AqD8H1a48ORSUVRV8V aC991q!!AaC1q1$5a$771Q441"R%!AqD&.1,rAqDy1}=>>?AaDD!$$a!fadd1fqsl&;qsQqTT!V|QqTTSTW>T qSUQqTT!V^?' ttAv'&),QTDF):S)@2ad7RPQSTPTWCT(U'VWW1Q44"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!HH*MPQQTUVWXUXQYQYzZ"R%2+Bq"u ,s2vae|;d2gaemKeTVhWXZ[W[m[Ri1n%byA~.07 1a4@BI"*QPQT/R RK1 '{1a4/02:2+ad2BCEMb[QR]S RKac((#$ %&&'  kAAE E ErVc`dttdtttttttzztdztttttt tdzdz zzt dt tztt dtzdz zdz z t tdz dz zgdidS)N 1st_linear linear_01rrr) r?r<r rrrrr r#rWrVrTrGrGs!aDD11!A$$&1--1Q44"qQT!V},d1ggd2hh.>tDGGAIaK?P?P.PQR.RRTWYZ\]Y]X]^_X_T`T``aab   rVcD0td}td\}}}}}dttidtttttttzzt ttdt ttt tz ggdddttttttttzzd zzt ttt td ttzttd tzd z zz ttd z d z zt ttdggd dd ttttttd zzztttd tttzzzt ttt tt zttt zzt ttt ttd z d z zgd dttttttd zzztttd ttzzt ttt ttd z d z zt ttt ttztttzzgd dtttd zdz tttd zdz zt ttt tz t ttt tzt ttt d tzzt ttt d tzz gd dtttd tttz tttzt ttt!t ttt ttzzdzz t ttt!t ttt ttz zdzz t ttt gdddtttd zdz tttttzdz 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 tzzgdddtttd tttz tttd 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 ttttttd zz t t)dt!|d zz |ttft tzgd!gdd"tttd zttd zz t ttdt d zd t ztzz td zzz gd!gdd#tttd zttz t ttt d zdz t tzd z z td zdz zgd d$tttd zttd zz t ttt ttzt ttt tt zgd d%tttd zttd zz t ttdt d zd t ztzz td zzz gd d&ttt|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 ttttd 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_bestrrrr^1nth_linear_constant_coeff_variation_of_parametersrdr.rwfact_02r)rrrfact_03rrfact_04fact_05rrfact_06TrErFr;fact_07fact_08fact_09rsrfact_10rrffact_11rrErFrKfact_12rfact_13fact_14fact_15fact_16fact_17grx)fact_18fact_19fact_20r)r rr?r<rr rrr r#rr*rr+r rrrr,r-r/r%rr)r>a0a1a2a3a4s rTr9r9s c A122NBr"R aDDK ddQqTT!A$$))A,,&&1Q44 R!b1f--.???  KddAaDDIIaLL1a')*1Q44"tAwwtBxx/T!WWQYq[0A0AAA31uQw<<OPPRTUVWXUYUY[\R]R]^999K$ttyy||AaDDAI%! !Q!AaDD&(@AAaDD"VQBZZ-"VQBZZ-788AaDD"S!Q$q\\/9R9RS%K.ttyy||AaDDAI%! !Q!A$$(>?1Q44CAaLL))2addBs1vvI3q66 4I+J+JK/K8ttyy||Qq 1Q4499Q<<?1#451Q44a  "QqTT26"2"2BqttR!A#X4F4F1Q44QSVWXYVYQYHZHZ[9KBttyyAs1Q44yy(!A$$))A,,6 qttS"c$s))R!V"4559:;; < < qttS"c$s))R!V"4559:;; < < qttRLL CKVttyy||Qq 1Q44! ! #4Q#671Q44a  "QqTTDacNN?";";BqttT"qQRs(^^?? @ @ qttSSDBKKa!899A=>?? @ @ qttSQRRU ^R!V%?? @ @ qttSQRRU ^R!V%?? @ @     YKp!A$$))A,,!a'(3q!t99q!A$$i@@"q&IIJqK|ddiillAo!a'1Q44BEAbDFNQT12334}KHddiillAo!$1Q44Qq2a46)AqDF2334IKRddiillAo!a'1Q44CFF##R!baRj%9%9:SK\ddiillAo!a'1Q44BEAbDFNQT12334]KhddiillBq!tGbAg-1a47"Q$>CtLL1Q44hqb2a4i"QT'.ABq!tG.KbQRTUQUg.U)V)V'VXYZZZ[[\iKx1q53z!A$$//00111Q44hrD1Q3LL'8!<<<==r!A$$XVZ[\]^_`]`[a[aVbVbdeMfMfHf?g?gh 1q!!1$q!t+1Q44a8Aa== 002233R!b1QQRST EUCUVWCW>W5X5XY ! !Q!AaDD&(1Q44CGGbQi/001OKKKO O OrVc ddlm}tdd}td}|tt}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_05r)'sympy.functions.special.error_functionsr"r)r(r<rr rrrrrrr)r"rr?ras rTrJrJ)s ::::::sT"""A A ! ! A$AaDDd11Q447l1nqqttQw&*11Q44"S1XX+/HQNN:;; 11Q44"tAwwqy.2c!A#hh;?Xa^^"CCAEFF 11Q44"tAwwqy.2c!A#hh;?Xa^^"CCAEFFH!fQh1QQqTT1W$q(11Q44$Qrr!A#wwYBqD 9:::;;R!dBQRSUSUVWXYVYSZSZQZN\_`bcd`d\e\eKeFfFf=g=gh sQqqttV|a11Q44"rr!uu*c1"gg-../ AAaDD k!mc!!A$$ii'!A#-11Q44RTAhq!nn$4455667 !QaaddA&1-..aaddA>1EI11Q44"y Bq1uaLLaRTAX\A-s1vv5q1u=tD F FFGJA2wwOPPQ1! ! !rVctd}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 t ttt ttttzzz gdt ttttdttz t tttdzzzdtz t tttzzt tttttt tzz t tttttt tzzgdtt ttttztttz t tttdzzztt tttzzt ttttttt zz t ttttttt zzgdt ttttzttdt 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) Nrqr>rMrrrrr) liouville_01 liouville_02 liouville_03 liouville_04 liouville_05 liouville_06 liouville_07 liouville_08r)r)r r?r<rr rrrrexpandrr.rr/)rq_ys rTrNrNTsps A sBaDD1Q44mmAoQqTT1a 0 0 22T!A$$]]A5Ea5GG1Q44QRT ]++,,-1S!A$$ZZ<A&&1Q44QRT ]++,,- addAq4!a#3#3A#55QqTT1 q8H8JJCQRSTQUQUPUJJVWZ[\]^[_[_W`W``hhjj1Q44QRT ]++,,- 1Q44A1Q44addA(:!::QqSaddA=NN1Q44$rBs1vvI~...//AaDD$rBs1vvI~:N:N1O1OP QqTT1a   1QqTT6$qttQ--*:#::QtAaDD!}}_L1Q44$rBsA2wwJ///00"QqTT4RQBZ;P;P2Q2QR !C!II+##Aq))1--1Q44R"Q$Y(() AaDD!}}QaddAq!1!1!!33d1Q44mmQ6Fq6HH#qQRtte**TUXYZ[\Y]Y]U^U^^1Q44QRT ]++,,- d4!Q<<1QqTT6AaDD!G#3T!A$$q\\1_"DDtAaDDSTVWRXGYGYY2:hqmmQT!V<<<>1Q44bQA%6q%8%:!;!;;;bQQRTUEVWXEXYZEZA[A[>[[\\]ddii1oo!A$$))A,,.1Q44bRTl*++, ddii1oo! !Q/!AaDDIIaLL.@1Q44bRTl*RAaC[899:ddii1oo! !Q/!AaDDIIaOO2CC !A$$))A,,1Q44bRTl*RAY6C!HHDEEF ddii1oo!A$$))Aq// 111Q44bdR1W,r#bd))|;<<=ddii1oo!A$$))Aq// 111Q44bdRT!WWHQJ%77"Sa^^:KKLLMddii1oo!A$$))Aq// 111Q44bdRAaC[02c!A#hh;>??@ddii1oo!A$$))Aq// 11AaDDIIaLL@1Q44bd3q66k)Bs1vvI51SVV CbQiORTUXYZU[U[R[[^`adefagag^gghhiddii1oo! ! a/1Q44d1ggd2rAv;&7&77a@@AA 1Q44d1ggd2rAv;///a8899; 1Q44Aq))JqttQ1,E,EE1Q44bQA!a5!I!I!III 23q1a4!A#:>U;;;; ;<=>>? MLLP P PrVcVNtt }tttd}dtttdzd|zztttzttz 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|dtttzzdttzzdz tttttdtzzttt zzdzgddd|dtttzzdttzzdttzz tttttdtzzttt zzdttzzgddd|dtttzzdttzztttzz ttttddz tdz zttdtzzttt zzttttzzz zgddd|dtttzzdttzzttz tttttdtzzttt zzttdz zdttzdz z gddd|dtttzzdttzzttz tttttdtzzttt zzdttzdz zttdz zgddd|dtttzzdttzzddttzzdttzzz tttttdtzzttt zzttzttdz zdttzdz z dzgddd|tttzttztdzz tttdtztdzzttttdzdz zttttdzdz zztt dz zzgddd|dtttzz dttzz d tzttzz dtt zz tttt ttzdtt zz ttdtzzzttdtzzzgddd!|dtttzz dtdtzzttzz ttttttdtzzzdtdtzzttzdz z tdtzttzdz z gddd"tttdd|zz ttztz ttzttttttdz z tttzztt zztt tzzttzzgddd#|tttztdzz dtzz tttttdzdz zttt zzgddd$|tttztz tdtzz tttttz tdtzdz z tdtzdz z tdzdz zttt zzgddd%|ttzdtzttzz tttttdzz ttzttzttzzgddd&|dttzzttdtzzz tttttdzdz z tdtzzttd'z ztdtzzzgddd(|dtttzzttztdztt zz tttttttdzdz zzztt zgddid)tttdd|zzdtttzzttzdtt zz tdztt zztttttttttdzd*z z tdz zzzzztt zgddd+|dtttzzdttzztdtzz tdzz tttttt ztdzdz ztt#ddzz ttz tdtzzzt#d,dzgddd-|dtttzz dttzzttt zz ttttttzttdtzzzdtzdztt zd.z zgddd/|tttzdttzz tz tdtzz tttt#d0d.tdz z ttd1tzzzttdz ztdtzzzgddd2|ttzttz tt z ttttttzttdz z ttzztt dz zgddd3tttdd|zz dtttzzttz ttz tttttttttdz zzzzzttzgddd4|ttztjz tdtzdz z ttttjtdtzdz z tttzztttzzgddd5tttdtttz tdtztjtdtzdz z zz ttttttt zztttzzd6tdtzzd7tdtzzzd8ztdtzzd9z zgddd:tttddtttdzztttzdtzz ttz ttz tttttdzzttttdz z zzttzzt tt&tdz zzzttzzgddd;|ttzttdz z tdtzdz zttttdtzd'z tttzzttdz zttzzgdddtttdtttz t)tz ttttttt zztttzztt)tzdz zgd?d@tttddAttzz t)dtzz tttttdBtzzttd,tzzzt)dtzdCz z gd?dDtttdtttz t)tz ttz ttttttt zztt)tzdz zttdz zttzzgd?dEttttttztttdz zz ttttttztttzztdzttdz zdz z tttdz zdz zgd?dFtttdttztdzttdz zz ttttttztttzztdzttdz zdz z tdzttdz zdz ztdzttdz zdz zdtdzzttdz zdz z dtzttdz zdz z gd?tttdttz ttdz z tttttt ztttd0zdz zttzzgd?tttdttzttdz dzz ttttttztttzztttdz zdz z dzgd?tttddtttdzztttzdtzz tttzz tttttdzttttzzdz dzzztttzzttzzt t&tzzttzzgd?t||||||z |||zzd||zz ||zz| t||||z| ztt| |t||dzzd|zz t|z zzd|zz zztt|| t||dzzd|zz t|z zzd|zz zzg||dGt||z||||z||z|||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)JNrrrrWr=)clszR L C E_0 alphaTromegar.undet_01rrrrundet_02reundet_03rrundet_04rfundet_05rundet_06undet_07undet_08rsundet_09undet_10ruundet_11undet_12undet_13undet_14undet_15undet_16rundet_17undet_18rundet_19rundet_20rUundet_21rrrundet_22undet_23undet_24undet_25iiundet_26undet_27undet_28undet_29undet_30rundet_311irundet_32undet_33undet_34r)undet_35undet_36undet_37undet_38undet_39undet_40r)rr<r?rrr(r)r rrrrr r rrrrr/rrr) r@f2rYrWr=RLCE_0alphar.s rTr'r'vsc QBA 11aB !A$$))Aq//AbD 1Q4499Q<</!A$$6:A A!!!A!"3TBBBAq!S% 7OOEIaDDv!A#g1Q44C!HHqBBAbL1Q3r64I1J,JCQRPRGG+SSVWWXXYv!e1Q44C!HHqBB1H,=sA2ww+FFJKKLv1QqTTYYq\\>!AaddF*Q.1Q44C1II3r77 2Q6778v&1QqTTYYq\\>!AaddF*RAY61Q44C1II3r77 2Qs1vvX=>>?'v21QqTTYYq\\>!AaddF*S1XX51Q44!A$$r'AbD.2c"Q$ii<"S!WW*+DqQqSz+QRSST3v>1QqTTYYq\\>!AaddF*SVV31Q44C1II3r77 2SVVBY>3q66"LMMN?vJ1QqTTYYq\\>!AaddF*SVV31Q44C1II3r77 2Qs1vvXb[@3q66"9LMMNKvV1QqTTYYq\\>!AaddF*a!CFF(lQs1vvX.EF1Q44C1II3r77 2SVV;c!ffQhF3q66RSSVWWXXYWvb1Q4499Q<<!A$$&A-1Q44A13qay{+;+;(;bQtAwwYq[AQAQ>Q(QSVXYWYZ[W[S\S\'\\]]^cvn1QqTTYYq\\>!AaddF*QqSQZ7"S!WW*D1Q44!CFFQsA2wwY.C1II=3qs88 KLLMovz1QqTTYYq\\>!Ac!A#hhJs1vv$551Q44bQqSk)Ac!A#hhJs1vv,=a,??#ac((3q66/RSBSSTTU{vFddii1oo"$qtt+a/#a&&81Q44SVVAXbdCGG(;;rBqDy#a&&>PPQQRGvR1Q4499Q<<!Q$&1,1Q44ad1fr#qb''z1223Sv^1Q4499Q<<!#c!A#hh.1Q44a#ac((1*,s1Q3xx{:QT!VCbaRjPQQR_vj1Q44i!A#c!ff*$1Q44"q!t)SVV+rAvs1vvo=>>?kvv1QqTT6kAc!A#hhJ&1Q44"q!tAv+s1Q3xx/2"9c!A#hh2FFGGHwvB1QqTTYYq\\>!AaDD(1a4QB<71Q44"q"q!tBw,//aR899:CvvNddii1oo"$q11~5!!AaddF*QsA2wwY61Q44CFFRAaC[0AaC!GS!WW3DR3GGHHIivt1Q4499Q<<!AaDD&(1,s1Q3xx71Q44"b))AaC/"SAYY,>"qs(CPQRSPSHHATTUUVuv@1Q44i#a&& 3r77*1Q44CFFb1Q3hA%66aRBCCDAvLddii1oo"$q11~5!MvX1Q44i!& 3qs88A:-1Q44#ac((1*,r#a&&y82c!ff9DEEFYvdddii1oo! ! ,s1Q3xx#ac((1*9L/MM1Q44baRj2c!ff94C!HH r#aPQc(({8RUX8XZ]^_`a^aZbZb7bcg7gghhievrttyyA1Q4499Q??!22QqTTYYq\\AAaCG AQ1Q44adb1b1Q3h<&7Q%??22PQRSPS8 CTVYZ[V\V\B\\]]^ sv@1Q44i#a&&("S1XXaZ/1Q44QqS"r#a&&y0B1Hc!ff3DDEEFAvLddiillQ1Q44a  !MvZ1Q4499Q<<#ad))+1Q44baRj3sQw<<>9::;[vhddii1oo! ! ,tAww61Q44baRj2c!ff94qay{BCCDivrddii1oo1Q44'$qs))31Q44C1II3qs88 3d1Q3iilBCCDsv|ddii1oo! ! ,tAww6Q?1Q44baRj1T!WW9Q;6"qs(CFF9JJKKL}vHddii1oo!$qQU|31Q44CFFRAY.Ac!a%jj1BBQs1q5zz\RS^STTUIvRddii1oo!$q!tC!HH}4AaDD"SVV)bQi/!Q$s1q5zz/"2DDq!tCPQTUPUJJWXGXX[\^_[_`cdehidi`j`j[jkl[llopqrtuquouvyz{~zwAwApABCpCCFGHIFIJMNORSNSJTJTFTUVFVVWWXSvv^ddii1oo!$s1q5zz11Q44CGGrAc"ggIaK'7Q&??@@A ddii1ooadd"C!HHQJ/1Q44CFFRAY.3q1u::a?!CDDEddii1oo!A$$))Aq// 11AaDDIIaLL@1Q3FQqSQ1Q44adAc!A#hhJqL1$455bdCFF8JJbSUVWSWiY\]^Y_Y_M__``a!!A$$))Aa..1a4! ! #44q!A#wqqtt|CUKK11Q441URQBD1a4!A#4F4FtAww4N0N,OQRSTQT,U(V(V%VV SQBa1fqsl++DGG334ac: ; ; ;<==>! !A#aaddii!nn$qs11Q4499Q<<'77!!A$$>C%PQ NN@RTT11Q44CAQq!tVac\(:(:477(B$B CQqS IJJJ SQBa1fqsl++DGG334ac: ; ; ;< c!E'!)nn ac%(lQqSU5[81< =>??@! AaDD!Q''!JqttQ,?,?*??!A$$FAOO1Q44"r!t)SVV+c!ffQh6778evvvz z zrVc-td\}}d}d}d}||z}td}dttidtttttz t ttt ttzgdd ttttzttz t ttt tzgdd tttttzt ttt ttzgdd ttd zd ztd zd ztttzz t tttt ttzgddtttttz ttz d z t ttt ttz d z gddtttd tttz zd z t t ttztttzt gddtttd zztttzttdzz d td zztttzz 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 ttztttzztttzz t tttt td zzdt ztzzdt zzd z t tttt td zzdt ztzzdt zzd zgdddtttztttztd ttd zzzt tttt tttz gddddttd ztttzttttttzzt tttzttttd z d z zttttd zd z z tttzt gdddttttztd ztttztttzz|d ztttztttzz t ttt!t t|d z td zzz d t"zzt ttt!t t|d z td zzzgdddtttttttzz t ttt ttz gddtd z tttztttzd tztttzz t ttt"t%t td zd tzz d zztd tzzz t ttt%t td zd tzz d zztd tzzgddttttttttzttz z t tttt ttzgddttttzd ttd zztttzzt tttt tz gddddttttttd zzzt ttdt ttd z d z zzgdd tttd ztd zd tzzd zztttzttztttzzt ttttd t td zzd tzzz  t ttttd t td zzd tzzz gdtttttzt ttt tt zgdtttd ttzztttd ttzztttzzt ttt"t!t ttd zz d z z t ttt!t ttd zz d z gdd tz tttztttd zzz t ttt tt ztz d ztd z z gdttt tttztzdtzttd zzz t tttd tt tdtd zzzd zzdz t tttdtt tdtd zzzd zzdz gdtttttttzz t tttdt ttzz gd!gd"ttttzd zttd zz t ttt td zz t td zzz gdt||d d#t&zt|zz z t t|d$tt&t d#|zz zzd$zt&z gt|d%tttt ttzz t ttt(tt tzzgd||||d zzz |t+||zz t ||d&t-t d'|zz z g||dd(tttttttz 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_17rrwgfffffI@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_27r)rr(r?r<rr rrr rr!rrrrrrr rr>rrr)rWrrr@rf1rbs rTr+r+~sK %..CAa A A A QB AaDDcddiillQqTT!1Q44CFF##$c ! ! nqtt#1Q44A cddiillSVV#1Q44c!ff%%&c ddAgkQTAXqttyy||331Q44R$q''\**++,!c*ddiill3q66!AaDD(1,1Q44CFFQ''(+c4ddiilla#add))m,q0A2!9s1Q44yy("--.5c>dd1a4i! ! $qttQw.1a4! ! 1DD1Q441"tAqtAvzA~$6777"Q$(CAEFF1Q441"tAqtAvzA~$6777"Q$(CAEFFH ?cLddAgkQqttVa!f_addiill::1Q44$r!Q$w2a/!B$6:;;;<< 1Q44bAg"Q&2-122335 McZAhqttyy||#d1qttQw;&7&771Q44b3s1vv;;.//001! [ch!a%jjQqTT"S1YYqttyy||%;;1SVV8c#add))a-0022SQqTTQ5G5G5IICPQRSPTPTIIUWYZZ[ictQqTT{QT#add))^AaDDIIaLL881a4AaDD >!A$$))TU,,;VV qttd2dAqD51a4<000111AbD8 9 9 qttT"T1a4%!Q$,///00 1 1  ucFddiillQqTT#a&&[(1Q44CFF##$GcP1uc!A$$ii! ! ,qs3qtt99}< qttR$r1a4!A#:>23qs88;<<< = = qttT"adQqSj1n-c!A#hh677 8 8 Qc`ddiillQqTT#add))^CFF221Q44RAY(()acj! ! nAaDD!G T!A$$ZZ771Q44RT##$! kcxddiillQ!q\)1Q44badU1Wo-../ycB!A$$'llAqD1Q3JN+q1v!}addiill.JJ qttd3q"q!t)ac/233444 5 5 qttT#aadQqS12233 4 4 CcTddiillQqTT!1Q44CGG$$% !ffS1Q44[[ 3q66#a!f++#5addiill#BB1Q44d2c!ffai<002233R!d2c!ffai<>P>PQR>R5S5ST 1uaddiill"Q!q\11Q44"S!WW*q.1,q1u5667 dd4!a== 1$qs1Q447{21Q44$q''$r#a1f++~'9":"::1<== 1Q44abQq!tVnq0111!3446ddiillSQqTT]]*1Q44Rc!ff-..//011"QqTT1W,1Q442#1*sQTz2334ddiillAeAgaddlN+1Q44.s1b;Ma;O6O3P/Q/QQThhjkkllm!ddiillR!W$1Q44C1II&&' 1! ?"QAAaDD)9)9%9911Q44,T"7J17L2L-M-MMNNO!! ddiillS1]]*1Q44c!ffWbQi!m455667cccg g grVcldtttttttztttttztttzzt ttt tttz  dtzzt ttt tttz gdddtzttzdzttz tttz ttdzz tttzzt tttttdzz tt tt tdzzzzgdgddddtzttttzzdttzttztttz tttzzt ttttztttztdzzttdzztgdgddtttttttzttdzz tttzz t ttttzttd zd z ztgdddtzttztdzttdzztttzzt tdzttzttd zd z ztgdd d tttttttzttdzz tttzz t ttttzttd zd z ztgd d tttdzttdzzztdzttzttttdzttdzzz z tttzz 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 tdztttzd tzttzztttz z dt tttttz td zz gdtttdzz tttzdz tz tttzzddDdtttzdz ttttdztttzz zzt ttttttdzzdttzz z t ttttttdzzdttzz zgdtdztttzttttztttzzt tttttt ztdzz  tzt tttttt ztdzz gdd dS)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;rIrErFrIruirH)rErFr;rKrrc g|]R}ttt|tttdzzdzzdztz SS)rr)r r?r<rrrs rTrz3_get_examples_ode_sol_1st_exact..sKGGG!1Q44!DAqD1---11455GGGrV)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_11r)r?r<r rrr rrr rrrrrrrWrVrTrr2sH  aDD!ffS1YYQAaDD !1!A$$))A,,!>>1Q44$r#a&&y//)AbD0112addDCFFOO3L3LMs1Q44x!|QqTT!QqTTAXqttQw$6qttyy||$CC1Q44R!Q$Y1"S"q!t__2D)E)EEFFGGH cAaDDQK1QqTT6CFF?S1YY#>! ! "LL1Q44A;QqTT*QT1AaDD!G;R@@A !A$$ii1S1YY;1q0!A$$))A,,>>1S1YY;1q*B//0c!A$$h!Q$1q.!A$$))A,,661a4!9qttQwqy("--. !A$$ii1S1YY;1q0!A$$))A,,>>1S1YY;1q*B//0QTAaDD!G^$$$1QqTT 1Q44$q!taddAg~:N:N3N(OQRSTQUQUQZQZ[\Q]Q]']]3q66 QtA!a1 $%% %eAaDDFmm 3S1XaZ 4AaDD!GAqDL ! !!6"# #%&tA!a1 ,<'='=%= AQ1qA%&&qtt+A- -!A$$' !Q$ >??&@ QqTT!ba!A$$'!Q$,./// /&1 1 %!Q--!aQqTT!ba!A$$'!Q$,6F1G1G.G!GHI  J !A$$s1tA!a1 ,--add2144q1qyA~EFFF C!HQJDQqTT1WQT\!1222 2 3 5  5667   Aaddiill"QqS1X-Aq8!<<1Q44"s1vv+q!t+,,-dd1a4iAaDDFQJ>1Q4499Q<<77GGwGGG 1vzQqTTYYq\\1a4!AaDD&=991Q44T"q!t)aAh"677788"QqTT1tBQRTUQUIXYZ]^_Z`Z`X`L`GaGaCa:b:bc 1uc!A$$ii!A$$))A,,q.QqTT"::1Q44$r#qb''z!Q$///"455r!A$$RQBZPQSTPT_@U@U7V7VWSMMQ Q QrVc!tt }tttd}dtttdzd|zztttzttz 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 tttdz t ttt tzgdd |dtttzzdttzzd z t ttt td tzzttt zzdzgdd |dtttzzdttzzd ttzz t ttt td tzzttt zzdttzzgdd |dtttzz d ttzz dtzttzz dtt zz t ttt ttzdtt zz t td tzzzttd tzzzgdd |dtttzzttztdztt zz t ttt tttdzd z zzztt zgdd |dtttzz dttzzttt zz t ttt ttzttdtzzzdtzdztt zdz zgdd tttdd|zz dtttzzttz ttz t ttt tttt tdz zzzzzttzgdd |dtttzzttztt tz z t ttt ttttzzztt zgdd |ttzdttz dzttz z t ttt tttdz dz ztttdzdz z ttzttttdz dz ztttdzdz z ttzzgdd tttd dtz z t ttt ttzzt tdzzztdztttdz zzzgdd tttddtttdzztttzdtzz tttzz t ttt tdzzttt td z z dtttzzdz zdtt tzzdz zzzdtdtztzzdz zdtztttzzzdtztt tzzz ttzztttttzd z zdtztttzzdz zdtztt tzzdz z zzdtztdtztzzdz zdtttzzz dtt tzzz ttzzttttzzz gdtttddtttdzztttztttzz t ttt ttt td z z zzdtdtztzzdz zttzztttttzd z zzzdtztdtztzzdz zttzzttttzzz gdtttttdttzzttz t ttt tdz z tdtzdz ztdtzzttttzdtttdzzdz z dttdzzdz ztdtzzzgdd ddS)NrrrrrrrTrrerrrfrurrsrUrr)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_15r)rr<r?rr rrrrr rrrrrr)r@rZrYs rTr r se  QBA 11aB !A$$))Aq//AbD 1Q4499Q<</!A$$6:AGaDD!A#g1Q44C!HHqBBAbL1Q3r64I1J,JCQRPRGG+SSVWWXXY!e1Q44C!HHqBB1H,=sA2ww+FFJKKLddiillQ1Q44a  !1QqTTYYq\\>!AaddF*Q.1Q44C1II3r77 2Q67781QqTTYYq\\>!AaddF*RAY61Q44C1II3r77 2Qs1vvX=>>?1QqTTYYq\\>!AaddF*QqSQZ7"S!WW*D1Q44!CFFQsA2wwY.C1II=3qs88 KLLM1QqTTYYq\\>!AaDD(1a4QB<71Q44"q"q!tBw,//aR899:1QqTTYYq\\>!AaddF*QsA2wwY61Q44CFFRAaC[0AaC!GS!WW3DR3GGHHIddii1oo"$q11~5!1QqTTYYq\\>!AaDD(3r771941Q44"q"s1vv+.QB78891Q44i!CFF(1*SVV++1Q44"s3q66A:q003s1vvz??13DD a&&SVVaZ**SQ!__Q->>AFGHHI ddii1oo!#1Q44bdR1W,q!tR#a&&(]/CCDDEddii1oo!A$$))Aq// 11AaDDIIaLL@1Q3FQqSQ1Q44adb1b1Q3h3qs88A.E#qbQRd)) TU .U+V&VYZ[^_`ab_bcd_d[e[eYefhYh&hklmnknorstuvsvowowkw&wz{|}z}BDECEFGCGHH{H'HJMNOJPJP&PPTVYZ\^abcdadefaf\fijklilmpqrstqtmumuiuvwiw\wz{|}z}~ACDBDEFBF~G~GzGHIzI\IYJTJ A#c!A#a%jj. TAaCjT)+,S!AYY;T78;AS??AB3qs88LMMNddii1oo!A$$))Aq// 11AaDDIIaLL@3qs88K1Q44rArAaCxL01S1QZZ<?BCFFJJbSTVX[\]^[^_`[`V`SaNadefgdghklmnolopqlqhrhrdrsuduNuwz{|w}w}M}}ABCFGHIJGJCKCKAKKLLM1q!$$q1v-A61Q44"qs(S1XXaZ/QqS9R#c!ff++=M CA NN 1 > Q{1}>-.1!A#hh=77889 o]]a a arVc&dtttdztttdzttttzztdzdz ttzzttttt dtzt tdtzzgdtdztttdzttttzztdzdzttzzttttt dtztzt tdtztzzgdtdztttdzttttzztdzttzzttttt dtzt tdtzzgdtdztttdzttttzzdtdzztd d z z ttzzttttt td d z d tzzt ttd d z d tzzzgdtdztttdzttttzztdzdz ttzzttttt d tdzdz zt td tdzdz zzgdtdztttdzdtztttzztdzdz ttzzttttt td dz tdzdz zt ttd dz tdzdz zzttz gdtdztttdzttttzztdztd dz z ttzzttttt td dz tzt ttd dz tzzgdtdztttdzd tztttzz dtzdzttzzttttdztt ddttzzt tddttzzzzgdttttdztttz dtd zzttzztttttt td dz tdzzt ttd dz tdzzzzgdtdz dztttdztdz tttzz dtdz dzzttzzttttdz tt d dtzdz zt td dtzdz zzzgdttttdtz tttzzttzttttt td dz tzt ttd dz tzzttz gdtdztttdzttttzztdztdzztdzz tdzz ttzzttttt ttdztttdztdzz zzt tttdztttdztdzz zzzgdd dS)N2nd_linear_besselrrrrrQrrur) 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_12r)r?r<rr rr,rr-rr rrrYrXrWrVrTrrs(aDDdAaDDIIaOO$q!A$$))A,,'771a4!8QqTT/I1Q44GAqMM)Bwq!}},<<==> dAaDDIIaOO$q!A$$))A,,'771a48QqTT/I1Q44GAaCOO+b1a.@@AAB dAaDDIIaOO$q!A$$))A,,'771a41+E1Q44GAqMM)Bwq!}},<<==> dAaDDIIaOO$q!A$$))A,,'772ad7AaDDF?AaDD:PP1Q44GAaDDFAaC0002gadd1fac6J6J3JJKKL dAaDDIIaOO$q!A$$))A,,'771a4!8QqTT/I1Q44GAq!tAv...GAq!tAv4F4F1FFGGH dAaDDIIaOO$qsAaDDIIaLL'99QTAXqttOK1Q44"WT"XXaZAa8882gd2hhqjRSUVRVWXRX>Y>Y;YY[_`a[b[bbccd dAaDDIIaOO$q!A$$))A,,'771a4!A$$q&=!A$$:NN1Q44GAaDDFA...GAaDDFA4F4F1FFGGH dAaDDIIaOO$qsAaDDIIaLL'99QqS1WaddNJ1Q44Ar'!QtAwwY"7"77"WQ$q'' =R=R:RRSTTU 11a!AaDDIIaLL01QT6!A$$;>1Q44BwqttAvq!t444r'!A$$q&!Q$:O:O7OOPQQR sQh! !Q(AaC11+==1Q3( 1Q44O1Q44!a%"WQ!a%8%8"82ga1q>Q>Q;Q"QRSSTddii1oo!AaDDIIaLL 001Q4471Q44"WQqTT!VQ///"WQqTT!VQ5G5G2GGaPQQR d1Q4499Q??"Qqttyy||^3q!tAqDyA~17LaPQdd6RR1Q44GDAJJ$q!tAqDy//0ABBBRPTUVXYUYPZPZ\]^bcdfgcghiklhlcl^m^m\mHnHnEnnoopq<<@ @ @rVc dttttdz ztttdztddz dtzz tttzzdttzzt ttt ttddz zzt tddz tddz ftddz ftztt dtd dz ftzzgdttdz ztttdztddz tztttzzttzt ttt dtz tddz zzt tddz dftddz fdtz ztt td dz d ftd dz fdtz zztdz tddz zz gdttdz ztttdztdtddz tzztttzzttzt 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)rrrrrrTrui`iQ rPi*rs-brrkTTrx)2nd_2F1_hyper_012nd_2F1_hyper_022nd_2F1_hyper_032nd_2F1_hyper_04r) r?r<rr r rr$rrrWrVrTrrAs{)aDD1gaddii1oo%1a1addiill(BBQqttVK1Q44A!QKqttAvqttAv.>1a 1(M(MMPRSXYaefgheieidijkdkcmopSqSqPqqrrs 1gaddii1oo%1a1Q4499Q<<(??!A$$F1Q44"a!eqttAv..uadd1fa[1Q446)QQRU/S/SS UQqTTE!GR=AaDD57*a!e 4 4 45781u!Q6GHIIJ1gaddii1oo%1qttAvax11(EE!L1Q44"a!equuQw//qttAvqkAbEE!G:qSTu0U0UU UQqTTE!GR=AaDD57*a!e 4 4 45781u"a6HIJJK1Q446{lDQqTT!V,Q.a!A$$q&k1A"1DDqvvaxPQzQRSTUVSWSWWY\_`cdefcgcghici_j^j _ ^ZAaDDF a!#Z$%'1!A$$A'?'?@1Q44QrUU2X1qttBw<qttAv!uQw>ORSTURVRVWXRXQZ\]`abc`d`def`f\g8h8h(h UAaDDGaeeVBY'!A$$q&A!QK @ @ @)ABCDqttAv;QR?VWXZV[V[\^V^B_`aab %  rVcdtttttdtttzt ttzt t dttdtzt tzz dtzzdttzz z tttfttzt t dttdtzt tzz dtzzdttzz z tttfttz gddtttdtttzdttz zt tdt dtdtzdttddzzz dt tzz z tttfzttzt tdt dtdtzdttddzzz dt tzz z tttfzttz gddtttdtttzt t dttdttzzz tttfttzt t dttdttzzz tttfttz gddtttdtttzt t dttdt!tzz z tttfttzt t dttdt!tzz z tttfttz gddtttdt#ttzt t dttdtzt#tzz dtdtdzz zz z tttfttzt t dttdtzt#tzz dtdtdzz zz z tttfttz gddgd d d S) N"2nd_nonlinear_autonomous_conservedrrFrrrrrErFrIrG)%2nd_nonlinear_autonomous_conserved_01%2nd_nonlinear_autonomous_conserved_02%2nd_nonlinear_autonomous_conserved_03%2nd_nonlinear_autonomous_conserved_04%2nd_nonlinear_autonomous_conserved_05r)r?r<rrrr r%rrrrrrr rrrrrWrVrTrrbs9aDDddii1ooAaDD )C!II5 x$rAbDRL01R47!CGG)CDDDr1Q44jQQSUXYSY Z Z x$rAbDRL01R47!CGG)CDDDr1Q44jQQSUXYSY Z Z  ..ddii1ooQqTT *QqttV3 tAwwx$qtaHQNN0B.B'BQs2wwY'N"O"O ORTVWXYVZVZQ[\\\^`cd^d e e tAwwx$qtaHQNN0B.B'BQs2wwY'N"O"O ORTVWXYVZVZQ[\\\^`cd^d e e  ..ddii1ooAaDD ) x$rAc"ggI~...QqTT ;;R!V D D x$rAc"ggI~...QqTT ;;R!V D D  ..ddii1ooQqTT * x$rAd2hhJ///"add<>!#c!Q$qtt)hq!nn*D&E&Ea&GGcRSffUUV ddiilladdad1Q44i!m453q!tAaDDy>>!#c!Q$qtt)hqmm*C&D&DR&GGcRSffUUV dQqTT1a4!9%%3q!tAaDDy>>!#c!Q$qtt)a-&8&8&::BQKHHI11!Q!q1Q446{1Q44/?!.C*C DDaHH2c!add1f+add*+++aAa!A$$q&kM!A$$'4IA4M0N0N.Nq.PPRTWZ[\W]W]R]^^_11!Q!q1vk/ ::A>>1Q44$q''$q"s1vv+"7"771=>>addDGGDBQK$9$991Q3?@@B11AaddAg1QqTT !91Q44 ?AqD1qLSTVWSWXYZ[X\X\^_X_S_D_A` aacdee1Q44a1a4j))* 111q!A$$ 22A66 qttd1ggXd1b3q66k?333A5 6 6 qttT!WWT!R#a&&[/22214 5 5 111a1 q!A$$q&z ::A>>C!qMM>!#aAaDD k!m3R#a&&[AAB 111a1 a/331Q44Bq!tGaK))* 111q1qtt 3Q 77;;SVV8c!A$$(mm+c!A$$ii7!AaDD&@"EEFddiilladdad1Q44i!m451Q44$q''$qHQqMM)9'943r"u::PQTWPXAXZ]^`ac^cZdZd@d;e;ehklnoqlqhrhrstvwswhw;w{CDEFG{H{H;H(HAx!}}SBZZqRU|a#h'>BrE &J!K!KcRTUWRWjjYZ\]Y]o!]`hijkl`m`m mm(npqrsuvrvpv(w#x#xxxyzq''$r!Xa]]""D!C2JJ,S*A3r"u::)M$N$NQTUWXZUZQ[Q[\]_`\`Q`$`cklmnocpcp#ppAx!}}SBZZqRU|a#h'>BrE &J!K!KcRTUWRWjjYZ\]Y]o!]`hijkl`m`m mmnpqrsuvrvpvwDGG)QT$qHQqMM))43r"u::C1H#bQSe**0T+U+UX[\^_a\aXbXbcdfgcgXg+gjrstuvjwjw*wwAx!}}SBZZqRU|a#h'>BrE &J!K!KcRTUWRWjjYZ\]Y]o!]`hijkl`m`m mmnpqrsuvrvpvwxxxyy z zz{||| @BCDEGHDHBHI IJJ 1Q44$q''$qHQqMM))43r"u::C1H#bQSe**0T+U+UX[\^_a\aXbXbcdfgcgXg+gjrstuvjwjw*wwAx!}}SBZZqRU|a#h'>BrE &J!K!KcRTUWRWjjYZ\]Y]o!]`hijkl`m`m mmnpqrsuvrvpvwxx xxy zq''$r!Xa]]""D!C2JJ,S*A3r"u::)M$N$NQTUWXZUZQ[Q[\]_`\`Q`$`cklmnocpcp#ppAx!}}SBZZqRU|a#h'>BrE &J!K!KcRTUWRWjjYZ\]Y]o!]`hijkl`m`m mmnpqrsuvrvpvwDGG)QT$qHQqMM))43r"u::C1H#bQSe**0T+U+UX[\^_a\aXbXbcdfgcgXg+gjrstuvjwjw*wwAx!}}SBZZqRU|a#h'>BrE &J!K!KcRTUWRWjjYZ\]Y]o!]`hijkl`m`m mmnpqrsuvrvpvwxxxyy z zz{|| | @BCDEGHDHBHI  IJJ 1Q44aa8Aa==(($#be** q3x0GRPRU/S*T*TWZ[]^`[`WaWabcefbfWf*fiqrstuiviv)vvAx!}}SBZZqRU|a#h'>BrE &J!K!KcRTUWRWjjYZ\]Y]o!]`hijkl`m`m mmnpqrsuvrvpvwxx xxy zq''$r!Xa]]""D!C2JJ,S*A3r"u::)M$N$NQTUWXZUZQ[Q[\]_`\`Q`$`cklmnocpcp#ppAx!}}SBZZqRU|a#h'>BrE &J!K!KcRTUWRWjjYZ\]Y]o!]`hijkl`m`m mmnAqD& T!WW9ad4!Xa]]"2 2D!C2JJ,S:QSVWYZ\W\S]S]9]4^4^adeghjejakaklmoplpap4ps{|}~tAtA4A!AAx!}}SBZZqRU|a#h'>BrE &J!K!KcRTUWRWjjYZ\]Y]o!]`hijkl`m`m mm!npqrsuvrvpv!wxxx yy z zz{|| | @BCDEGHDHBHI  IJJ 1Q44aa8Aa==(($#be** q3x0GRPRU/S*T*TWZ[]^`[`WaWabcefbfWf*fiqrstuiviv)vvAx!}}SBZZqRU|a#h'>BrE &J!K!KcRTUWRWjjYZ\]Y]o!]`hijkl`m`m mmnpqrsuvrvpvwxx xxy zq''$r!Xa]]""D!C2JJ,S*A3r"u::)M$N$NQTUWXZUZQ[Q[\]_`\`Q`$`cklmnocpcp#ppstuvyABCDEyFyFvFtFGJKMNPKPGQGQtQSWYZ[^_abd_d[e[eYec(Y 2JJXSSbe**QT/S*,4QqMMR:t::<=adFCEFtAwwYPQSTPTUYZ[\]_ghijk_l_l\lZlnrtuvyz|}zwAwAuADEHKDLuLNQRTUWRWNXNXtXoYoYbe**QT/oQqMMn"["$%a!A&6$6s2b5zz$A43rRTu::XY\_X`I`befhikfkblblHlCmCmpstvwytypzpz{|~{pCCKLMNOCPCPCP%P[PSTUVXYUYSY[YVZVZQZF[[ \ \\]^^ ^abdefgijfjdjak  kll%m.!5>d1Q447lQz!A$$22221Q442r!Q$w{+,,-Ceei i irVctd\}}}dttidttttzttdzzttdzzdttzz dtzz gddddttttzttdzzttdzzdttzz dtzz gdd d t td zt tttzd td zzttdzzzdtdzzdzttd zzz dgdd dttttttt tzz dztdzz z t tzgdgddtttdzt tttgdgddt ttttdzttzt tttttd ztdd zzgddttt|ttzz|t|tzzz t tttt ||zz|tt||zzzzt| tzz||zz t|| ft |tzzt| tzzdfgddtttdtzttzztttdz zz t ttttdzdz zttdz zgddddtzztttzdzdttt zz t ttt tdtzdzz dzgddtdztttzttz 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!tttttz tttttzzt tttttzt ttttt zgdttttttttz zt tttttzt tttgdtttttz tttttzzt ttttt zt tttttzgdttttttttzzt tttt ttttt zgdd"d#S)$NrAr lie_group_01rrT)rErFrKrL lie_group_02r lie_group_03rrrr lie_group_04rw 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_20r)rr?r<rr rrrrrrr rr r%)rrXrYs rTrr s, gGAq!aDDm! ! nadd1f%1q11QqTT69!A#= $ m! ! nadd1f%1q11QqTT69!A#=mAj1q)))AadF1Q447N:a1fqj!A$$PQ'=QQSTUUm*ddiilladdQs1vvXo1!Q$66Q?+m6ddiillAo1Q44 ~7mB11q!tAaDDy))1Q44C1IIx1~~55667CmLddiillQqttV#aAaCj01Q44bS!a%[1SAE^^3C%CS!AYY$NPQTUPU$V 1qb $cAaCiaRT*D133445MmXddiillQqS1X%#q!te** 41Q44"q!tAv+sAqD5zz1223Ymb1Q3w11&*QsAaDD5zz\91Q44R1q\A-..//0cmldAaDDIIaLL!AaDD(1a4A1I+>>1Q44"s1vv+s2a4yy0112mmxd1Q447lQz!A$$22221Q44BAI''(ymB1Q44llQqS1X%#q!te** 41Q44adUR!Q$q&[1223CmL1Q44llQqTT#a&&[(3qs8831Q44c!ffWrHS1XXc#a&&kk5I1,M,M'MNOOPMmV1Q44llQqTT#a&&[(3qs88A:51Q44CQLL3q661A5667Wm`QqTT! nqtt#aAh.1Q44"qQx-#a&&0!3445amjQqTT! nqtt#aAh.1Q44BSVV,-../kmtttyy||AaDD QqTTYYq\\!A$$%671Q44CFF##R!baRj%9%9:um@ddiilladdiillQqTT121Q44CFF##R!b\\2 ttyy||AaDD QqTTYYq\\!A$$%671Q44CGG$$b1r#a&&y&9&9: ddiillAaDDIIaLL1-.1Q44 b1r#qb''z223Smmmq q qrVc&dtttttdtttzz ttttt tzt ttzzgdtttddtzttzzttttt dtddz z tzzt tdtddz z tzzzgdddS)N2nd_linear_airyrrrr)2nd_lin_airy_012nd_lin_airy_02r) r?r<rr rr*rr+r rWrVrTrr{ s&aDDddii1oo!A$$&1Q44F1II6!99 4556 ddii1oo!AaDD(1Q44FA!QK<>222RAaDDF |A~8N8N5NNOOP    rVc4jtdd}tdd}dtdD\}}}}}dtdD\}}} } } d tdD\} } }}}d tdD\}}}}}d tdD\}}}}}td }d ttidtttddtttzzt tttttdtzzzgddtttddtttzz dttzzt ttttttzzttzgddtttdttz t ttttt ztttzzgddtttdtttdzdtttzz t tttttdtzzzttdtzzzgddddtttdzdtttzz dttzzt tttttdz zttttddzzzgdddt tttddtttzzttz dt tttttdtdzzzttt tdd zzzzgdddt tttdt tttdzd t tttzz dttzz t ttttdtzzttt dtdzzzztttdtdzzzzgddd!tttdtttdz dtttdzz dtttzzt tttttdtzzztttzzttdtzzzgddd"tttddtttdzztttdzdtttzz dttzz t ttttt ztttzztttd tzztttdtzzztdtzzzgddd#tttd|dzttzz t ttttt| tzzttt|tzzzttt|tzzzttt|tzzzgddd$tttdd|ztttzz dttzz t ttttt|t|dzdzz zzttt|t|dzdzzzzzgddd%tttdd|ztttzzd&|dzzttzz t ttttd'|ztzzttd|ztzzzgddd(tttdt tttttzzttdzzzttdzzzgddd)tttddtttzzdttzzt tttttzztdtzzgddd*dtttdzdtttdzztttzttz t tttttzztt ztttdz zzgddd+tttddtttdzz d&tttzzd,ttzz t tttttttzzzztdtzzgddd-tttdd|ztttzz |dzttzzt tttttzzt|tzzgddid.tttddtttdzzt tttttzzttdzzzttdtzzzgddd/tttddtttdzz t tttttzztttd tzzztttdtzzzgddd0tttddtttdzzdtttdzz d&tttzz d1ttzzt tttttzztdtzztttzztdtzzzgddd2d1tttdzd3tttdzz dtttzzdttzzt ttttt zttt dz zztttdz zzttttddzzzgddd4tttdd,tttdzz d5ttzzt tttttzztdtzztttzztdtzzzgddd6tttddtttzz dttzzt ttttdtzzttdtzzzttzgddd7tttdtttz ttzt ttttttdzdz zttttdzdz zzttdz zgddd8tttddtttdzzdttzzt tttttdtzztttdtzzztttdtzzztttdtzzzgddd9tttddtttzz d:ttzzt ttttdtzzttdtzzztdtzzgddd;tttddtttdzzdttzzt tttttzztttdzztttzztttdzzzgdddtttddtttdzztttzt ttttttzzttzztttzzttzzgddd?tttdtttdzttzt tttttdtzdz ztttdtzdz zztt dz ztttdtzdz ztttdtzdz zzttdz zzgddd@tttddtttdzzttzt tttttttd dzzztttttddzzzztttttd dzzzztttttddzzzzgdddAtttddtttzzdttzz t tttt|tzztt!|tzttt#|tzzttt#|tzzzzztt!|tzttt#|tzzttt#|tzzzzzgddBdCtttddtttzz ttzt tttt|tzztt|tzzztt| tzzztt!| tzttt#| tzzttt#| tzzzzzgddBtttddDtttdzz dEtttzzttzt tttt| tzztt| tzzztt|tzzztt|tzzztt|tzzzgddBtttddtttdzz dtttdzzd tttzzdFttzz t ttttdtzzt$tt|zzztt!|tzttt#|tzzttt#|tzzzzztt!|tzttt#|tzzttt#|tzzzzzgddBtttd dtttdzz dtttdzztttdzdtttd zz ttzt tttttzztt|zztttzzttt#|zz tt$tzzttt#|zzzttt!|zzzt&t(tzzttt#|zz t*tzt,zttt#|zzzttt!|zzzgddBt tdtttttztttzdt tttttdtddztzdz zzttdtddztzdz zzgdt |tttttztttzdt ttttttt|z zztttt|z zzgdt t.tttttztttzdt tttttttt.z zzttttt.z zzgdt t0tttttztttzdt ttttttt0 tzzztttt0tzzzgdttttt2ttzzt ttttt tt2 zzttttt2 zzzgdt dGttttzttzdt tttttdz ztttdz zzgdt dGttttzttt ttttt dz ztttdz zzgdt ttttdt tttzz dttzzdt ttttttzzttzgdt ttttdt tttzz dttzzdt tttttzztdtzzgdt ttttdt tttzzdttzzdt ttttttdzzttttdzzztt zgdttttdttzzt ttttdtzzttdtzzzgddHdIS)JNrTrr)realcZg|](}ttdzdtzzdz |)S)rrrrorps rTrzO_get_examples_ode_sol_nth_linear_constant_coeff_homogeneous.. s2GGG&A1q!44GGGrVrcZg|](}ttdzdtzz dz|)S)rrrrorps rTrzO_get_examples_ode_sol_nth_linear_constant_coeff_homogeneous.. s2GGG6!Q$1*q.!44GGGrVcvg|]6}ttdzdtdzzz dtzzdz|7S)rdrrrorps rTrzO_get_examples_ode_sol_nth_linear_constant_coeff_homogeneous.. s@YYY1vadSAXoQ&>&BAFFYYYrVcZg|](}ttdztdzz dz|)S)rrrrorps rTrzO_get_examples_ode_sol_nth_linear_constant_coeff_homogeneous.. s3MMMqvadQTkB&6::MMMrVcTg|]%}ttdztz dz|&S)rrrorps rTrzO_get_examples_ode_sol_nth_linear_constant_coeff_homogeneous.. s.III1vadQhlA66IIIrVr%nth_linear_constant_coeff_homogeneouslin_const_coeff_hom_01rrrlin_const_coeff_hom_02rlin_const_coeff_hom_03lin_const_coeff_hom_04rsrrrlin_const_coeff_hom_05rrlin_const_coeff_hom_06rrlin_const_coeff_hom_07rlin_const_coeff_hom_08lin_const_coeff_hom_09lin_const_coeff_hom_10lin_const_coeff_hom_11lin_const_coeff_hom_12rfilin_const_coeff_hom_13lin_const_coeff_hom_14lin_const_coeff_hom_15lin_const_coeff_hom_16relin_const_coeff_hom_17lin_const_coeff_hom_18lin_const_coeff_hom_19lin_const_coeff_hom_20rUlin_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_33rxlin_const_coeff_hom_34rr2ru)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_48r)r)rrr?r<rr rrrrrrr rrrrrC7C8C10C9r rr>)rrrrrrrr6r7r8r9r10r11r12r13r14r15r16r17r18r19r20r21r22r23r24r25rs rTrr s" sT"""AsAGGeAhhGGGBBBGGeAhhGGGBBCYYPUVWPXPXYYYCc3MME!HHMMMCc3IIaIIICc3 AA;aDDjddii1oo!A$$))A,,.1Q44bRTl*++,j ddii1oo!A$$))A,,.1Q4471Q44"r#a&&y.#a&&0112 jddii1oo!$1Q44CGGbQi/001j ddii1oo! !Q/!AaDDIIaLL.@1Q44bRTl*RAaC[899:!j,! !Q"QqTTYYq\\/1AaddF:1Q44C!HHr#aA.>*?*?'??@@A-j811a1QqTTYYq\\>1AaDD8!<<1Q44C2Q< 0111BsA2tAww{;K7L7L4LLMMN9jD1Q44AaddAq!1!11BtAaDD!}}4DDq1vM1Q44C!HHr#qb!d1gg+.>*?*?'??"SBQUVWQXQXLIYEZEZBZZ[[\EjPddii1oo! !Q/!AaDDIIaOO2CC !A$$))A,,1Q44bRTl*RAY6C!HHDEEF Qj^ddii1oo!A$$))Aq// 11AaDDIIaOOC !A$$))A,,1Q44 1Q44CGGbQi/2c477(1*oo3E3tTUwwWXy>>HY3Y[^_abc_c[d[d2ddeef _jlddii1oo1QqTT )1Q44Ca OO+bT!WWQY.??"SaQR^^BSSVXY\]abc]d]def]fYgYgVgghhimjxddii1oo!AaDDIIaLL 001QqTT691Q44C1tAqD1H~~#5 6777"SAQPQTTUXDVAW=X=X:XXYYZyjDddii1oo!AaDDIIaLL 002ad71Q44<?1Q44C1QKK"S1QZZ-7889EjPddii1oo1Q44bdR1W,r!Q$w6778Qj\ddii1oo!A$$))A,,.1Q4471Q44"r!t)SAYY.//0]jh! !Q!AaDDIIaOO"33addiillBQqTTI1Q44"r!t)S!WW,r#ac(({:;;<ijtddii1oo!A$$))Aq// 11Bqttyy||OCa!fL1Q44"q"r!t)},c!A#hh6778uj@ddii1oo!AaDDIIaLL 001a4!9<1Q44"r!t)S1XX-../AjjLddii1oo!A$$))Aq// 111Q44bdR1W,r#bd))|;<<=MjXddii1oo!A$$))Aq// 111Q44bdRT!WWHQJ%77"Sa^^:KKLLMYjdddii1oo!A$$))Aq// 11BqttyyA4FF 1Q4499Q<<QqTT'"1Q44"r!t)SAYY."r!t)S1XX1EEFFG ejr11a 2addii1oo#55!A$$))A,,F1Q44O1Q44CGGbaRTl2RAaC[@2c!HUVXYNNJZF[F[C[[\\]sj~ddii1oo!A$$))Aq// 11BqttG;1Q44"r!t)SAYY."r!t)S1XX1EEFFGjJddii1oo!A$$))A,,.1Q4471Q44"S1XX+3qs88 3SVV;<<=KjVddii1oo! ! ,qtt31Q44"S4771---3qay{3C3C0CCS1XXMNNOWjbddii1oo!A$$))Aq// 11AaddF:1Q44 3tAwwqy>>Bs47719~~--3tAwwqy>>0AABs4PQ77ST9~~DUUWWX cjpddii1oo!A$$))A,,.AaDD81Q44"S1XX+3qs88 3S1XX=>>?qj|ddii1oo!A$$))Aq// 11AaddF:1Q44"r!t)S477^^3rBqDy#aQi..6PPQQR}jHddii1oo!A$$&1Q44"S477^^+bQtAwwY.??QG"SQSTUQUYY,VWWXIjTddii1oo!A$$))Aq// 111Q44bdRAaC[02c!A#hh;>??@Uj`ddii1oo!A$$))Aq// 11AaDDIIaLL@1Q44rBqDy#a&&00BAIs1vv3EEFFGajlddii1oo! !Q/!A$$61Q44"Sa1---3tAwwqy{3C3C0CCS!AYYN c$q''!)A+ CQ ! $4$4!4 4c!A#hh >?@@A mjzddii1oo!A$$))Aq// 11AaDD81Q44C$Qx!|"4"4 45553qd1ggPQkARAR?R;S;S8SS S4a1 %%% & & &')+C$tAww{2C2C0C,D,D)DEFFG {jJddii1oo1Q4499Q<</!AaDD&81Q44 3r!t99 s2b66!8}}3r"vvax==(82c"R&&(mm;K(KLL bffQh--2c"R&&(mm+bRVVAX.>> ? @AAB! Kj\ddii!nnq11~-!41Q44 3r!t99 r#bd))|#bRTl2 bffQh--2c"R&&(mm+bRVVAX.>> ? @AAB! ]jjpddii!nns1Q4499Qq>>11D114EE!L1Q44CAJJCAJJ6CAJJFCPSTUPUJJVY[\_`cde`e\f\fYffggh ddii1oo!A$$))Aq// 11Aaddii1oo4EE1Q4499UV<<WZ\_`ab_c_cZcc1Q44 3qs88 S3ZZ-  bggai..Bs2c7719~~-3r#wwqy>>0AA B C bggai..Bs2c7719~~-3r#wwqy>>0AA B CDDE !ddii21Q4499Q??!22QqttyyA5FF A$$))Aq//addii1oo-.01!51Q44"r!t)S3ZZ/b2a4iQr#wwY5O3P 19c!BsGG)nn $4%&)!BsGG)nn355;=19c!BsGG)nn:T8U q52:s1RWW9~~ %9&'*1RWW9~~866778! a1Q4499Qq+++addiill:A>>1Q44bQA%6q%8%:!;!;;;bQQRTUEVWXEXYZEZA[A[>[[\\] QqTTYYq1%%%! ! 4a881Q44bQtAwwY//"S477^^2CCDDE addii!A&&&115q991Q44bQtBxxZ002c!DHH*oo3EEFFG QqTTYYq1%%%! ! 4a881Q44bd1ggXaZ002c$q''!)nn3DDEEF ddii1oo!A$$&1Q44C488 ,,,r#aaRj///AABBC 1Q4499Q??"QqTT)1--1Q44C!HHr#ac(({2334 1Q4499Q??"AaDD))1Q44C1II3qs88 3445 11a1T!A$$]]?2QqttV;Q??1Q44"r#a&&y.#a&&0112 11a1T!A$$]]?2QqttV;Q??1Q44"r!t)S1XX-../11a1T!A$$]]?2QqttV;Q??1Q44"S477^^+bQtAwwY.??aRHIIJddii1oo!A$$&1Q44C!HHr#ac(({2334M jjjn n nrVc dtttttz ttttz ztttz ttttz zttttz ztttzz t t ttt ttttttz ztz z gddttttzttz tttttz zz t t tt tt ttttz dz dz zt ttttz dzdz z gddtttttz zttttttz zz tttttz ztttzzt t ttttt tz ttttz zdz zttt tz ttttz zdz zgddttttttz z dttttz z zt tttttt tz dtzzzt tttttt tzzgddtttttz zttzttttzz t ttt dtt tz z tzgddd d S) NrTrrrFrrx)dep_div_indep_01dep_div_indep_02dep_div_indep_03dep_div_indep_04dep_div_indep_05r) r?r<rr rr rrrrr rWrVrTrr s?aDDdd1fS1a[[ AaddF3qttAv;;$6QqTT!V$DaddiiPQll#RR3q662AaDDQqTT!V$4Q$6 7 77889! ! nqtt#aAaDDF m33q663r77SQqTT!Vq%9%9!%;;c#add1f++PQ/>R>RST>TTUUVAaDDF madd3qttAv;;..3qttAv;;qttyy||1KK3q662QqTTE!G S1a[[ 8 ::S!A$$q\\#aPQddSTf++=UVW=WWXXYddiillQqTT!V#aAaDDF m31Q44T"s1vv+...25677AaDD!DcRSffDUDUBU9V9VWAaDDF madd"Qqttyy||^31Q44ac!ffo122334 9!!% % %rVc |dttdtttddttzztdzz ztttttdzdtzzdzz t ddz gdidS) Nlinear_coefficientslinear_coeff_01rrrsrurr)r?r<rr rrrWrVrTrr s*aDDddiilla!AaDD&j1q5111Q44QTAaCZ!^,x1~~=>>?   rVcXdtttttttttz zdtzz ttttzzt ttt t  tt t t dzzzgdddttzt tttz ztttttzzdtzt tttz ztttzz t tttt dt tttz zz gddtdzzttzttdzztttdzzdtdzzz tttzzt ttt dt ztdtd zzt d t zzdz ztz gddttttdztttzz ztttzttz t tttt dtt ttz dzzz gdd tttztttz tttzz t ttt ttdttdztdzz zz ttttz zgdttttzttz tttttz zz t ttdtztt tzzgdtdzttdzzdtzttztttzz t ttttt tzz t ttttt tzzgdttdzttttdztdzz ztttzz tttzzt ttt tt tdt zz z zt ttt tt tdt zz z zgddd d S)NrrrTrxrrrrrr)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_08r) r?r<rrr rrrrrr rWrVrTrr s1aDD ddaAaDDF mac)4!a==881Q44#b''(A2c2#'ll?";";;<<= &&!fS1Q44[[ 1Q44! ! #44qs3q1v;;qttyyQR||7SS3qtt99b1S1Q44[[=0112&&1fQqTTkAaDD!G#q1qy1QT6'91Q4499Q<<&GG1Q44QrTHR1WSBZZ-?$@$@$BBCCAEFFG &&4!a!AaDD&()))1Q4499Q<<7!A$$>3qtt99b1T1"QqTT'A+%6%6#66778&&!A$$h!add(AaDDIIaLL003q662DQqTT1WQT\)9$:$: ; ;;d1Q446llJKKL&& ! ! nqtt#aAaDDF m31Q441T"Q$ZZ(()&& dQqTT1Wnqs1Q44x! ! 441Q44$q"q&z***++R!d1b1f:6F6F-G-GH&& ddAg4!a!Q$///!AaDD&8!A$$))A,,FF1Q44"T1"a!B$h-000112addBtQBAbDM?R?RrGrrJr$r'rNr+rAr rrrrrrrrrr)rms rTrhrh: s :23()%%$ $ ( ( " "%%$$( (. .? ?$ $$ $//77 ,!,"1#1$=%=&,',($)$**+*,@-@.5/50C1C2.3.L6 rV)F)__doc__sympy.core.functionrrsympy.core.mulrsympy.core.numbersrrrr sympy.core.relationalr r sympy.core.singletonr sympy.core.symbolr r$sympy.functions.elementary.complexesrr&sympy.functions.elementary.exponentialrrr%sympy.functions.elementary.hyperbolicrrrr(sympy.functions.elementary.miscellaneousrr$sympy.functions.elementary.piecewiser(sympy.functions.elementary.trigonometricrrrrrr r!rr"r#sympy.functions.special.hyperr$sympy.integrals.integralsr%r&sympy.polys.rootoftoolsr' sympy.corer(r)sympy.functionsr*r+r,r-r.sympy.integrals.rischr/sympy.solvers.oder0r1sympy.solvers.ode.oder2r3sympy.solvers.ode.singler4r5r6r7r8sympy.solvers.ode.subscheckr9sympy.testing.pytestr:r;rr<r=rr>r?r@rrrrrrr%r&r(r'rrXrYrrrrrrrXrbror{r}r]rjrrrrrrrrrrrrrrrr rrrrrr"r%r(r,r0r5r7r:r<r?rBrErHrKrQrZrbrhrVr_rgrDr>rArGr9rJrNrr$r'r+rr rrrrrrrrrrrhrWrVrTrhs  <32222222333333333333********""""""........99999999GGGGGGGGGGKKKKKKKKKKKKAAAAAAAA::::::[[[[[[[[[[[[[[[[[[>>>>>>>>//////;;;;;;;;******''''''''HHHHHHHHHHHHHH77777722222222GGGGGGGG::::::::::::::433333-------- F3KK F3KK U3ZZ F3KK HSMM HSMM+277+;+;'BBBBC ''  1a  * " .''''JJJJ111888MMMM`<>4"*"*"*J@@@ G G G[[[ RRRfff<<<666>>>RRRCCC---OOOQQQ>>>:::]]]III Q Q QJJJ666AAA TTTeeeJJJQQQ @ @ @666777 KKK444AAA666777999 6 6 6 > > > E E E ; ; ;777t)))X###L2   HHHX   VVVt'''T000fTTTnQQQjDDDNpppf]]]@eeePAAAH@///dkkk\sssl$yyyx &&&P   555prV