K irdZddlmZmZddlmZddlmZddlm Z ddl m Z ddl m Z mZmZmZmZmZmZddlmZdd lmZdd lmZmZmZmZmZmZmZm Z dd l!m"Z"m#Z#dd l$m%Z%dd l&m'Z'm(Z(m)Z)ddl*m+Z+m,Z,m-Z-m.Z.ddl/m0Z0m1Z1m2Z2m3Z3ddl4m5Z5m6Z6m7Z7m8Z8ddl9m:Z:m;Z;mZ>m?Z?m@Z@mAZAddlBmCZCmDZDddlEmFZFmGZGmHZHddlImJZJddlKmLZLmMZMmNZNmOZOmPZPmQZQmRZRmSZSmTZTmUZUddlVmWZWddlXmYZYddlZm[Z[m\Z\ddl]m^Z^m_Z_ddl`maZaddlbmcZcmdZdddlemfZfddlgmhZhddlimjZjdd lkmlZlmmZmdd!lnmoZodd"lpmqZrmsZtmuZvmwZxdd#lymzZzdd$l{m|Z|dd%l}m~Z~dd&lmZdd'lmZdd(lmZmZmZmZdd)lmZdd*lmZdd+lmZdd,lmZdd-lmZdd.lmZmZmZmZdd/lmZdd0lmZdd1lmZmZdd2lmZdd3lmZmZdd4lmZdd5lmZdd6lmZdd7lZdd8l/mZdd9lXmZdd:lmZmZmZdd;lemZddlmZmZdd?lmZmZmZdd@lmZmZddAlmZmZddBlmZddClmZddDlmZddElmZddFlmZddGlmZddHlmZddIlmZmZmZmZmZmZmZddJlmZddKlmZmZmZddLlmZddMl!mZddNlmZmZddOlmZddPlmZddQlmZe%jZeZe)dR\ZZZe)dSdTU\ZZZZZZedVZedWZdXZdYZdZZd[Zd\Zd]Zd^Zd_Zd`ZdaZdbZdcZddZdeZdfZdgZdhZ diZ djZ dkZ dlZ dmZednZdoZdpZedqZdrZedsZdtZduZdvZdwZedxZedyZedzZed{Zed|Zed}Zed~ZedZ edZ!dZ"dZ#edZ$edZ%edZ&dZ'dZ(dZ)dZ*dZ+edZ,edZ-dZ.dZ/dZ0dZ1edZ2dZ3edZ4dZ5dZ6dZ7dZ8dedzzdedzzz dedzzz dedzzz dedzzz dz Z9dedzzdedzzz dedzzz dedzzz dedzzz dedzzzdzZ:dedzzdedzzz dedzzzdedzzzdezz dz Z;dZdZ?dZ@dZAdZBdZCedZDdZEdZFdZGdZHedZIdZJed„ZKedÄZLdĄZMdńZNdƄZOdDŽZPdȄZQdɄZRedʄZSed˄ZTd̄ZUd̈́ZVed΄ZWdτZXdЄZYedфZZd҄Z[dӄZ\edԄZ]edՄZ^edքZ_edׄZ`ed؄ZadلZbedګedۄZced܄Zdd݄ZedބZfed߄ZgdZhdZidZjdZkdZldZmdZndZoedZpdZqdZredZsedZtdZuedZvdZwdZxedZydZzdZ{dZ|dZ}dZ~dZdZdZdZdZdZedZdZedZedZedZedZdZedZdZdZedZedZd Zed Zd Zed Zed ZdZedZedZedZdZdZdZdZdZdZdZdZdZedZdZdZdZedZed Zed!Zd"Zd#Zd$Zd%Zd&Zd'Zed(Zed)Zed*Zed+Zeed,Zed-Zeed.Zd/Zd0Zd1Zd2Zd3Zed4Zd5Zd6Zed7Zd8Zd9Zed:Zd;Zed<ZĐd=ZŐd>ZƐd?Zǐd@ZedAZɐdBZʐdCZːdDZ̐dEZ͐dFZedGZedHZАdIZѐdJZҐdKZӐdLZԐdMZՐdNZ֐dOZאdPZؐdQZedRZڐdSZېdTZܐdUZݐdVZސdWZߐdXZdYZedZZed[Zd\Zed]Zd^Zd_Zed`ZedaZedbZdcZeddZedeZdfZdgZdhZdiZdjZdkZedlZedmZednZedoZedpZdqZdrZedsZdtZeduZedvZedwZedxZedyZedzZd{Zd|Zd}Zeed~ZedZedZedZ dZ dZ dZ dZ dZdZedZedZedZedZedZdZdZdZdZdZdZdZdZedZedZedZdZ dZ!dZ"dZ#dZ$dZ%edZ&dZ'dZ(edZ)dZ*dZ+dZ,edZ-edZ.dZ/edZ0edZ1edZ2edZ3dZ4dZ5dZ6dZ7edZ8edZ9dZ:edZ;dZ<dZ=dZ>edZ?eedZ@dZAedZBedZCedZDedZEeedZFeedZGeedZHeedZId„ZJedÄZKeedĄZLedńZMedƄZNdDŽZOedȄZPdɄZQedʄZRd˄ZSd̄ZTd̈́ZUed΄ZVedτZWdЄZXdфZYed҄ZZdӄZ[eedԄZ\eedՄZ]dքZ^dׄZ_d؄Z`dلZadڄZbdۄZced܄Zdd݄ZedބZfd߄ZgdZhdZidZjedZkdZledZmeddZndZoedZpedZqedZredZsdZtedZudZvdZwdZxdZydZzedZ{edZ|dZ}dZ~ededZedZedZedZdZdZdZedZedZedZy7(z Tests from Michael Wester's 1999 paper "Review of CAS mathematical capabilities". http://www.math.unm.edu/~wester/cas/book/Wester.pdf See also http://math.unm.edu/~wester/cas_review.html for detailed output of each tested system. )Qask)refine)product) EulerGamma)N) DerivativeFunctionLambdaSubsdiffexpand expand_func)Muligcd)AlgebraicNumberEIRationalnanoopizoo)EqLt)S)DummySymbolsymbols)rfbinomial factorial factorial2) bernoulli fibonaccitotient partition) conjugateimresign)LambertWexplog)asinhcoshsinhtanh)ceilingfloor)MaxMinsqrt) Piecewise) acosacotasinatancoscotcscsecsintan)besselj) DiracDelta) elliptic_e elliptic_f)gamma polygamma)hyper)assoc_legendre chebyshevt)polylog)idiff)And)hessian wronskian)MatMul)continued_fraction_convergentscontinued_fraction_iteratorcontinued_fraction_periodiccontinued_fraction_reduce) factorint) primerange)ZZ) legendre_poly)apart)Polyfactorgcd resultant)limit)O)residue)series)ImageSet) FiniteSet IntersectionIntervalUnion)combsimp) hyperexpand) powdenestpowsimp)radsimp) logcombinesimplify) sqrtdenest)trigsimp)solveNstirling) Heaviside)CiSierf)zeta)XFAILslowSKIPtooslowraises) partitions)mpimpc)Matrix GramSchmidteye) BlockMatrixblock_collapse) MatrixSymbol ZeroMatrix) Commutator)PolyRing) FracField) solve_lin_sys)Sum)Product) integrate)laplace_transforminverse_laplace_transformLaplaceTransformfourier_transformmellin_transformlaplace_correspondencelaplace_initial_conds)rsolve)solveset solveset_reallinsolve)dsolve)Equality)islice takewhile)fps)fourier_series)minimumx y zz i j k l m nTintegerfgctttttttttttzttt tztttttt k(sJyNreijklmg/mnt/ssd/data/python-lab/Trading/venv/lib/python3.12/site-packages/sympy/utilities/tests/test_wester.pytest_B1rvsR aAq!Q ')Aq!*< < aA  #,Q1a#; << Ad2hJ648CD ED 7T!W tAw &D 4<>FFtFT UYZ ZZ Zrc ttddtdzztddzdtdzk(sJy)NZ"rrr)rmror8rrrrtest_C19rs9 8R"T!W*_1a89 :a$q'k II Ircddtdzz}t|tddzdztdz|tddzz }t|tdk(sJy)NNrrrr)r8rrro)insiders rtest_C20rs\BtAwJF FAaGOa/47:VQq!W_L MD D>_R0 00 0rc ttddtdzztddztdtdzk(sJy)Nrrrrr)rorr8rrrrtest_C21rsE OR"T!W*_1a$@A BDG $ %% %rc ptddtdzz tddtdzz zddtdzz tddtdzz zzdzdtdzz d tdzd z z }tddz ttdd z dz z }||k(sJy) Nrrrrrr 0Hr)ror8r/rs rtest_C22rs a!DG)mSQtAwY%771qay=#b T!W*OK;QZ(,.tAwJO= >D 719s47Q;') )D 4<0QTWXY[\T]0] ]] ]rctd)Nz:discretize a PDE: diff(f(x,t),t) == diff(diff(f(x,t),x),x)r"rrrtest_D13rJLs Z [[rc^ttdtdtzzdtzzk(sJy)Nrrr)r!rrrrtest_F1rLVs' a8q!a%y!a%( (( (rcvtttdttdz ztdz zdz k(sJy)Nrrrr)rr"rrrrtest_F2rNZs1 x1~ &!QU)QU*;A*= == =rctdtzttztdtzdz ztdtzk(sJyNrr)rirr#r$rrrtest_F3rQ^s< AqD9Q<'*QqS1W*== >)AaC. PP Prc tdtzttztdtzdz tdtfztdtzk(sJyrP)rirr#rrrrrtest_F4rScsG AqD9Q<''!A#'Aq!9*EE F)TUVWTW. XX Xrctttddzttz t tz t dtzddtzzz t tdzz k(sJyNrr)rHrrr8rr#rrrtest_F5rVhsY Qq!W R (9Q< 79QqS>!aPQc(;RS\]^S_abSb;b bb brctdDcgc]}|j}}ddidddddidddddig}||k(sJycc}w)Nrr)rrr)rr)rcopy)ppartTest partDesireds rtest_F6r\msV",Q-0Q0H0q6qQztest_G15..sACC2Irr2) rr8r*limit_denominatorrrfrcf_ccf_irrrtest_G15rusj DGMMO $ 6 6r :aBi GG G  -tDaM/BC DR H "b   rcVttttdgdk(sJy)Nr) rrrri$rrrrr)rfrrtrrrrtest_G16rws! tBx# $(I II Irc0tddddgdgk(sJy)Nrrrr)rrrr)cf_prrrtest_G17rzs 1b>a. .. .rctddddggk(sJtdggjtjt ddz zk(sJy)Nrrr)rycf_rrrHalfr8rrrtest_G18r~sL 1a=aSE !! ! !;   166DGAI#5 55 5rctddd}ttd|z dz td|z dzz }tt |ddd|zd|zd |zd |zgk(sJy) NsTrpositiverrrrrrr)r rtr.rfr)rits rtest_G19rskTD1A s1Q3x!|c!A#hl+ ,B r1 1ac1Q31bd"; ;; ;rcntddd}td|zgg|t|dzdzzk(sJy)NrTrrr)r r|r8rs rtest_G20rs;TD1A !A#=AQTAX. .. .rcZtddd}t|d|dzdzd|zggk(sJy)NrTrrr)r ryrs r test_G20brs7TD1A 1adQh QqSE7 ** *rctddtzztdtdzzk(sJtddtzztdtdzzk(sJyrP)rorrkrrrtest_H1rsR AadF xAE 3 33 3 Qq!tV QU 4 44 4rcLtddtzzdtdzzk(sJy)Nrr)rlrrrrtest_H2rs% 1q!t8 AE ** *rc4dttdzzzdk(sJyNr2rrrrrtest_H3rs !QU)  !! !rctdtzdz }t|tusJ|jddk(sJ|jddtzdz k(sJy)Nrrrrrrr)r]rtyperargs)exprs rtest_H4rsZ !A#( D :   99Q<1   99Q<1Q37 "" "r@rr~r.rr<Qrrrr'S46rFrr[Vc@ttttdk(sJyNr)r^p1p2rrrrtest_H5rs r2q>Q  rc~ttttztttztk(sJyr)r^rrrorrrrtest_H6rs' vb1f~vb1f~ .! 33 3rcdtztdzztdzzdtdzztdzztdzzz dtdzztd zztd zzzd td zzz dz}d tdzztdzztdzzdtdzztdzztdzzzdtd zztdzztdzzzdtdzztzz}t||tttdk(sJy)NrrrrrrrrrrrrrrrrrrrPrrryzr^rrs rtest_H7rs AaeAqD2ae8AqD=A- -!R%1 QT0A AAaeG Ka OB AqDAae bAgadl1a4/ /"QT'!R%-12D Dr!R%xPQz QB r2q!Q 1 $$ $rcdtztdzztdzzdtdzztdzztdzzz dtdzztd zztd zzzd td zzz dz}d tdzztdzztdzzdtdzztdzztdzzzdtd zztdzztdzzzdtdzztzz}dtdzztdzztdzzdtd zztdzztdzzz dtdzztdzztdzzz}t||z||zttt|k(sJy)Nrrrrrrrrrrrrrrrrrrrrrrrrr)rrros rtest_H8rsR AaeAqD2ae8AqD=A- -!R%1 QT0A AAaeG Ka OB AqDAae bAgadl1a4/ /"QT'!R%-12D Dr!R%xPQz QB 1b5A aebAgadl1b5002ae8AqD=A3EEA rAvrAvq!Q '1 ,, ,rctdd}d|tdzzz|tdzzz }d|tdzzzd|tzzz}t|||tzk(sJy)NrF)zerorrrr)rrr^)rrrs rtest_H9rsbsA 1q1u:AE "B 1q1u:!Q$ B r2;!Q$  rcdtdzzdtdzzztdzztz dz }tdzdtdzzz tzdz}t||tdk(sJy)Nrrrrr)rr_rs rtest_H10rsd 1a4!AqD&1a4 ! #a 'B A!Q$ Q B RQ 1 $$ $rc\tttzttztdk(sJyNr)r_rrorrrrrtest_H11rs! R!VR!VQ '1 ,, ,rctdzdz }tdzdtzzdz}t||z tdz tdzz k(sJy)Nrr)rro)r.dens rtest_H12rsG Q$(C Q$1*q.C CG QQ // /rctttdz ttdz dzz ttdz dz k(sJyrU)ror.rrrrtest_H13rs: SVaZC!HqL1 2c!A#hl BB Brctdzdz}t|}|ddtzzdtdzzzdtdzzzdtdzzzd td zzzd td zzzd tdzzzdtdzzzdtdzzzdtdzzzdtdzzzdtdzzzd tdzzzd tdzzzd tdzzzdtdzzzdtdzzzdtdzzzdtdzzztdzzk(sJt|t}|ddtzzdtdzzzd tdzzzd tdzzzd!td zzzd"td zzzd#tdzzzd$tdzzzd%tdzzzd%tdzzzd$tdzzzd#tdzzzd"tdzzzd!tdzzzd tdzzzd tdzzzdtdzzzdtdzzzdtdzzzk(sJt|ddtzdzzk(sJy)&Nrrritriri<rihri.riririrrrrrrrrri|i\ iKipiGi`ii1)rrr r])rYepdeps rtest_H14rs Q2 A B !bd(SAX%QT 1DAI=ad J 1*QTz"$*1a4K028A+>@Fq"u M B,2&(-ae 46;ArEkBDI!R%KP q"u*ArEz"%(2I.131b59<=b5A BB B r1+C 2A:QT )E!Q$J6q!tC A+q!t $&-adl35D 4<rrCrr])rs rtest_H26rPsyA3q619$qQ{2R78A !E "Aw3q619'T;U7U VF T>>rctdzdtzzdztdzdtdzzzdtzzdzz }dtdzdzz dtdzz z dtdzz z}t||k(sJy)Nrrrrr)rr[rs rtest_H31rosu A!a!Q$1a4-!A#-12A QUQJa!e$qAE{2A 8q==rctd)Nzp[A*B*C - (A*B*C)**(-1)]*A*C*B (product of a non-commuting product and its inverse)r"rrrtest_H32rus K LLrctdd\}}}t|t||t|t||zt|t||zjjdk(sJy)NzA, B, CF) commutativer)r rdoitr)ABCs rtest_H33r {sriU3GAq! q*Q* + Q 1a( ) * Q 1a( ) *+/46&&(a @@ @rcxtttddztddtdz z k(sJy)Nrrrrr)rCrrr8rrrtest_I1r s3 r!Ar({ Q47] 33 33 3rcVtdtdzdz td k(sJy)Nrrrr)r8r>rrrtest_I2rs' SVQ CF7 ** *rctttztdtzdz tzdz zdtzdz k(sJy)Nrrrr2)r>rrrBrrrtest_I3rs: qt9sAaC!GR<>* *rAgk 99 9rc ttttttzzt tdz tttzzzt j tdtzdz k(sJy)Nrr2r)rr>rrrBrrrrrtest_I4rsT #bQrTl#c"Q$s1R4y.&99199Q< HRRSGVWK WW Wrcttdzdz tdzdz ztdzdz ztdz z tzdk(sJy)Nrrrrrr)rBrrrrrtest_I5rsD 1QAa!Q$q&(1R4/25 6! ;; ;rctd)NzHassuming -3*pirrBrrrtest_I7rs9 qs8CF?c!fai!CFAI+5 55 5rc~tdtzttz dtdtzzdz k(sJy)Nrrr)r>rrrrtest_I8rs0 qs8CF?aAaCj1n ,, ,rctdtzttz ttdzdttdzzz k(sJyrrrrrtest_I9rs; qs8CF?c!fai!CFAI+5 55 5rctttdzdzttdzz t tdzttdzzdz z t usJyNrrr )rqrCrr>rBrrrrtest_I10r"sN SVQY]SVRZ/CFAIA 4IA4MN OSV VV Vrhangsctttdzdzttdzz t tdzttdzzdz z tddk7sJy)Nrrr r)r`rCrr>rBrrrtest_I11r%sU #a&!)a-#a&"*,QSVQY1F1JKQPQ RVW WW Wrctttdzdzttdzz t tdzttdzzdz z t}|t usJyr!)r rCrr>rBr)ress rtest_I12r(sT A A A *s1vqy3q619/Dq/HI1 MC #::rc8tdtddk(sJy)Nrii)r%rrrrtest_J1r*s R=AeSM )) )rcttttdzttttdzt ttdzz tz k(sJyNr)r rFrrrGrrrtest_J2r-sF  1ad#Q 'Jq!Q$,?*QPQSTPTBU,UWX+X XX Xrctd)NzDJacobi elliptic functions: diff(dn(u,k), u) == -k**2*sn(u,k)*cn(u,k)r"rrrtest_J3r/s d eercXttdddttzk(sJy)Nr2rr )rHrr8rrrrtest_J4r1s" 2q?bbk )) )rctdtddtd tdtzdz z t z ttdz k(sJy)Nrrrr)rIrr/r8rrrrrtest_J5r3sG Q!Q SVGd1gbjl$:Z$G#dSTg,$V VV VrcftjddjtddsJy)Nry??z0.04157988694396212z0.24739764151330632)mpmathrDr/rrrrtest_J6r6s* >>!V $ ' ',ACX(Y ZZ Zrcntttddtdz dtdzz k(sJy)Nr rr)rorDrrrrrtest_J7r8s. GAbGRT* +r2q5z 99 9rcttddt}tttz t tz t t tzdz z }tt||z dk(sJyNrrr) rDrrrBr>r8rror)rYros rtest_J8r;sW!AA QCF DAaL(A KNA% &! ++ +rcrtdtjttdt k(sJy)Nrr)rDrr rrrtest_J9r=s) 1a=  a gamO 33 3rctdd\}}t||dd|zttzt ||z dz dzz t | |z dzdz z k(sJy)Nzmu, nuTrrrr)r rKr8rrH)munus rtest_J10rAsk Xt ,FB "b! $2d2hub2gq[1_7M(MeVXUX[]U]`aUacdTdNe(e ee erc ttddtttdd t dtdzz zdtdzzdz zk(sJy)Nrrrr)rorKrrr8rrrtest_J11rCsT N1a+ ,!Aq'$q1a4x.9PRSTUWXTXRX[\R\9]0^ ^^ ^rcttdtdtztdtzz tdtzdk(sJy)Niriir)rorLrrrrtest_J12rEs? JtQ'!A#jq.A*AAJtUVDWW X\] ]] ]rcHtddd}t|dd|zk(sJy)NrTF)rnegativer2)r rLrs rtest_J13rHs)TE2A a a '' 'rcttjtjgtddgtdz}t |t ttz k(sJyr)rJrr}rrrjr<)rYs rtest_J14rJsC qvvqvv1a 1a40A q>T!WQY && &rctd)NzcF((n+2)/2,-(n-2)/2,R(3,2),sin(z)**2) == sin(n*z)/(n*sin(z)*cos(z)); F(.) is hypergeometric functionr"rrrtest_J15rLs D EErctd)Nz&diff(zeta(x), x) @ x=0 == -log(2*pi)/2r"rrrtest_J16rNr;rc ntttdzdz ttdz dz zt tt ttdz tzz tddfdtt ddzttt tttdzk(sJy)Nrrrrrr) rrrrErr rr r rrrtest_J17rP s QAqy\*a!eQY"77!A$tJqSTuDUWX?Y:YY\]_`bc[d eijklmnoprsmtkuiux|~HIJKLIMOP~QSTVWyXjX XX Xrctd)Nz define an antisymmetric functionr"rrrtest_J18rRs @ AArctdd\}}t|t|zzt| t|zk(sJt|t|zzt|t|zk(sJy)Nzz1, z2Tr )r r+rr*)z1z2s rtest_K1rVsa Xt ,FB b1R4i=RVGbf, ,, , b1R4i=BrFRVO ++ +rc tdtdz ttdtdzdz zzdk(sJy)Nrrrrr)absr8rrrrtest_K2rYs8 q47{QtAd1gIN333 4 99 9rctdd\}}ttd|t|z zt|zzz dt |dzt|z |zdzzz k(sJy)Nza, bTrealrr)r rorXrr8rbs rtest_K3r_ s^ 6 %DAq C1qs7QqS=)* +qadacAg\6I1J/J JJ Jrc tddtzzjdtdttt ddzzk(sJy)NrrTr r)r/rrr=rrrrtest_K4ra&sA q1Q3w<  t  ,A4!Q=0H HH Hrc4tdd\}}t|t|zzjdt d|zt d|zt d|zzz ttd|zzt d|zt d|zzz zk(sJy)Nzx, yTr[r r)r rCrrrBr>r1r2rrs rtest_K5rd*s 6 %DAq q1Q3w<  t  ,QqS3qs8 QqS <2tAaCy[#ac(T!A#Y"6728 99 9rctttzttdzzttttzz tttzttz k(sJtttzttdzzttttzz ttk7sJyr,)r8rrrXrrrrtest_K6rf0s !CFAI QA /4!9T!W3D DD D !CFAI QA /47 :: :rctddd}tt|zttdzzttttzz }t |}|t|k(sJy)NrTF)r\rGr)r r8rrXrro)rrsexprs rtest_K7ri5sV$/A !CFAI QA /D TNE DG  rctdd}ttd|z dt|z z dk7sJtddd}ttd|z dt|z z dk(sJy)NrTr rrF)r rGr ror8rs rtest_K8rm<slT"A D1I$q' ) *a // /TE2A D1I$q' ) *a // /rcttdd}ttd|z dt|z z dk(sJy)NrTrrrrkrls rtest_K9rpC6d#A D1I$q' ) *a // /rcttdd}ttd|z dt|z zdk(sJy)NrT)rGrrrkrls rtest_K10rsHrqrcDtddtddzz dk(sJy)Nii;rrrr8rrrrtest_L1rvQs$ 91a( (A -- -rcDtddtddzz dk(sJy)Nrel^|E!orrrrurrrtest_L2rxUs$ <9qAw. .! 33 3rc tdtddzdtddzzdzddtddzdtddzzzz dz dk(sJy)Nrrrrrr)rorrrrtest_L3rzYs] Q!QZ!Qq!W*,q01a1aj1a1g:6M3NNQRR SWX XX Xrctttdzttttdzzzttz dk(sJyr:)rqr>rrBrrrtest_L4r|]s< CFAIAs1vqy 003q69 :a ?? ?rctttddtztdz zt ttz dk(sJy)Nrrrr)r/rCrrrr0rrrtest_L5r~as: s1Q719r!t#$ %c!f 5 :: :rctttdz tdz zt ttz j tj tdidk(sJy)Nrrr)r/rCrrr0r subsrrrtest_L6rfsL AaC"Q$J 5Q= 0 6 6q 9 > >1v F! KK Krc ttdttzdztdtzdttzzdzz dk(sJy)Nrrrr)ror/r8rrrrtest_L7rjsE C47QacAd1gIo.A)BCD E JJ Jrc 4tdtzdttzzdzttdttzdzz zdttzdzddttzdzz zzdttzz dz dk(sJy)Nrrrr)rorr8rrrtest_L8ros QqS1T!W9_q(DGQtAwY],CD DG)a-1aQi!m, -.01$q' :<=> ?BC DD Drctdd}tdd|z zt|zt|zt |t zdz zt dztd|z zz dk(sJy)NrTr rrr)r rorHryr>rrls rtest_L9ruscT"A AAJuQx'Q/AbDF ;b!eDQK>OO PTU UU Urcvttddz tddzttdz dzdk(sJyrP)rrrrrtest_M1r}s5 Aq>! hq!n ,1q!0D DD Drctdtdzzdtdzzz dtzzdz t}td|DsJy)Nrrr!rc3TK|] }|jdj"yw)Tr N)ris_real).0rs r ztest_M2..s!;!qxxx%--;s&()rrall)sols rtest_M2rsG 1QT6Bq!tG#bd*R/ 3C ;s; ;; ;rc ttdzdtdzzz dtdzzz dtdzzzdtzz dz ttd tdzd tdz tdd z ttdzd z ztdd z ttdzd z z td d z ttdzd z ztd d z ttdzd z z k(sJy) Nrrrrr$rgr(?gr(?grrrrBrs rtest_M6rsr x1q!$ %:?1+FQQrT!Aq'\ Qs1R4!Q<00 0F GG GFsA Bctttdzdtdzzz dtdzzzdtdzzz dtdzzzd td zzz d td zzzd tzz dztdtd tzttddtd zz d zzd z z dtd tdttd dtd zzzzzd z z dtd tzttddtd zzd zzd z z dtd tdttd dtd zzzz zd z z dtd tzttddtd zzd zzd z zdtd tdttd dtd zzzz zd z zdtd tdttd dtd zzzzzd z zdtd tzttddtd zz d zzd z zhk(sJy)Nrrrr\rrrrrrr)rrrrr8rrrrtest_M7rsD uQTAadF]R1W,r!Q$w6QTAC1HL AqDq5 " # DGAIdDaQi00145 5a 77 DGDaQ47] 3334 4Q 66 DGAId4QtAwY/!34 4Q 66 DGDaQ47] 3334 4Q 66 DGAId4QtAwY/!34 4Q 66 DGDaQ47] 3334 4Q 66 DGDaQ47] 3334 4Q 66 DGAIdDaQi00145 5a 77 '   rc Btd}tdd}ttd|zdt|zzdz|z |tj t td|zdt|zz dz td|zdt|zzdz k(sJy)NrrTr rr) rr rr.rRealsrer/r8)rrs rtest_M8rss AT"A C!HqQx'!+a/AGG <#a!eaQi'(*CA$q' 0A,B1,DE FF Frctd)Nz6solveset(exp(2-x**2)-exp(-x),x) has complex solutions.r"rrrtest_M9rs V WWrclttttz ttd gk(sJy)Nr2)rrr.rr-rrrtest_M10rs' Q!Q Xb\M? 22 2rcftttztz ttddk(sJyr)rrrerrrtest_M11rs& AqD1Ha Ib!$4 44 4rcttdzttdzdzdzztdtzdzztdtdz tdz t t dtdzzt t dtdzztt t dtdzzz tt t dtdzzzgk(sJy)Nrrrr2r)rrrrBr>rrr/r8rrrtest_M12rs !a%#a&!)a-!++C!HaK7 ; BqD"Q$ 3s1tAw; aAQK(8&8 Qs1tAw; aAQK(8&8!8@  rc td}tttt tz tt t ||tzttddzz tjk(sJy)Nrrr) rrrBrr>rdr rrrIntegersrs rtest_M13rsU c A Q#a&! ,1R4"QqRSW*CT9UWXWaWa0b bb brc td}tttdz tt t ||t zt dz ztjk(sJy)Nrrr) rrrCrrdr rrrrs rtest_M14rsF c A Q!Q '8F1adRTk4JAJJ+W WW Wrctd}ttttj z t fdttt|d|ztztdz ztjtt|d|ztzttddzztjttt|d|ztzttddzztjtt|d|ztztdz ztjfDsJy)Nrc3@K|]}j|ywr)dummy_eq)rrgots rrztest_M15..s:1s||A:srrr) rrrBrrr}anyrhrdr rrr)rrs @rtest_M15rs c A 3q6AFF? #C : hva1R"Q$/<1Q3r6BqAwJ./< > hva1R"Qq!W*!45qzzB1Q3r6BqD=)1::6 8)9: :: :rctd}tttt tz tt t ||tztjk(sJy)Nr) rrrBrrCrdr rrrrs rtest_M16rsA c A CFSVOQ '8F1adOQZZ+P PP Prcztttttz tt dk(sJyr)rr<rr=rerrrtest_M17rs( a47*A .)A, >> >rctttttz tt t t ddz dz k(sJy)Nrrr)rr:rr=rer8rrrtest_M18rs: a47*A .)D$q'A+qr.)rr^s rtest_M27r s$A$A T$q!Aq'zA~./!34q8! 3s3r7Q;'((1Q72 3 44 4rctdtzttdz dz zdtdzzz ttjtgdk(sJy)Nrrrr)assume)gOqg0eg?)rrr.rr\rrrtest_M28r*sD 1sAE19~-!Q$6!&&) LPp pp prctd}tt|dz dz tjt ddk(sJy)Nrrrdomainr2r)r rrXrrrerps rtest_M29r/s7 A CAJN177 3yQ7G GG Grcttdtzdzttdz z ttddk(sJy)Nrrr2)rrXrrerrrtest_M30r4s9 QqS1WAE 2A 6)B:K KK Krctdttz tt dz tdz z tt t ddt ddk(sJy)Nrrrr)rrXrr6rerrrrtest_M31r;sO SVc1"q&!a%&88! < !BPQ(TUVWYZT[@\ \\ \rcttdtdzz ttt tdzdz z ttddk(sJy)Nrrrr2rr6rrerrrtest_M32rBsF QAXq)3rAqD!8+< BiPQTUVWTWPWFX XX Xrctdd\}}td|zd|zz t|zz dtzzj||t dk(sJy)Nx yTr[rr)rr)r rr as_real_imagrercs rtest_M35rXsT 5t $DAq QqS1Q3Y1_qs*88:Aq AYvEV VV Vrctttdzttzdz tttddk(sJy)Nrr r)rrrrerrrtest_M36r]s9 AaD!GadNQ&! -2q1A AA Arctttztzdz dtztzdtzzdz tdtzztzdz gtttt t dzdtfk(sJy)Nrrrrr)rrrrrerrrtest_M37rdsm QUQY]AaC!GacMB$6AaC! b8HJAqRS TA261a.! "" "rc1|td\}}}t|||gtj}t d|}|j \1}}}}} } } } } }}}}}}}}}}}}}}}}}}} }!}"}#}$}%}&}'}(})}*}+},}-}.}/}0}1}2}3}4}5g| | z|z || z|z z| |z|z ||z|z z| |z|z ||z|z z|z| || z|z z || z|z z| |z|z ||z|z z| |z|z ||z|z z| |z|z ||z|z z|z | |z|z ||z|z z| |z|z ||z|z z|z| |z|z ||z|z z||z|z z ||z|z z|| z|z z|| z|z z ||0z|z ||0z|z z | |1z|z ||1z|z z| |z|z ||z|z z| |0z|z ||0z|z z||2z|z ||2z|z z |dz|3z|dzz d|z|z|3z|dzz z |dz|3z|dzz z|| | |z|z ||z|z z|| z|z z || z|z z| |z|z ||z|z z|| z|z z|| z|z z ||1z|z ||1z|z z | |2z|z ||2z|z z|4 ||4z|z z||4z|z z|dz|4z||zz z |5 ||5z|z z||5z|z z|dz|5z||zz z ||z|z ||z|z z ||z|z ||z|z z ||z|z ||z|z z | z| ||z|z z||z|z z || z|z || z|z z | |||z|z ||z|z z || z|z z|| z|z z ||z|z ||z|z z |||z|z z||z|z z | ||z|z z||z|z z ||z|z ||z|z z | |z|z ||z|z z||z|z z||z|z z || z|z z|| z|z z | |0z|z ||0z|z z||1z|z ||1z|z z ||z|z ||z|z z ||z|z z||z|z z ||0z|z ||0z|z z | |2z|z ||2z|z z|3 ||3z|z z||3z|z z|dz|3z||zz z ||z|z ||z|z z | |1z|z ||1z|z z||2z|z ||2z|z z |dz|4z|dzz d|z|z|4z|dzz z |dz|4z|dzz z|5 ||5z|z z||5z|z z|dz|5z||zz z || | |z|z ||z|z z| ||z|z z ||z|z z| |z|z ||z|z z|||z|z z ||z|z z||z|z ||z|z z || z|z z || z|z z| | z|z || z|z z| | z|z || z|z z| |z|z ||z|z z|z| |z|z ||z|z z|z | | z|z || z|z z| |z|z ||z|z z||z|z z ||z|z z| |z|z ||z|z z||z|z z||z|z z || z|z z || z|z z||0z|z ||0z|z z | |1z|z ||1z|z z| |0z|z ||0z|z z||2z|z ||2z|z z |3 ||3z|z z||3z|z z|dz|3z||zz z | |z|z ||z|z z||1z|z ||1z|z z | |2z|z ||2z|z z|4 ||4z|z z||4z|z z|dz|4z||zz z |dz|5z|dzz d|z|z|5z|dzz z |dz|5z|dzz z| || ||z | | | |! ||$z | |z|#z | |"z |' |0|1 |(||z |+z| |*z||)z||z|z ||z|z z |&z |.zd|0z|1|2||3z|z ||3z|z z |-|0|2 | |3z|z ||3z|z z||z| |z | |z|z ||z|z z|%z |1|2 | |4z|z ||4z|z z|| z|z || z|z z |,z|1 |2||4z|z ||4z|z z ||5z|z ||5z|z z | |5z|z ||5z|z z| | | |||z ||||| z ||"z||z|#z | |0 |1|" | z ||$z|||%z|z z||%z|z z |0|2 | |3z|z ||3z|z z|( |#|+z | z |$ |*z ||)z |||&z|z z ||&z|z z|.z |0d|1z|2||4z|z ||4z|z z ||'z|z ||'z|z z |-z |0 |2||3z|z ||3z|z z | |5z|z ||5z|z z|, |1|2 | |4z|z ||4z|z z||5z|z ||5z|z z |||| | | |z| |)| z| |*z | | |)z | | z| |z |+z| ||,z|z z||,z|z z |0|1 | |0 |2||3z|z ||3z|z z |!|$|*z|#|+z | z| |"z|'||-z|z z ||-z|z z|0 |1| |&z||.z|z z||.z|z z |0|1d|2z| |3z|z ||3z|z z| |4z|z ||4z|z z||5z|z ||5z|z z |%|1 |2||4z|z ||4z|z z | |5z|z ||5z|z z}6i|5d|4d|3d|2d|1d|0d|-d|,d|*d|)d|(d|'d|%d|$d|"d|!d| di|d|d|d|d|d|d|d|d|d|d|d|d|d|d|d|d| d| d| d| d| d|d|d|d|d|&||z |.z|#|+|||z |.z||+i }7t |6||7k(sJy)Nza, b, czk1:50rr r)r rrY to_domainrgensr)8rr^rrringk1k2k3k4k5k6k7k8k9k10k11k12k13k14k15k16k17k18k19k20k21k22k23k24k25k26k27k28k29k30k31k32k33k34k35k36k37k38k39k40k41k42k43k44k45k46k47k48k49systemrs8 rtest_M38r&isi GAq! 1ay" % / / 1F GV $D 4899 1RRRRRc3S#sCc3S#sCc3S#sCc3S#sC( 2a!B$q&( "S&(QsU1W,( /0bfQh3q.@2.E( HJsQrTRSV|VWXZVZ[\V\G\( 3q1S57( BsF1HquQw.( 123q1S570BR0G( JKCPQTUVYTYZ[T[I[(  3q1S57R( #$CAcE!G!3aeAg!=#a!G!B$q&!PSTUWSWXYSY!Y(  #a!C%' ( 2c6!8aeAg- ( 12r#vax!C%'/A ( EF2c6!8aPSeTUgCU ( #a!C%' ( a48AqD=1Q3q59QT>9AqDHQTMI ( LN ( QSs ( 3q1S57QrT!V#ad1f, ( 01bfQh3q.@1R46.IAbDQRF.R (  #a!C%'( 2c6!8aeAg-( 14taeAg~#a/G!Q$s(TUVWTW./X(  quQw3q 1a48QqS>1( 45R46AbDF?( EFbDFQrTRSVO(  "Q2a"(  #dQrT!Vmad1f4( 78d1fqtAvo( IKs( MP(  #a!C%'AbDF"QrT!V+( ./sU1WquQw->( ADaeAg PQRUPUVWPW@W(  quQw3q ( #$C%'AcE!G"3( 78RVAX#a5G!C%PQ'5QTUVYTYZ[T[5[^_`b^bcd^d5dghikgklmgm5m(  3q1S57( cE!GaeAg-( 01uQw3q/@1S57/JQsUSTW/T(  #a!C%'( 2c6!8aeAg-( 14taeAg~#a/G!Q$s(TUVWTW./X(  #a!C%'( 2c6!8aeAg-( 01uQw3q/@( CDQ$s(1a4-RSTURUVWRWX[R[\]_`\`R`B`cdfgcghkcklmoplpcpBp(  quQw3q 1a48QqS>1( 47( :=( @AbeAg"Q>N( qtAv "Q!( "#2a!B$q& 0!( 36"Q,2a2G!( JK3qSTUXSXYZSZIZ]^_a]abc]cIcfghjfjklflIl!( " 2a!B$q&#( ""R%'AbDF*#( "./BsF1HquQw,>,D#( "HIbfQhQRSVQVWXQXFX[^F^#( $ 2a!B$q&%( $"S&(QsU1W,quQw63q@%( $DE"S&(QsUSTWBTWXY\W\]^W^B^abcfafghahBhklmokopqkqBqtuvxtxyztzBz%( & #a!C%''( &2c6!8aeAg-'( &12r#vax!C%'/A'( &DES57QsUSTWCT'( ( quQw3q 1a48QqS>1)( (562c6!8aeAg3E)( (HIuQwQRSVQVWXQXGX)( * 3q1S57+( *!D1S57NQsU1W4q!tCx1~E+( , 1SA !Ac !Q$&Ac!Q$6-( ,9;-( ,=@-( ,CF-( ,HKSy-( . /( ./( .$/( ./( ."Ci/( .+.c)9/( .=@4#:/( .ILt/( .NQ/( 0 1( 01( 09s?1( 0&)D3J1( 014c 1( 0<=S57QsU1W;Ls;RUX;X1( 2 33( 23( 23( 2C%'AcE!G+3( 2.13( 2363( 299( 8@C9( 8EH9( 8JMr9( : c ;( :9s?;( :%(;( :+.;( :03;( :69D2I;( :@CSy;( :KNPQRUPUVWPW-Z[\_Z_`aZaJa;( < =( <T=( <2c6!8aeAg%=( <),t=( <.13Y^=( <?BTCZ=( <JMs=( > aeAg #a#%?( >(+?( >-/sF?( >58?( >:;3q1S579J?( @ #a!C%'CA( @#&A( @(+A( @-.cE!GaeAg,=A( @ABr#vax!C%PQ'?QA( B C( BC( B4C( BCAcE!G+C( B./sU1WquQw->C( BACC( BEGC( BIKC( BNPCC( D E( DdRiE( DE( D"BhE( D*-s E( D69DE( DG( FADtG( FGJTG( FLOG( H #a!C%'I( HI( H!$c I( H,/9r>I( H=@4#:I( HHKQsUSTW}WXY\W\]^W^G^I( J K( JK( J4#:#a'!C%'1K( J47K( J9@VK( JGHRVAXPQRUPUVWPWEWK( L 3q1S57M( LcE!GaeAg-M( L03M( L69DM( L;>M( L@AuQw3q?PM( N 3q1S57O( FR QQQ #Q(+Q03Q8;Q@CQ QQQ #Q(+Q03Q8;Q@CQ QQQ!$Q),Q14Q9QS TAq!R 1aAq!"4aAq1vQ ! 3Qb4719n9M8M NaAq1vQ ! 3QR$q'!)^8L M$q'!Q!Aa1 4a$rDGAI~:N9N O$q'!Q!Aa1 4ab4719n9M N  P PP PrcPtttzttzkDsJyr)rrrrrrtest_N1r+s q"ur1u}  rctdd}t|dz|z dzdkDdusJt|dz|z dzdkDdusJy)NrTr[rrrFr rrps rtest_N2r.sV$A q!tax!|a D (( ( q!tax!|a E )) )rc tdd}tttd|t|dt |dksJy)NrTr[r2r)r rrOrrXrps rtest_N3r0s9$A s2b!9bAh'Q! 55 5rcttdd\}}td|dzzd|dzzkD||kD|dkDzdusJy)NrTr[rrr-rcs rtest_N4r2sH 5t $DAq qAv!Q$Q1q5 1 2d :: :rctdd\}}}t||dzz||dzzkD||kD|dkDz|dkDzdusJy)Nzx y kTr[rrr-)rrrs rtest_N5r4sSgD)GAq! qAv!Q$Q1q5 1QU ; < DD Drctdd\}}}}t|||zz|||zzkD||kD|dkDz|dkDz|dkDzdusJy)Nzx y k nTr[rr-)rrrrs rtest_N6r6s`.JAq!Q qAv!Q$Q1q5 1QU ;q1u E F$ NN Nrcbtdd\}}t|dkD|dkD||dz k\zdusJy)NrTr[rrr-rcs rtest_N7r8s< 5t $DAq q1uq1ua!e, - 55 5rctdd\}}}tt||t||z||k\||k\z||k\zsJy)NrTr[)r rrr(s rtest_N8r:sTgD)GAq! r!Qx"Q("Q16"a1f- // /rc td}tt|dz dkDtjt t t dddt dtdk(sJy) Nrrrrr2FTr)rrrXrrrhrgrrps rtest_N9r<s[s A CAJN177 3uXrc2uVZ=[-5aT-B8D DD Drc td}|dz |dz z|dz z|dz z|dz z}tt|dktjt t t dd d t ddd d t ddd d k(sJy) NrrrrrrrrT)rrrrrrhrgr)rrYs rtest_N10r>ss A QQQQ'Q/A F1IM!'' 2eHbS!TSWrrrts rtest_N15rKs? 5>DAq S1c!fqj!A% &! +Q8rctd\}}t|dzt|dz dzzt|dzzdk|tj y)NrHrrr)r rr>rBrrrIs rtest_N16rMsF 5>DAq adc!fqj1_ %c!fai /! 3Q@rctttzdkDttz dkfttftttkk(sJyr)rrrrXrrrtest_N17rOs: QUQYA *QF 3A CC Crctdtzddtzf}tt|j |j tdk(sJy)Nrr rr)rrr8rdotHMs rtest_O1rUsCAr1Q3 A quuQSSz" #tBx // /rcrtdjtdtdgdgdggk(sJy)N)rrr)rrrrr r)rcrossrrrtest_O2rX$sI *  # #F9$5 6&2$CE$CD#BG;H HH HrctdNzfThe vector module has no way of representing vectors symbolically (without respect to a basis)r"rrrtest_O3r[+ = >>rcddlm}m}|d}|}|j\}}}|j\}}} |||z| zz|||z| zdzzz||dz| dzzzz} |j | j d|dzz|dzz| zd|z| dzzz|z||z|zzd|z|dzz| dzz|| zz |zzk(sJy)Nr) CoordSys3DDelrrrr ) sympy.vectorr^r_ base_vectors base_scalarsrWr) r^r_rdeloprrrrrrFs rtest_O4re1s,3A EEannGAq!annGAq! 1Q3q5 A!AzN"QA1 %66A ;;q>   R1WQT\!^ac!Q$h%>$AAaCE$IQqSQRTUQUXVWYZVZ]]^_`]`M`bcLc$c cc crctdrZr"rrrtest_O5rg:r\rc Htgdtgdtgdg}t|tdgdgdggttddgtd dgtd dggttd d gtd dgtddgggk(sJy)N)rrr)rrr)rrrrrrrr?iiiaiiiY)rrr)Ls rtest_O10rkCs  F9-vi/@AA q>f c c c&#$% Ryk Ryk bzl&,-% s|n }o s|n&./ 0 0 0 0rcbtdddjdtddddgk(sJy)Nrc ||z Srr)rrs rrqztest_P1..Ts QUrr2rr)rdiagonalrrrtest_P1roSs; !Q* + 4 4R 8F 1r2h=  rctgdgdgdg}|jd|jd|tddgddggk(sJy)N)rrr)rrr)rrrrrrr)rrow_delcol_delrSs rtest_P2rsXsT 9i01AIIaLIIaL AA ! !! !rc Dtgdgdgdgdg}|ddddf}|d }|}|dd d df }|d  }t||ggj}| |g||ggttfd t}|tgd gdgdgdgdgdgdgk(sJy)N)rrrr)rrrr)rrrr)r*r,rrrr))rrr)rrrr))rr)rrctSr)r)rowssrrqztest_P3..ns {40r)rrr)irGrrr)iiirrrv)rrrrrzr|)rrrrr{rG)rrrrrz)rrurrviry)rrTr~ ValueError) rA11A12A21A221A222A22r rxs @rtest_P3r`s  A AaC1H+C !C C ac1Q3hKrrBr rSs rtest_P6rsQQ a&#a&!# $A 66!Q<6SVGc!fW#5$'FSVG#4#67 77 7rc tttggttgdgdgztgdgdgzz}|tttdzztdtzdz zztdtzd z ztd tzd zzztd tzd zztdtzdz zzggk(sJy)N)rrr)rrr)rr)rFrryrrrrrrrrrrr)rrrrrSs rtest_P7rsAx &)   & (5(7!8 8 9A AEQ!a[0!QqS1W+1Q38 2LAaC"H1Q38 4678 88 8rctddtzgdtzdgg}|jtjdk(sJy)Nrr rr)ordr)rrnormrInfinityrSs rtest_P8rsCBqD Aq  A 66ajj6 !Q && &rc @tdd\}}}t|||zz d|z d|z gd|z |||zz d|z gd|z d|z |||zz gg}t|jd|dz|dzz|dzzt |t |zt |zz k(sJy)Na b cT)nonzerorfror)r rr]rrXrr^rrTs rtest_P9rsgt,GAq!AaC!A#qs#1a1gqs#1ac1ac7#% &A !&&- QTAqD[1a4%7#a&Q-A:N$O OO Orc tdddtzzgtddtzz dgg}|jtdtddtzzgddtzzdggk(sJy)Nrrrrrr)rrrrRrSs rtest_P10rswA!G 1qs7Q! "A 33&1aAaCj/qs7A,() )) )rctttgdttzggjdtdzdz z ttdgdtz ttz ggzk(sJyNrrr2)rrrinvrrrtest_P11rsr Aq6qs8 !ce1a4!8 fq"g?A!tQqSk>K7L(L LL Lrc tttgdttzggjd}t t |}t |||z dt |tttzt gdtggdk(sJy)NrADJF)evaluater2)rrrrr^tuplerR)rTrs rtest_P11_workaroundrs AAaC!"&&u-A E!H A !QqS5 )VAv 1qb q ?8#.% %% %rctdtt}tdtt}tdtt}t||gttt|gg}t |j t|j d|j z|z|j zgttt|j ggk(sJy)Nrrrr2)rrrrrr)rrrr s rtest_P12rs ua #C ua #C ua #Cc3Z A&,. /A !## +CEEz#~cee7K/L0:1a0@#%%/H/J#K KK Krc tdtdz tdz gtdz tdzdtzz dztdzdtzz dz gtdz tdzdz dtdzzdtzz dzgg}|j\}}}t|tgdtdz dd gtdz tdz dggk(sJt|tdtdz tdz gd d td z gd d td z ggk(sJy) Nrrrrrrr)rrrrrrr)rrLUdecompositionro)rTrjU_s rtest_P13rs(Aq1u=Q1qs QadQqSj1n=Q1q1a42a4"(<=? @A!GAq! A;&"3#$q5!"3#$q5!a%"3"56 66 6 A;&1a!eQU"3#$!QU"3#$!QU"3"56 66 6rctgdgdgdgdg}|j\}}|tgdgdgdgdgk(sJy) N)rrrrr)rrrrr)rrrrr)rrrrr)rrr2rr)rrrrr2)rrrrr)rrrrr)rrref)rTrrs rtest_P14rsW! "A 668DAq ))))+, ,, ,rcVtgdgdgdg}|jdk(sJy)N)r2rrr )rr rr)rrrrr)rrankrSs rtest_P15rs-! "A 668q==rctdtdzdgdtdzdtdzgg}|jdk(sJy)Nrrrrrr)rr8rrSs rtest_P16rsI47A47BtAwJ') *A 668q==rc tdd}ttd|ztd|zgddt|dzz zt|zddt|dzzz t|zgg}|j dk(sJy)NrJTr[rr)r rrBr>r)rJrTs rtest_P17rs$A QqS3qs8 AQO SV #a!SVQY-&7Q%?@BCA 668q==rctgdgdgdg}|jtdgdgdgdggtdgdgdgdgggk(sJy) N)rrr r)r rrr)r2rrrrrrr)r nullspacerSs rtest_P18rsy  !A ;;=VaS&'S&'S&'S%*+$aS&(T&'S&'S%*+ , ,, ,rc ptd}tgd|tttg|dztdztdztdzg|dztdztdztdzgg}|j |dztdzztz|dztdzztzz |dztztdzzz |dztztdzzz|dztdzztzz|dztztdzzz |dztdzztzz |dztdzztzz|dztztdzzz|dztztdzzz |dztdzztzz |dztztdzzz|tdzztdzzz|tdzztdzzz |tdzztdzzz |tdzztdzzz|tdzztdzzz|tdzztdzzz tdztdzztzz tdztztdzzztdztdzztzztdztztdzzz ttdzztdzzz ttdzztdzzzk(sJy)Nw)rrrrrr)r rrrrdet)rrTs rtest_P19rs A%qQ1%Aq!tQT1a4(Aq!tQT1a4(* +A 557q!tAqDy{q!tAqDy{2QT!VAqD[@1a46!Q$;NAad1%'(!tAvad{356T!Q$Yq[ACDa419Q;OAa1%'(!tAvad{356T!Q$Yq[ACDa46!Q$;O1a41%()Avad{367q!tVAqD[ADEQT6!Q$;O1a41 %()Avad{ 367T!Q$Yq[ ADEa46!Q$; O Aad1 % ()!tAvad{ 3 67q!tVAqD[ A DEQT6!Q$; O  rctd)Nz'Matrix minimal polynomial not supportedr"rrrtest_P20rs G HHrctgdgdgdg}|jtjtdzdtdzzz dtzz dzk(sJy)N)rrr)r rr)rrrErrrr)rcharpolyras_exprrSs rtest_P21rsX  A ::a= "adQq!tVmac&9A&= == =rcpd}dtz t|z}|jt dz|ik(sJy)Nrr)rr eigenvals)rrTs rtest_P22r!s6 A QAA ;;=aR!VQK '' 'rc tgdgdgdgdgdg}|jtddtddtd dtd dtd dik(sJy) N)rrrrr)rrrrr)rrrrr)rrrrr)rrrrr1r23z sqrt(3) + 2z -sqrt(3) + 2rrrrSs rtest_P23r'sn  A ;;= # # # -! .1   rctgdgdgdgdgdgdgdgdg}|jtd d td d td d td dtdd tdd tdd ik(sJy)Nc@rFrrqrr}rFrrrr=1rrrrrrrrv;r|rFrrrrrr}rrvrrrrrrFrrrrcqrrrr|rrrr0rz10*sqrt(10405)z100*sqrt(26) + 5101000rz-100*sqrt(26) + 510z-10*sqrt(10405)1020rrSs rtest_P24r6s??@@@@@@B CA ;;= # Q  & 1  ! a & 1  rc ttgdgdgdgdgdgdgdgdg}t|jd }tgd }t ||D]\}}t ||z d krJy) NrrrrrrrrT)multiple)gMcdr&gyP9?@@rgx5F27ߏ@g@gMcd@g-q=)rrsortedrziprX)MFev_1ev_2rrs rtest_P25rIs 6CCCCCCCCE F GB ",,,- .D  8 9DD$"11q5zE!!!"rc td\}}}}}t| | | | | ddddg gdgdgdgdgdgdgd gd g }|jd td dtddik(sJy)Nza0 a1 a2 a3 a4r) rrrrrrrrr) rrrrrrrrr) rrrrrrrrr) rrrrrrrrr) rrrrrr2r2rr) rrrrrrrrr) rrrrrrrr2r2) rrrrrrrrrF)error_when_incompletez-1/2 - sqrt(3)*I/2rz-1/2 + sqrt(3)*I/2)r rrr)a0a1a2a3a4rTs rtest_P26r\s !12BBB"rcB3bS1q1q999999999; A ? # ? # ? #  Q AAqy!Q* +   Q AAqy!Q* +      rctdNzQGeneralized eigenvectors not supported https://github.com/sympy/sympy/issues/5293r"rrrtest_P28r , --rctdrr"rrrtest_P29rrrc tgdgdgdgdgdg}|j\}}|tgdgdgdgd gd gk(sJy) N)rrrrr2)rrr rr)rrr2rr )rr2rrr)rr2rr2r)r2rrrr)rrrrrr)rrrrrr)r jordan_form)rTrJs rtest_P30r s^##### % &A ==?DAq ((((( *+ ++ +rctd)Nz!Smith normal form not implementedr"rrrtest_P31rs A BBrc (tddgddgg}t|jtj tt tdzt t dzgt t dzt tdzggk(sJy)Nrr r)rr.rewriter>rorrBrSs rtest_P32rsBA A q6>>#  ' ' )VaAh3q6 5J67Ah3q65J5L.M MM Mrc td\}}tgddddd|zggddd|zd|dzzdgg}t||zjtj tdd |zd t ||zz|z zd |z|zd t ||zzz dt ||zz|z d|z zgdd t ||zzdz d |zt ||zzd |zzdt ||zzgddt ||zz|z d|z z d t ||zzd zt ||z|z gddt ||zzd|zt ||zzt ||zggk(sJy) Nzw t)rrrrrr)rrrrr rrrrrr)r rr.rr>rrB)rrJrTs rtest_P33rs 5>DAq&qq!A#&&BqD!AqD&A&( )A qs8  C ' ' )V BqD1S1X:a< !A#a%!C!H*"4bQqSk!mac6IJ 3qs8aAc!A#h1!4qQqSzJ Qs1Q3xZ\AaC RAaC[1_s1Q3xzJ Bs1Q3xK!C!H QqSJ 5L.M MM Mrctdd\}}}t|dddddgd|ddddgdd|dddgddd|ddgdddd|dgddddd|gg}t|tt|t|ddddgdt|ddddgddt|dddgdddt|t|t| dz gddddt|t|gdddddt|ggk(sJy)NrTr[rrr)r rrBr>rs rtest_P34rs3gD)GAq!Aq!Q"Aq!Q"Aq!Q"Aq!Q"Aq!Q"Aq!Q" $ %A q6Vc!fc!faAq9QAq!4CFAq!4As1vs1vAwqyAAq#a&#a&9Aq!SV4 67 77 7rcrtdz tgdgdgdgz}t|tdk(sJy)Nr)rrr)rrr)rrrr)rrrBrrSs rtest_P35rs< 1VY ! !A q6SV  rcdtddgddgg}t|tddgddggk(sJy)Nrrrrrr)rr8rSs rtest_P36rsMQB A 7fq!f !f&' '' 'rctgdgdgdg}|tjztdtdddggdgdgk(sJy)N)rrrrrr)rrrrrr)rrr}rrSs rtest_P37rsU  A aff9AaGQ(7(7 9: :: :rctgdgdgdg}tt5|tjzJ#1swYyxYw)Nr)rrr)rr~AssertionErrorrr}rSs rtest_P38rsE iA   166  t s :Actd)zQ M=Matrix([ [1, 1], [2, 2], [3, 3]]) M.SVD() z,Singular value decomposition not implementedr"rrrtest_P39r!s L MMrc $tdd\}}t|t|z|t|zg}|j t||gtt|| t|zgt||t|zggk(sJy)NrHTr[)r rr>rBjacobian)rrJrTs rtest_P40r$s 5t $DAq#a&!CF(#$A ::faVn %#a&1"SV)1D25a&1SV81D1F*G GG Grc tdd\}}t|dzt|z||ftdt|zd|zt |zgd|zt |z|dz t|zggk(sJy)NrHTr[r)r rPrBrr>rIs rtest_P41r&s 5t $DAq 1a4A;!u %Ac!fH!CF 1K23A#c!f*q!teCFl1K1M*N NN Nrctttttgtj dk(sJyr)rQr>rrBrorrrtest_P42r( s- c!fc!f%q ) 2 2 4 99 9rc d}tdd\}}t|t|z|t|zg}||||gtt|| t|zgt||t|zggk(sJy)Nct|Dcgc]}|j|jc}jScc}wrrr r~)rTYvs r __my_jacobianztest_P43..__my_jacobians-A.qqvvay{{./111.s"=rHTr[)r rr>rB)r.rrJrTs rtest_P43r/s2 5t $DAq#a&!CF(#$A Aa5 !Vc!fqbQi-@.1!fqQx-@-B&C CC Crc d}tdd\}}||dzt|z||ftdt|zd|zt|zgd|zt|z|dz t|zggk(sJy)Nc t|Dcgc]}t||c}}t|Dcgc]}|jj|c}Scc}wcc}wrr+)rr,r-Vs r __my_hessianztest_P44..__my_hessiansG *1DAJ* +Q/ /00+/s A"ArHTr[r)r rBrr>)r3rrJs rtest_P44r4s1 5t $DAq 1SV aV ,/0Qx1Q3s1v:,F-.qSQZ!Q$s1v,F8H1I II Ircd}|ttttgtjdk(sJy)Nc ttdt|Dcgc]+}t|jj t |-c}}|j Scc}wr)rrlenr~r rr)r,r-rrTs r__my_wronskianz test_P45..__my_wronskian"sH E!SV4DEqF1IKK$$Q*E FFs0A$r)r>rrBro)r8s rtest_P45r9!s4 3q63q6*A . 7 7 9Q >> >rc tddd\}}}td|d}t||dft||df|d|dz f|z z dz|d|dz f}|jt d) Nzi j nTrxnrrrzUnknown result)r rrrr#)rrrr;Sms rtest_R1r=*sgtd;GAq! dAq !B bAhR1X1a!e}5a77! ;aAE] KBGGI . //rc6td\}}tddd\}}td|d}td|d}t||df|||dfzz |z d z|d|dz f}t||}t||}t ||f||ft j y) Nzm bzi nTrr;rynrrr)r rrr rrr) rr^rrr;r?rf1f2s rtest_R2rB4s 5>DAq 5$ 6DAq dAq !B dAq !B R1X"QT( "Q & *Q1q5M:A aB aB b"X1vagg.rctddd\}}d|ztd|z|dzz}t||dtf}|j }|j }|d|ztd|z|zk(sJy)Nn kTrr2rr)r r"rrrri)rrskr<r~T2s rtest_R3rGAs~ 5$ 6DAq 'hqsA&* *B R!Q B  A B "q!A#q)) )) )rctd)NzIndefinite sum not supportedr"rrrtest_R4rIMs& < ==rcbtddd\}}}}}d|zt||z||zt||z||zzt||z||zzz}t||dtf}|j }|t ||z|zt |t |zt |zz k(sJy)Nz a b c n kTrr2r)r r"rrrr#)rr^rrrrEr<r~s rtest_R5rKcsKEMAq!Q 'HQUAE*a!eQU+,,4QUAE,BC DB R!Q B  A !A#a% )A,y|";IaL"HI II Irctddd\}}td|dzd}t||df||dz dfz |d|dzf}|j|d ||dzdfzk(sJy) NrDTrgnrrr)rr)r rrr)rrrMr<s rtest_R6rNms 5$ 6DAq dAE1 %B R1X1q5!8 $q!QUm 4B 779D Bq1uaxL0 00 0rctddd\}}t|dz|d|fj}|j|dz|dzdzzdz k(sJy)NrDTrrrrr)r rrr])rrr~s rtest_R7rPts\ 5$ 6DAq AqD!Aa A 88:Aq1uqj* ** *rctddd\}}t|dzt||z|d|f}|j}|j ||dzzd|dz zzk(sJyNrDTrrr)r rr"rrirrr<r~s rtest_R8rTysi 5$ 6DAq QT(1a. 1a) ,B  A ::<1a!e9QQZ/ // /rctddd\}}tt||dz |z |d|dzf}|jj d|dzzdz |dzz k(sJyNrDTrrr)r rr"rro)rrr<s rtest_R9rWsl 5$ 6DAq XaQ  !Aq!a%= 1B 779   AAJNQU#; ;; ;rctddd\}}}}tt||t|||z z|d|f}|j}|j j t }|t ||zt |t ||z|z zz k(sJ|t||z|j t k(sJ|j t}|t||z|k(sJy)Nzn m r kTrr)r rr"rrirr#)rrrrr<r~rFT3s rtest_R10rZsD4@JAq!Q Xa^HQA. .Aq :B  A   i (B 1q5!9Q< !a%!)0D#DE EE E !a%#++I6 66 6 H B !a%# ## #rctddd\}}t||t|z}t||d|f}|j }|td|zk(sJy)NrDTrrr)r r"r&rr)rrrEr<r~s rtest_R11r\s] 5$ 6DAq !Q ! $B R!Q B  A !A#  rctddd\}}tt|dz|d|f}|j}|t|t|dzzk(sJy)NrDTrrrr)r rr&rrSs rtest_R12r^sY 5$ 6DAq Yq\1_q!Qi (B  A ! Yq1u-- -- -rc8tddd\}}tt|tz|d|f}|j }|j t tdz dz ttd|zdzzdz dttdz zz z k(sJyrV)r rrBrrror?r>rSs rtest_R13r`s 5$ 6DAq S1X1ay !B  A ::<3qs8A:AqsQwKM(:Ac!A#hJ(GG GG Grctddd\}}ttd|zdz tz|d|f}|j }|j t|tzdzttz k(sJyrR)r rrBrrrorSs rtest_R14rbsk 5$ 6DAq S!A#'1 1ay )B  A ::<3qs8Q;s1v- -- -rc tddd\}}tt||z ||dt|dz f}|j }|j t |dzk(sJy)NrDTrrrr)r rr"r5rror&rSs rtest_R15rdsb 5$ 6DAq Xa!eQ !Qac !3 4B  A ::<9QU+ ++ +rctddd}td|dzz d|dzz z|dtf}|jt dt dzdz zk(sJy)NrTrrrrr)r rrrryr)rr<s rtest_R16rfs[TD1A Qq!tVa1f_q!Rj )B 779Q"a%') )) )rc tddd}tttd|dzz d|dzz z|dtfdz dksJy) NrTrrrrg3@V瞯<)r rXfloatrrrs rtest_R17rksZTD1A uS1a4!AqD&1a*56#$ %', -- -rctddd}tdd|z|dzzz |dtf}|j}|j t ddz dz t dzdz zk(sJy)NrTrrrr)r rrrror/rrr<r~s rtest_R18rnslTD1A Q1QT ]Q2J 'B  A ::rrSs rtest_R20rrs 5$ 6DAq Xa1 1bz *B  A ::<1qs8C1QK/1AAJqL@ @@ @rc tddd}tdt||dzzt|t|dzzzz |dtf}|j }|j dk(sJy)NrTrr)r rr8rrrorms rtest_R21rtshTD1A QQAY47T!a%[#89 :Q2J GB  A ::<1  rc tddd\}}ttt|t|dzt|d|zz zz ttz |zzttz||z zz|d|zt f|dt f}|j }|dttdztdzzdtdzzz dtztzz dzz k(sJy) NrDTrrrr2rr)r rr#rrrrr8rSs rtest_R23rvs 5$ 6DAq S)A, ! a !ac'0B BCA#Q3!a%.)+,ac2,8:;Q EB  A 41QT AadF*QqSU2Q677 77 7rctddd\}}tt|d|zdz z |d|f|dtf}|j t dz k(sJy)Nzm kTrrr)r rrrrr)rrr<s rtest_R24rxsW 5$ 6DAq WQ!a[1a) ,q!Rj 9B 7791  rctddd}tt|dz |ddf}|jj dt dzt dzzdz k(sJy) NrTrrrrii)r rrHrror8r)rPrs rtest_S1r{s]TD1A qsaAY 'B 779   3tAw;r1u#4T#9 99 9rc|tddd\}}t||d|fjt|k(sJyNrDTrr)r rrr#rrs rtest_S2r s< 5$ 6DAq 1q!Qi % % '9Q< 77 7rctddd\}}tt|z|d|fjj t||dzzdz zk(sJyrV)r rrrror~s rtest_S3rsV 5$ 6DAq 1a4!Q # ( ( * 3 3 5QAYq[9I II Irctddd\}}tdd|z z|d|dz fjj|k(sJyr})r rrror~s rtest_S4rsM 5$ 6DAq 1qs7Q1a4L ) . . 0 9 9 ;q @@ @rctddd\}}td|zdz d|zz |d|fjjt |t j zttt |dzzz k(sJyrR) r rr gammasimprHrr}r8rr~s rtest_S5rs{ 5$ 6DAq QqS1WqsOaAY / 4 4 6 @ @ B !aff* tBxa!e 4 5 67 6rctddd\}}ttdzdtzt|tz|z zz dz|d|dz fj j td|zzdz tdzdz z k(sJyrR)r rrr>rrror~s rtest_S6rs 5$ 6DAq AqD!A#c!B$q&k/)A-1a!e} = B B D M M OAaC1 q!tax( )* )rctddd}t|dzdz |dzdzz |dtf}|j}|j t ddk(sJy)NrTrrrr)r rrrrorrrzr~s rtest_S7r&s\TD1A !Q$(QTAX&Ar 3B  A ::<1Q7 "" "rctddd}tddd|zdzz z |dtf}|j}|j dt z k(sJy)NrTrrr)r rrrrorrs rtest_S8r.sUTD1A Q!axZ!Q ,B  A ::<1R4  rctddd}tdd|dzzd|zdz z z|dtf}|j}|j t dk(sJy)NrTrrr2r)r rrrror8rs rtest_S9r7saTD1A bAE]AaC!G,,q!Rj 9B  A ::<47 "" "rctddd}t||dzzdztz||dzzdztz z |dtf}|j }|j dk(sJy)NrTrrrr2)r rrrrrors rtest_S10rAshTD1A !QU)a-!#aQi!ma&781a* EB  A ::<2  rctddtz ztztttk(sJtdt t z t dzz t dt jk(sJyr)r`rrrr>rrr}rrrtest_T1rJsQ !ac'Aq" % ** * !c!f*ad"Aq )QVV 33 3rcntdtzdtzzdtz zttdk(sJy)Nrrr)r`rrrrrtest_T2rOs. !Q$A+1%q" - 22 2rctttttttzz ttdk(sJyr)r`r/rrBrrrrtest_T3rSs. QQ#a&)1b 1Q 66 6rc ttttt ztt tdtdzztdzz zz ttz tz tttd k(sJy)Nr rr)r`r.rrrrrtest_T4rWsq #aQBiaR3r!Q$wA+?!?@Aq6 %),Q 00 0rcjttttzttttztdzz dzztttdzdttdtdzzttzzzz ttt ddk(sJy)Nrrr)r`rr/r.rrrrrtest_T5r\s !CF(3qQx!Q$/223s1a4c#a1fSVm,--<.8/40012B 8; QqT7<<1  $$ $ 88:  rctd}tttt|tt tf}|j ttt ttt ttzt|tt|ttzz k(sJy)Nr)r rrrrrr r )rr~s rtest_U6rs~ A!A$AaD!A$(A 66!9 !A$ 1Q4##a!gj1q.A&AA CC CrcNtdd\}}tt||y)Nzp tTr[)r r r)rYrJs rtest_U7rs! 5t $DAq  1aMrctdd\}}t||z|z}t||z ||| t||zzdz|t||zzdzz k(sJy)NrTr[r)r r>rNrB)rreqs rtest_U8rsg 5t $DAq QqSAB RA A2c!A#h;?Qs1Q3xZ!^"D DD Drc tdd\}}tt|||tt|||z}|jt||t |dz|dzz}|j j }|||zttt ||||dz|dzzzdzk(sJy)NrTr[r) r r rrrrr]r r )rrsus2s3s rtest_U9rs 5t $DAq a1gq D1a!, ,B 1a!AqD1a4K. )B    B !a%j1q11adQTkBB1D DD Drcttdzdztdzdz tdzzz tdtddk(sJy)Nrrrrr2r)rbrrrrrtest_U10rs= AqD1H1q1q512Ar :aAh FF Frctrr"rrrtest_U11rs  rctd)Nz0External diff of differential form not supportedr"rrrtest_U12rs : <Aa#O OO Orcttt tdkfttdk\ftttdz dz tdkftdzdz dfk(sJy)NrrT)rr9rrrrtest_V2r! s_ YAE{QQK8! AaQ'!Q$q&$ 8 99 9rctdtdzdzz tjjdtdzdzz k(sJy)Nrrr)rrr rorrrtest_V3r& sA Q1q\! $ ) ) + 4 4 6!QTAX, FF Frctdtztddtzzz ttdtzt dz k(sJy)Nrrr)rrr8r0r/rrrtest_V4r* s= QT$q1a4x.(! ,ad CF0B BB Brctdtzdz dzdtzdz tddzz tjddtzzdtdzzz ddtzdz tddzzz k(sJy) Nrrrrrir-)rrrrorrrtest_V5r. s qsQwlAaC!Ga1g#55q 9 B B D 2a4Z"QT' !AqsQw1a&8$8 9 :; :rc ztddtttzzdtt tzzz z tt dt dtttzzt dz t dtttzzt dzz zdtzz k(sJy)Nrrrrr)rr.rrr8r/rrrtest_V6r5 s a3qs8aQBqD k12A 6$r( #ac( T"X% &Qs1Q3xZ$r(-B)C C;EFHd;L LM Lrc$tttdzttdzz }|j tt ddzttdzdttzz zdt tzdz zk(sJy)Nrrrr)rr2rr1rorr3r1s rtest_V7r< sj 47A:d1gqj( )B ;;=AaAhJa!QtAwY)??!DG)A+M MM Mrctd)Nz'Integrate with assumption not supportedr"rrr test_V8_V9rA s 1 33rctdddttzzdttzzz tt dt tdz zdzdz k(sJy)Nrrrr)rr>rrBr/rCrrrtest_V10rR sO QAc!fH qQx/0! 4Ac!A#hJN8KA8M MM Mrc.tdddttzzdttzzz t}t |}t |dt ttdz dzttdz dzz tddzk(sJy)NrrrT)forcerr) rr>rrBr]rnr/rCrrr2s rtest_V11rV s~ 1a!CF(lQs1vX-. 2B B r & #ac(Q,QqSA.1a8 9 :; :rctdddttzzdttzzz t}|dt tdz dzz k(sJy)Nrrrrr2r)rr>rrBrCrs rtest_V12r] sK 1a!CF(lQs1vX-. 2B S1X\" "" "rc  tdddttzzdttzzz t}|j dt dzt t ddttdz zdzzdz zdz k(sJy)Nrrrrrr)rr>rrBror8r=rCrs rtest_V13rb sw 1a!CF(lQs1vX-. 2B ;;=Ad2hJtDHaAaCj1n,Eb,H'II"L LL Lrc tttttdztdzz t}|j ttttdztdzz zttttzzzttttz zz dtzz k(sJyr,)rr/rXrrrors rtest_V14rj s 3s1a4!Q$;'(! ,B KKMQs3q!tq!t|#455AE l+-.s1q5z\:<=aC@ @A @rc tttttz zt}t |ttztdztdzztttz zzdz z dk(sJyNrr)rrr;rrors rtest_V15rs sY 1T!A#Y; "B B!A#A1d1Q3i 77:: ;q @@ @rcttdtztdtzzttdtzt dtzzdz t dtzt dtzzdz z k(sJy)Nrrrrr)rr>rrvrBrwrrrtest_V16rx sk S1Xb1g%q )R!WS1X-=a-?2ac7RPQRSPSWCTVXBX-X XX Xrc tttttt tztttt ttzz ttdzt tdzz z t}t |ttt tz ttt tzz dz z dk(sJyr)rr rrrrors rtest_V17r} s D1qM!A$&qT$qtQ-'(*+A$'AaD!G*;=>? AB B!A$1+!qt 4Q66 71 << rrrrrrtest_W6r sN T!aAaCj.)!+aAb!Hrc!e-D Ea PP Prctdd}ttttdz|dzzz tt tf}|j t t| z|z k(sJy)NrTror)r rr>rrrorr.)rrs rtest_W7r sYd#A 3q61a4!Q$;'!bS" 6B ;;=BsA2wJqL (( (rctd)Nz1Integrate with assumption 0 < a < 1 not supportedr"rrrtest_W8r s ; ==rcJtdtdzzdtztdzztdzztdzzz tt tf}|j}|dtzt t d dz dzt t ddz dzzzk(sJy) Nrrrrrr rg?)rrrrrr8rs rtest_W9r s 1QT61q51a4rrrrrtest_W15r# s> Sq]3qtAv;.Aq :a2h FF Frctdtzdztdtztdtztddftddk(sJy)Nrrrr2r#)rrrZrrrrtest_W16r& sM a!eaZ a 33M!Q4GGQZ !$%b"I .. .rctdd\}}tt| tzt d|tzztdt fd|t |dz|dzz dzzz k(sJy)Na bTrorrr)r rr.rrDrr8r]s rtest_W17r) so 54 (DAq S!AYwq!A#.Q  !1T!Q$q!t)a-%8#8 9 :: :rcttdttz dztdtfddtzz k(sJyr)rrDrrrrrrtest_W18r+ s4 gamAo)Aq": 6!QrT( BB Brc ttttddt dtzzztdt ft ddz dz k(sJy)Nrrrr)rrvrrDr8rr>rrrtest_W19r- sH RU71aQqS k22Q2J ?CFQJPQ> QQ Qrc :ttdztddtdzz ztddftdz dz t ddz t ddz ztdz dz dt dzz t ddzzd zt dzd z zk(sJy) NrrrrrrlrgUUUUUU'@r)rrrMrrryr/rrrtest_W20r0 s ad71aQi001a) < UF2I"c "T!WQY .!eVAX#a& 3q619 ,t 3SV ;A = > >? >rc ttttdzt ddtdzz ztddfdz dksJy)Nrrrrg5?rh)rXrrrrMrrrtest_W21r2 sP q1a41a!e9 551ayAB  $ %% %rc Ntdd\}}tttdt tdk\tdkfd}t ||t |z|d|ftd|dkfttd| ttd|zdfk(sJy)Nzt uTr[rr)rTr) r r rr9rOrr>rBr7)rJurs rtest_W22r5 s 5t $DAqq)QAFAF 34i@AA QqT#a&[1a) , AE c!Qi.3s1ay> )4012 22 2rc"tdd\}}ttttdztdzzz t||ftt tf}|j t jt |zt |zzk(sJyNr(Tror)r rrrrcollectrcancel)rr^rs rtest_W23r:# sp 54 (DAq 9Q1q!t _q!Qi81rc2, GB ::b> "rc!ebdl 22 2rctdd\}}ttttdztdzzz tt tft||f}|j t t | |zzk(sJyr7)r rrrrr8r)rr^rs r test_W23br<* sc 54 (DAq 9Q1q!t _q2#rl;aAY GB ::b>R!a[ (( (rc tdd\}}ttt|dz|dzz|ddf|ddf}|tdtdzdz z j dk(sJy)NrTr[rrrr)r rr8r0ro)rrrs rtest_W24r>1 sr 5t $DAq 9T!Q$A+.Aq :Q1I FB $q'E!H$a' ' 1 1 3q 88 8rc ^tdd\}}}tt|t|ztdt|dzt|dzzt|dzzz z |dtdz f}t||dtdz f}|t|zdz z j dk(sJy)Nza x yTr[rrr)r rrBr8rro)rrri1i2s rtest_W25rB; sgD)GAq!  As1v d1s1vqyQ23q619<<== Ar!t  B 21bd| $B AaK ! ! #q (( (rc tdd\}}ttt||dzz |ddf|ddftdd k(sJy) NrTr[rrr2rrr)r rrXrrcs rtest_W26rDF sV 5t $DAq Ys1q!t8}q!Qi8QZ !$%b"I .. .rc td\}}}tttdtd|dt|z z t|z z zftd|dt|z z zftd|f||z|zdz k(sJy)Nrrrr)r rrrr)rr^rs rtest_W27rFL sgGAq! YyQ1a!A#g!m3D,EF !1aQqSk24AY #$Q3q57 ++ +rctdd\}}tdtd||z dzz z |ddd |d zzd |d zzz d |d zzd|d zzz z|dzd|dzzz zdzt|dzzk(sJy)Nv cTr[rrrrx0rrrrrrr rcr8ra)r-rs rtest_X1rLS s 5t $DAq 1T!qsQh,''qA 6 adFBq!tG qAvqAv .AqAv > BQq!tW L MN Mrctdd\}}tdtd||z dzz z |dd}d|dzz j|dd|dz |dzz dzt|dzzk(sJy) NrHTr[rrrrrIrK)r-rs1s rtest_X2rOY s 5t $DAq $qAaC!8|$$aA 3B b!eG  A!q  )adU1a4Z!^a1g-E EE ErcTttjttjz j}t tj}|ttdzdz zdtdzzdz zt tdzzk(sJ||k(sJy)Nrrrrr)rBrrcr>rCra)rNrs rtest_X3rQ_ s{ a&--/#a&--/ ) 1 1 3B QB QT!Va1fRi'!AqD'1 11 1 8O8rc8ttttz j}|tdz dz tdzdz z t tdzzk(sJtttttz j|k(sJy)Nrrr)r/rBrrcra)rNs rtest_X4rTf sw SVAX   B !Q$q1a48#a1g- -- - vc!fQh ' ' )R // /rc td}tdd\}}}}ttt |t zt t |t zzt||tztdt fzt |dy)Nrza b c dTr[rrrI) r r rcr rrrrr)rrr^rrs rtest_X5rVl si  A.JAq!Q 4!A#?QqsV #i!A#Aq &B B arc"tddd\}}tt||ztzt|tzt|tzzz tddtdz| |zdz ||zdz zzt tdzzk(sJy)Nr(F)rscalarrrrIr)r rcr.rrar]s rtest_X6rY s 5eE :DAq 3Aqy>C!Hs1Q3x$77qA FdQBqDFQqSUN#a1g- ./ .rctdtttdz zz tddtdzddtzz z tddztdzdz z tdzd z ztd zd z z t tdzzk(sJy) Nrrrr rriri vriu)rcrr.rrarrrtest_X7r[ s 1aQ!n%q!Q /1r7Q!W3D aH4!tCx4 "#Q$u*4-/0!tG|4<>?1g4F FG Frc tdd}ttt||tdzdz ddt|tt ddzz z |tt ddzz t ddzd z z|tt ddzz t d dzd z zt |tt ddzz dz|tt ddzfzk(sJy) NrTr[rrrrIrrr)r rcr8rArrrarps rtest_X8r^ s $A 4A<r!tAv 3 d1r!Aq'z>" "a"Qq!W*nqAw%>r%A A AaG^a1g %c ) *,-q2a1g:~.AAr!Aq'z?,S T TU TrcDtttztdddtttzztdzttdzzdz ztdzttdzzdz zttdzttdzzzk(sJy)NrrrIrrrr)rcrr/rarrrtest_X9r` s 1a4qA &!aAh,Ac!fai9I*I qDQN1 + Ac!fai0+1 12 1rc ,td\}}ttt|tt ||zz|ddtt |t|z|t|zt |z zt |dzzk(sJyNzz wrrrIr rcr/r2r1rarrs rtest_X10re s 5>DAq 3tAw<#d1q5k"22A!q A QL3q6 !Ad1gId1g$5 5!Q$ ? @A @rc td\}}ttt|t ||zz|ddtt |t|z|t|zt |z zt |dzzk(sJyrbrcrds rtest_X11rg s| 5>DAq 3tAwa!e,-qQ! < QL3q6 !Ad1gId1g$5 5!Q$ ? @A @rctdd\}}}tt||zt| |zz|dd|dz |zt|z d|d|zz|dz zdz z t |dz dzzzk(sJy)Nza b xTr[rrrI)r rcr/r.ra)rr^rs rtest_X12ri sgD)GAq! 3q619S!AY&a1 5 UQJs1v qA!Ga!e#4Q#66AEA:F G HI Hrcttdtdzzdzttdtdtzt dtz ttfzk(sJy)NrrrI)rcr8rrrarrrtest_X13rk sG $qAvz"A" 2d1gai!AaC!R/6Q QQ Qrctdddtzzz tdtztztttk(ddt t t tzz tdtz ttfzk(sJy)Nrrr)rcrr"rrr8rrarrrtest_X14rm so !A!H*Xac1--QUa !$%tBxQ'7$81QqS1b'?$J KK Krz*https://github.com/sympy/sympy/issues/7164ctdd\}}tt| |z ||tf}t ||tdd|dzz d|d zz zd |dzz z d |z zt |d z|tfzk(sJy) NrTr[rrIrrrrrr )r rr.rrcra)rrJe1s rtest_X15rp s 5t $DAq 3r719q!Rj )B 2qR1 % adFQq!tVOa1f $qs *Qq2wB-@ @ AB ArcLttttzttzdddttzdzdz z t tdztdztzztdztdzzzttdzzztdzzttzk(sJy)NrrrIrrr)rcr>rrrarrrtest_X16rr s 3q1u:q1ua 0AQ 1 4D adQT!Vmad1a4i'!AqD&01a47A >5? ?@ ?rc 2tttttz t dt zddt zdz zzt dt zztdt zzzt tdt zzz t dtfk(sJy)Nr2rr) rr/rBrrrr%r#rrrrtest_X17rt s} s3q6!8}  R!GA!aL 1Q3 /AaC 8!IacN:J KaQRTVZX YY Yrc ,td}ttt t tzt dt j|zzt tdd|ztzzt|zzt|z |dtfk(sJy)Nrrrrr) rrr.rrBrrr}rrr#rrjs rtest_X18rv s{. c A sA2ws1v~  AqM#a1gail+ +AqD 01 =1bzJ KK Krctd)NzTSolve using series not supported. Inverse Taylor series expansion also not supportedr"rrrtest_X19rx s 4 55rctd)Nz'Symbolic Pade approximant not supportedr"rrrtest_X20rz s G HHrctdd}tddd}ttt| |f}|jjdk(sJ|j jj |j jdddk(sJ|j jj |j jd|d|ztz d|zz|z t|tztz|z zk(sJy ) z Test whether `fourier_series` of x periodical on the [-p, p] interval equals `- (2 p / pi) sum( (-1)^n / n sin(n pi x / p), n = 1..infinity )`. rYTror)rrrr r2N) r rranformulabnr variablesrrB)rYrrs rtest_X21r" s d#AdD1Aq1qb!*%A 44<<1   44<<  QTT^^A. 2a 77 7 44<<  QTT^^A. 2 1R2'AAbDF1H - .. .rctd)NzFourier series not supportedr"rrrtest_X22r4 s< < ==rctdd}tdd}td}tt|dz |z||\}}}|||dz|dz dzzz k(sJy) NrJTrorr[rrr)r rr>rJrrrdrs rtest_Y1rU sgd#A$A AQUAI15GAq! 1a41q51*$% %% %rctdd}tdd}td}t||dz|dz dzzz ||d }|t||dz zk(sJy) NrJTrorr[rrrro)r rr>)rJrrrs rtest_Y2r] scd#A$A A!!QTQUQJ%6"7AMA Aq1uI  rctdd}tdd}td}tt||zt||zz||d\}}}|||dzd |dzzz z k(sJy) NrJTrorr[rrrr)r rr2r1rs rtest_Y3re sqd#A$A AQqS $qs) 3QDIGAq! 1a4!AqD&=! !! !rctdd}td}ttdt|z ||d\}}}|d|z t dt|z|z z k(sJy) NrJTrorrrrr)r rrxr8r.)rJrrdrs rtest_Y4rm sad#A AAd1gI1tDGAq! !c"T!W*oa'' '' 'rc tdd}td}td}td}ttt |||d||zdt |d z t |dz z zz ||d ||i}| |dz||zz||d zz ||zt t||||d z dt| z|z z dtd |zz|z z}|d k(sJt|||d }t|||d d gi}||dztd|zzdt|zzdz td |zz||dzd zzz k(sJt|||}|jt |d z jt |dz }|ddt|d z zz t |d z zdt|dz zdz t |dz zzt|t |zzk(sJy) NrJTr[rrr,rrr)nocondsrr )r r rrr rur r r.rrrrr8r>)rJrrr,rdDYfyfs r test_Y5_Y6rt s/ $A A A A0adAq1AAaD1H"#Yq1u%5 !a%8H%H"I2J !1d 467V =A QT!A$Y1Q4!A$&j1q.A1a)HH #qb' !  BqD k!m $ 6M6 q!A$ B r1q1a&k 2B !Q$s1Q3x-!CF(*Q.BqD 91adQh<H HH H "2q! ,B IacN # + +IacN ;B Qs1q5z\ 9QU++ 3q1u: 9QU++ , Ay|   rc ltdd}tdd}td}tddtd tzt |t|zz ztdt fzz||\}}}|dtd tzt | tz|zz|z tdt fzd|z zk(sJy) NrJTrorr[rrrr2)r rrrrurr.)rJrrrdrs rtest_Y7r s d#A$A AAc2')A!G2D*D+,a*'6%6!678!=GAq! #rAgc1"Q$q&k)!+aBZ881Q3> >> >rcRtdttttk(sJyr)rrrrErrrtest_Y8r s Q1 %A 66 6rcttdtdzzttt t tt dz tdzzdz zdz k(sJy)Nrrrr)rr.rrr8rrrrtest_Y9r sQ c"QT'lAq 1 HS"a%1Q' ' ) *+ *rc 6ttttdttzzttj dt dzztdzzdzdt dzztdzzdt dzztdzzzdzz k(sJy) NrrFrrrrrr)rrXrr.rr9rrrrtest_Y10r s c!fSCF^3Q : A A C AXad]R "RU(1a4-"RU(1a4-"?""D E FG Frz*https://github.com/sympy/sympy/issues/7181ctd\}}tdd|z z ||\}}}|ttt|zzk(sJy)Nx sr)r rrr?rrrdrs rtest_Y11r sI 5>DAq  1q5 1a0GAq! 3r!t9  rctd\}}ttd||dzz ||\}}}|d|dz z t|dz zt| dz dzz k(sJy)Nrrrr)r rrDrHrs rtest_Y12r sn 5>DAqwq!}QT11a8GAq! QU E!A#J&uaRTAX6 66 6rctdNzz-transform not supportedr"rrrtest_Y13r  9 ::rctdrr"rrrtest_Y14r rrc ,td}t|tdzd|tdzzz |tzdz |t|dd|dtij tdzttdz zzdzk(sJy)Nrrrr)r rrrrors rtest_Z1r s A 1QU8a!a%j(1Q4/!3QqTaD!QqT1% ''/xzQTAq1uI5E5I JK Jrc td}t|td|tdz zd|tdz zz z |t|dd|ddidtz dtzzk(sJy)Nrrrrrrr)r rrrs rtest_Z2r sr A 1Q41Qq1uX:!AE( 23QqTAaD!QqT1;M N1uq!t|  rc td}tddz tddz z tztddz tddz z ztddz tddz ztztddz tddz zzz}t |t|tdz |tdz zz |t|dd|ddi}||k(sJy)Nrrrrr)r rr8rr)rexpectedrs rtest_Z3r s A1a$q'!)#a'1a$q'"*)<=1a$q'!)#a'aad1f)<=>H 11q5Aa!eH,-qtadAqtQ5G HC (??rc Htd}td}t|td|z|tdz zz |tdzzz d|tzz z |tdz z|d|tdz zz zd|tdz zz z |tdz zz dzz |t|dd|ddd|zz|dzzd|zz i}||tdzz|tdzzd|zz dz ztdz|dzzzd|zztz |dz dz|dzzz z dk(sJy)Nrrrrrr)r r rr)rrrs rtest_Z4r s[  A A qtAAaC(1qs83a!Q$h?!a%HA!H %q1qs8|4Qq1uX=>@ABC Q4!A$1Q4!ac'AqD.1q5!9: r rrrr#)C1C2rrf0r@ const_dictresults rtest_Z5r s2W FB AaD!Q !AaD& (3qs8 3B QqT C 377 B a5Bs1Q3xK2!8S1X"55 55 5 41q> "B1r!u~&J UZZJrN + 0 0Z^ DF aRAaC[]S1XaZ/ // /  rcLtdd}td}tt||ddt|zztd|zz }t |||\}}}||dzt t|||zdt t|||zzd|dzdzz z k(sJy)NrJTrorrr)r r rrBrr)rJrrrdrs rtest_Z6r s d#A A AaD!Q !AaD& (3qs8 3BAq)GAq! A&qtQ22 qtQ* *+-.1q\: :; :r(__doc__sympy.assumptions.askrrsympy.assumptions.refinersympy.concrete.productsr sympy.corersympy.core.evalfrsympy.core.functionr r r r r rrsympy.core.mulrsympy.core.intfuncrsympy.core.numbersrrrrrrrrsympy.core.relationalrrsympy.core.singletonrsympy.core.symbolrrr (sympy.functions.combinatorial.factorialsr!r"r#r$%sympy.functions.combinatorial.numbersr%r&r'r($sympy.functions.elementary.complexesr)r*r+r,&sympy.functions.elementary.exponentialr-r.r/%sympy.functions.elementary.hyperbolicr0r1r2r3#sympy.functions.elementary.integersr4r5(sympy.functions.elementary.miscellaneousr6r7r8$sympy.functions.elementary.piecewiser9(sympy.functions.elementary.trigonometricr:r;r<r=r>r?r@rArBrCsympy.functions.special.besselrD'sympy.functions.special.delta_functionsrE*sympy.functions.special.elliptic_integralsrFrG'sympy.functions.special.gamma_functionsrHrIsympy.functions.special.hyperrJ#sympy.functions.special.polynomialsrKrL&sympy.functions.special.zeta_functionsrMsympy.geometry.utilrNsympy.logic.boolalgrOsympy.matrices.denserPrQ!sympy.matrices.expressions.matmulrR sympy.ntheory.continued_fractionrSrsrTrtrUryrVr|sympy.ntheory.factor_rWsympy.ntheory.generaterXsympy.polys.domains.integerringrYsympy.polys.orthopolysrZsympy.polys.partfracr[sympy.polys.polytoolsr\r]r^r_sympy.series.limitsr`sympy.series.orderrasympy.series.residuesrbsympy.series.seriesrcsympy.sets.fancysetsrdsympy.sets.setsrerfrgrhsympy.simplify.combsimprisympy.simplify.hyperexpandrjsympy.simplify.powsimprkrlsympy.simplify.radsimprmsympy.simplify.simplifyrnrosympy.simplify.sqrtdenestrpsympy.simplify.trigsimprqsympy.solvers.solversrrr5rtru'sympy.functions.special.error_functionsrvrwrxrysympy.testing.pytestrzr{r|r}r~sympy.utilities.iterablesrrrsympy.matricesrrr&sympy.matrices.expressions.blockmatrixrrsympy.matrices.expressionsrrsympy.physics.quantumrsympy.polys.ringsrsympy.polys.fieldsrsympy.polys.solversrsympy.concreterrsympy.integralsrsympy.integrals.transformsrrrrrrrsympy.solvers.recurrrsympy.solvers.solvesetrrrsympy.solvers.oderr itertoolsrrsympy.series.formalrsympy.series.fourierrsympy.calculus.utilrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr rrrrr r$r'r+r0r3r5r7r:r=r?rArCrHrJrLrNrQrSrVr\r^rarcrgrirkrurwrzr~rrrrrrrrrrorrrrrrrrrrrrrrrrrrrrrrrrrrrrr r rrrrrrrrr"r%r(r*r-r/r1r3r6r8r;r=rArCrErHrJrLrNrPrRrVrYr_rardrfrirmrprsrvrxrzr|r~rrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr&r)r+r.r0r2r4r6r8r:r<r>r@rBrDrFrKrMrOrUrXr[rergrkrorsrrrrrrrrrrrrrrrrrrrrrrrrrrrrr rrrrrrrrr!r$r&r(r/r4r9r=rBrGrIrKrNrPrTrWrZr\r^r`rbrdrfrkrnrprrrtrvrxr{rrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr r rrrrrrrrrr!r#r&r)r+r-r0r2r5r:r<r>rBrDrFrLrOrQrTrVrYr[r^r`rergrirkrmrprrrtrvrxrzrrrrrrrrrrrrrrrrrrrrrrrrrs)++!#("44ZZ  EE  >CC:(((2>D/:%#34--,-.0&>>% )&)DD,25*80,' :=??7EE033N?,&(-+%DDD(DD$*'#/' ::  ' 1a=$71aAq SM SM< N# 7K  :  2 .&$ :   [JJ1 % :: DB# XXNNHHPP9933PP ^^\\)>QQYYcc# . BFFNN  J/6 << / ++5 +"# 2X1b53q!t8#bAg-1b582=2X1b52ae8#bBh.ArE9Bq"uHDrIq"uHr!R%x"QT'!Bq!tG+bd2R74% -% -0 CC ) Y9:7 DD KK K) 8 " P $MM KK  LL @4++:X<<jj66--66 WgXX*Yff*W[:, 4f _^^( ' EEHHXBB , :KK I9 ; 00 0.4Y@;;LKKDD VVEE<jjG  FFXX 3 55cc XX :QQ ??SS/ ;5IH3Q= 4qqH L]UU VVYY W B" <3~ P** 66;; EE OO 66 // D Hl P ; II9 A DD 0 H>> d>> 0  !%4LL&78' P))LL%K 6 ,  , II>( &"& $ .-- -- +CCM M77$'':  NNGN :CI?00 / /**>>*JJ1+ 00<  $ $..HH..,,* - 355AA88 : 8 J A 7 **##  ##4 370 C >,C <<XXLL 7 7 & @ 00CEE2G <<:FF (( 33  P 9 GC;; MM N  3 3 N;# MMAAA YY==88==>> MM MM QQ ) ==HH:: > =GG . : CRR ??% 233 )9 9) ). +N F 0 (/GU2 A A IIRKK 23 B4 B@  Y Y KK6 5 5 I I.$>>@&"(B ? ?77+ G 23477;; ;; K   < <* ; ;r