*h6SrSSKJr SSKJr SSKJr SSKJrJ r J r SSK J r SSK JrJrJrJr SSKJr SS KJr SS KJr SS KJrJr SS KJr SS KJrJrJ r SSK!J"r"J#r#J$r$ SSK%J&r&J'r' SSK(J)r)J*r* SSK+J,r,J-r-J.r. SSK/J0r0J1r1 SSK2J3r3J4r4J5r5 SSK6J7r7J8r8 SBSjr9"SS\5r:"SS\5r;"SS\5r<"SS\5r="SS\5r>"S S!\5r?"S"S#\5r@"S$S%\5rA"S&S'\5rBS(rC"S)S*\5rD"S+S,\5rE"S-S.\5rF"S/S0\F5rG"S1S2\F5rH"S3S4\F5rI"S5S6\F5rJ"S7S8\5rK"S9S:\K5rL"S;S<\K5rM"S=S>\5rN"S?S@\5rOgA)Cz|This module contains various functions that are special cases of incomplete gamma functions. It should probably be renamed. ) EulerGamma)Add)cacheit)DefinedFunctionArgumentIndexError expand_mul)fuzzy_or)IpiRationalInteger)is_eq)Pow)S)Dummyuniquely_named_symbol)sympify) factorial factorial2RisingFactorial) polar_liftre unpolarify)ceilingfloor)sqrtroot)explog exp_polar)coshsinh)cossinsinc)hypermeijergc jURSR(aAU(a(SUS'UR"U40UD6[R4$U[R4$U(a1URSR"U40UD6R 5up4OURSR 5up4UR U[U--5UR U[U-- 5-S- nUR U[U--5UR U[U-- 5- S[-- nXV4$)NrFcomplex)argsis_extended_realexpandrZero as_real_imagfuncr )selfdeephintsxyrims d/home/mande/repo/quber/.venv/lib/python3.13/site-packages/sympy/functions/special/error_functions.pyreal_to_real_as_real_imagr8s yy|$$ $E) KK..7 7!&&> ! yy|""4151>>@1yy|((* ))A!G tyyQqS1 11 4B ))A!G tyyQqS1 1AaC 8B 8Oc^\rSrSrSrSrSSjrSSjr\S5r \ \ S55r Sr S rS rS rS rS rSrSrSrSrSrSrSrSSjrSrSrSrU4Sjr\r Sr!U=r"$)erf1aZ The Gauss error function. Explanation =========== This function is defined as: .. math :: \mathrm{erf}(x) = \frac{2}{\sqrt{\pi}} \int_0^x e^{-t^2} \mathrm{d}t. Examples ======== >>> from sympy import I, oo, erf >>> from sympy.abc import z Several special values are known: >>> erf(0) 0 >>> erf(oo) 1 >>> erf(-oo) -1 >>> erf(I*oo) oo*I >>> erf(-I*oo) -oo*I In general one can pull out factors of -1 and $I$ from the argument: >>> erf(-z) -erf(z) The error function obeys the mirror symmetry: >>> from sympy import conjugate >>> conjugate(erf(z)) erf(conjugate(z)) Differentiation with respect to $z$ is supported: >>> from sympy import diff >>> diff(erf(z), z) 2*exp(-z**2)/sqrt(pi) We can numerically evaluate the error function to arbitrary precision on the whole complex plane: >>> erf(4).evalf(30) 0.999999984582742099719981147840 >>> erf(-4*I).evalf(30) -1296959.73071763923152794095062*I See Also ======== erfc: Complementary error function. erfi: Imaginary error function. erf2: Two-argument error function. erfinv: Inverse error function. erfcinv: Inverse Complementary error function. erf2inv: Inverse two-argument error function. References ========== .. [1] https://en.wikipedia.org/wiki/Error_function .. [2] https://dlmf.nist.gov/7 .. [3] https://mathworld.wolfram.com/Erf.html .. [4] https://functions.wolfram.com/GammaBetaErf/Erf TcUS:Xa/S[URSS-*5-[[5- $[ X5eNr*rrr+rr rr1argindexs r7fdiff erf.fdiffs< q=S$))A,/)**483 3$T4 4r9c[$z( Returns the inverse of this function. erfinvrAs r7inverse erf.inverses  r9cfUR(aU[RLa[R$U[RLa[R$U[R La[R $UR(a[R$[U[5(aURS$[U[5(a [RURS- $UR(a[R$[U[5(a-URSR(aURS$UR[5nU[R[R 4;aU$UR!5(a U"U*5*$gNrr?) is_NumberrNaNInfinityOneNegativeInfinity NegativeOneis_zeror. isinstancerHr+erfcinverf2invextract_multiplicativelyr could_extract_minus_signclsargts r7evalerf.evals( ==aee|uu  "uu ***}}$vv c6 " "88A;  c7 # #55388A;& & ;;66M c7 # # (;(;88A;   ( ( + Q//0 0J  ' ' ) )I:  *r9cJUS:d US-S:Xa[R$[U5n[US- [S5- 5n[ U5S:aUS*US--US- -X-- $S[R U--X--U[ U5-[[5-- $Nrr*r?) rr.rrlenrRrrr nr4previous_termsks r7 taylor_termerf.taylor_terms q5AEQJ66M Aq1uadl#A>"Q&&r**QT1QU;QSAA))AD0!IaL.b2IJJr9cZURURSR55$Nrr0r+ conjugater1s r7_eval_conjugateerf._eval_conjugate"yy1//122r9c>URSRSLaggNrTr+r,rms r7 _eval_is_realerf._eval_is_reals 99Q< ( (D 0 1r9c>URSRSLaggrr)r+ is_imaginaryrms r7_eval_is_imaginaryerf._eval_is_imaginarys 99Q< $ $ , -r9cbURSn[URUR/5$rj)r+r is_finiter,r1zs r7_eval_is_finiteerf._eval_is_finites) IIaLa&8&89::r9cnURSRSLaURSR$grr)r+r,rSrms r7 _eval_is_zeroerf._eval_is_zeros1 99Q< ( (D 099Q<'' ' 1r9cnURSRSLaURSR$grr)r+r,is_extended_positiverms r7_eval_is_positiveerf._eval_is_positive1 99Q< ( (D 099Q<44 4 1r9cnURSRSLaURSR$grr)r+r,is_extended_negativerms r7_eval_is_negativeerf._eval_is_negativerr9c SSKJn [US-5U- [RU"[R US-5[[ 5- - -$Nr uppergammar*'sympy.functions.special.gamma_functionsrrrrPHalfr r1r}kwargsrs r7_eval_rewrite_as_uppergammaerf._eval_rewrite_as_uppergammas=FAqDz!|QUUZ1%=d2h%FFGGr9c [R[- U-[[5- n[R[-[ U5[[ U5-- -$NrrPr rr fresnelcfresnelsr1r}rr[s r7_eval_rewrite_as_fresnelserf._eval_rewrite_as_fresnels@uuqy!mDH$ HSMAhsmO;<>>r9c [US-5U- U[[RUS-5-[[5- - $Nr*rexpintrrr rs r7_eval_rewrite_as_expinterf._eval_rewrite_as_expints6AqDz!|aqvvq!t 44T"X===r9c "SSKJn U(aWU"X[R5nU[RLa-[R [ U*5[US-*5--$[R[ U5[US-*5-- $)Nr)limitr*) sympy.series.limitsrrrOrQrR_erfsrrP)r1r}limitvarrrlims r7_eval_rewrite_as_tractableerf._eval_rewrite_as_tractablesl- QZZ0Ca(((}}uaRyadU';;;uuuQxQTE ***r9c :[R[U5- $r)rrPerfcrs r7_eval_rewrite_as_erfcerf._eval_rewrite_as_erfcsuutAwr9c 6[*[[U-5-$rr erfirs r7_eval_rewrite_as_erfierf._eval_rewrite_as_erfisr$qs)|r9cDURSRXUS9nURUS5nU[RLaUR USUS:XaSOSS9nXR ;a&UR(aSU-[[5- $URU5$)Nrlogxcdirr-+dirr*) r+as_leading_termsubsrComplexInfinityr free_symbolsrSrr r0r1r4rrr[arg0s r7_eval_as_leading_termerf._eval_as_leading_termsiil**1d*Cxx1~ 1$$ $99Qdbjsc9BD   T\\S5b> !99T? "r9c>SSKJn USnU[R[R4;aUR SnUR U5upU *n U R(a[X- 5n [U 5V s/sH:n [RU -[SU -S- 5-USU -S--SU --- PM< sn U"SXz-- U5/-n [R[US-*5[!["5- [%U 6-- $[&[(U]WXX45$![[4a Us$f=fs sn fNrOrderr*r?)sympy.series.orderrrrOrQr+leadterm ValueErrorNotImplementedError is_positiverrangerRrrPrrr rsuperr; _eval_aseries)r1rdargs0r4rrpointr}_exnewnrfs __class__s r7rerf._eval_aseriess1,a QZZ!3!34 4 ! A  1 B~~qt}#Dk+)]]A% 1Q37(;;q1Q37|aQRd?RS)+.3AagIq.A-BCuuQTE 48 3sAw>>>S$-a@@ 34    +sD#AD;#D87D8r?r)#__name__ __module__ __qualname____firstlineno____doc__ unbranchedrCrI classmethodr] staticmethodrrgrnrtrxr~rrrrrrrrrrrrrrr8r/__static_attributes__ __classcell__rs@r7r;r;1sJXJ5B  K  K3 ;(55H==N?>+ #A*-Lr9r;c\rSrSrSrSrSSjrSSjr\S5r \ \ S55r Sr S rSS jrS rS rSrSrSrSrSrSrSrSr\rSrSrg )ri aR Complementary Error Function. Explanation =========== The function is defined as: .. math :: \mathrm{erfc}(x) = \frac{2}{\sqrt{\pi}} \int_x^\infty e^{-t^2} \mathrm{d}t Examples ======== >>> from sympy import I, oo, erfc >>> from sympy.abc import z Several special values are known: >>> erfc(0) 1 >>> erfc(oo) 0 >>> erfc(-oo) 2 >>> erfc(I*oo) -oo*I >>> erfc(-I*oo) oo*I The error function obeys the mirror symmetry: >>> from sympy import conjugate >>> conjugate(erfc(z)) erfc(conjugate(z)) Differentiation with respect to $z$ is supported: >>> from sympy import diff >>> diff(erfc(z), z) -2*exp(-z**2)/sqrt(pi) It also follows >>> erfc(-z) 2 - erfc(z) We can numerically evaluate the complementary error function to arbitrary precision on the whole complex plane: >>> erfc(4).evalf(30) 0.0000000154172579002800188521596734869 >>> erfc(4*I).evalf(30) 1.0 - 1296959.73071763923152794095062*I See Also ======== erf: Gaussian error function. erfi: Imaginary error function. erf2: Two-argument error function. erfinv: Inverse error function. erfcinv: Inverse Complementary error function. erf2inv: Inverse two-argument error function. References ========== .. [1] https://en.wikipedia.org/wiki/Error_function .. [2] https://dlmf.nist.gov/7 .. [3] https://mathworld.wolfram.com/Erfc.html .. [4] https://functions.wolfram.com/GammaBetaErf/Erfc TcUS:Xa/S[URSS-*5-[[5- $[ X5e)Nr?rarr*r@rAs r7rC erfc.fdiffos< q=c499Q<?*++DH4 4$T4 4r9c[$rFrUrAs r7rI erfc.inverseus r9cUR(agU[RLa[R$U[RLa[R$UR (a[R $[U[5(a [R URS- $[U[5(aURS$UR (a[R $UR[5nU[R[R4;aU*$UR5(a SU"U*5- $gNrr*)rMrrNrOr.rSrPrTrHr+rUrWr rQrXrYs r7r] erfc.eval|s ==aee|uu  "vv uu c6 " "55388A;& & c7 # #88A;  ;;55L  ( ( + Q//0 04K  ' ' ) )sC4y=  *r9cvUS:Xa[R$US:d US-S:Xa[R$[U5n[ US- [S5- 5n[ U5S:aUS*US--US- -X-- $S[R U--X--U[U5-[[5-- $r`) rrPr.rrrbrRrrr rcs r7rgerfc.taylor_terms 655L Ua!eqj66M Aq1uadl#A>"Q&&r**QT1QU;QSAA!--**QT11Yq\>$r(3JKKr9cZURURSR55$rjrkrms r7rnerfc._eval_conjugaterpr9cxURSRSLagURSRSLagg)NrTF)r+r,rwrms r7rterfc._eval_is_reals9 99Q< ( (D 0 99Q< $ $ , -r9Nc JUR[5RSSUS9$N tractableT)r2rrewriter;r1r}rrs r7rerfc._eval_rewrite_as_tractable#||C ((4((SSr9c :[R[U5- $r)rrPr;rs r7_eval_rewrite_as_erferfc._eval_rewrite_as_erfsuus1v~r9c V[R[[[U-5--$r)rrPr rrs r7rerfc._eval_rewrite_as_erfisuuqac{""r9c [R[- U-[[5- n[R[R[-[ U5[[ U5-- -- $rrrs r7rerfc._eval_rewrite_as_fresnelssIuuqy!mDH$uu HSMAhsmO$CDDDr9c [R[- U-[[5- n[R[R[-[ U5[[ U5-- -- $rrrs r7rerfc._eval_rewrite_as_fresnelcsIuuQwk$r("uu HSMAhsmO$CDDDr9c [RU[[5- [ [R //S/[ SS5/US-5-- $r)rrPrr r'rr rs r7rerfc._eval_rewrite_as_meijergsCuuqbz'166(Bhr1o=NPQSTPT"UUUUr9c [RSU-[[5- [ [R /S[R -/US-*5-- $r)rrPrr r&rrs r7rerfc._eval_rewrite_as_hypersAuuqs48|E166(QqvvXJA$FFFFr9c SSKJn [R[ US-5U- [RU"[R US-5[ [ 5- - -- $r)rrrrPrrr rs r7r erfc._eval_rewrite_as_uppergammasFFuutAqDz!|QUUZ1-Ed2h-N%NOOOr9c [R[US-5U- - U[[RUS-5-[[ 5- -$r)rrPrrrr rs r7rerfc._eval_rewrite_as_expints?uutAqDz!|#aqvvq!t(<&>> from sympy import I, oo, erfi >>> from sympy.abc import z Several special values are known: >>> erfi(0) 0 >>> erfi(oo) oo >>> erfi(-oo) -oo >>> erfi(I*oo) I >>> erfi(-I*oo) -I In general one can pull out factors of -1 and $I$ from the argument: >>> erfi(-z) -erfi(z) >>> from sympy import conjugate >>> conjugate(erfi(z)) erfi(conjugate(z)) Differentiation with respect to $z$ is supported: >>> from sympy import diff >>> diff(erfi(z), z) 2*exp(z**2)/sqrt(pi) We can numerically evaluate the imaginary error function to arbitrary precision on the whole complex plane: >>> erfi(2).evalf(30) 18.5648024145755525987042919132 >>> erfi(-2*I).evalf(30) -0.995322265018952734162069256367*I See Also ======== erf: Gaussian error function. erfc: Complementary error function. erf2: Two-argument error function. erfinv: Inverse error function. erfcinv: Inverse Complementary error function. erf2inv: Inverse two-argument error function. References ========== .. [1] https://en.wikipedia.org/wiki/Error_function .. [2] https://mathworld.wolfram.com/Erfi.html .. [3] https://functions.wolfram.com/GammaBetaErf/Erfi TcUS:Xa.S[URSS-5-[[5- $[ X5er>r@rAs r7rC erfi.fdiff.s9 q=S1q))$r(2 2$T4 4r9c:UR(agU[RLa[R$UR(a[R$U[R La[R $UR(a[R$UR 5(a U"U*5*$UR[5nUbU[R La[$[U[5(a[URS-$[U[5(a'[[RURS- -$[U[5(a5URSR(a[URS-$gggrL)rMrrNrSr.rOrXrWr rTrHr+rUrPrVrZr}nzs r7r] erfi.eval4s ;;AEEzuu vv ajjzz! 9966M % % ' 'G8O ' ' * >QZZ"f%%|#"g&&!%%"''!*,--"g&&2771:+=+=|#,>& r9c US:d US-S:Xa[R$[U5n[US- [S5- 5n[ U5S:aUSUS--US- -X-- $SX--U[ U5-[ [5-- $r`)rr.rrrbrrr rcs r7rgerfi.taylor_termRs q5AEQJ66M Aq1uadl#A>"Q&%b)AqD0AE:AC@@14x9Q<R!899r9cZURURSR55$rjrkrms r7rnerfi._eval_conjugate_rpr9c4URSR$rjrsrms r7_eval_is_extended_realerfi._eval_is_extended_realbyy|,,,r9c4URSR$rjr+rSrms r7rerfi._eval_is_zeroeyy|###r9c JUR[5RSSUS9$rr r s r7rerfi._eval_rewrite_as_tractablehrr9c 6[*[[U-5-$r)r r;rs r7rerfi._eval_rewrite_as_erfksr#ac({r9c B[[[U-5-[- $r)r rrs r7rerfi._eval_rewrite_as_erfcnsac{Qr9c [R[-U-[[5- n[R[- [ U5[[ U5-- -$rrrs r7rerfi._eval_rewrite_as_fresnelsqrr9c [R[-U-[[5- n[R[- [ U5[[ U5-- -$rrrs r7rerfi._eval_rewrite_as_fresnelcurr9c U[[5- [[R//S/[ SS5/US-*5-$rrrs r7rerfi._eval_rewrite_as_meijergys9bz'166(Bhr1o5FANNNr9c SU-[[5- [[R/S[R-/US-5-$rrrs r7rerfi._eval_rewrite_as_hyper|s6s48|E166(QqvvXJ1===r9c SSKJn [US-*5U- U"[RUS-*5[[ 5- [R - -$r)rrrrrr rPrs r7r erfi._eval_rewrite_as_uppergammasAFQTE{1}j!Q$7R@155HIIr9c [US-*5U- U[[RUS-*5-[[5- - $rrrs r7rerfi._eval_rewrite_as_expints:QTE{1}qA!66tBx???r9c ,UR[5$rr r!s r7r"erfi._eval_expand_funcr$r9c0URSRXUS9nURUS5nXR;a&UR(aSU-[ [ 5- $UR(aURU5$URU5$)Nrrr*) r+rrrrSrr r{r0rs r7rerfi._eval_as_leading_termsuiil**1d*Cxx1~   T\\S5b> ! ^^99T? "yy~r9c>SSKJn USnU[RLaURSn[ U5Vs/sH&n[ SU-S- 5SU-USU-S---- PM( snU"SXq-- U5/-n [*[US-5[[5- [U 6--$[[U];XX45$s snfr)rrrrOr+rrr rrr rrrr r1rdrr4rrrr}rfrrs r7rerfi._eval_aseriess,a AJJ  ! A"1X'%AaC!G$1q1Q37|(;<%'*/!$*:);J@!-L B Br9rc|\rSrSrSrSr\S5rSrSr Sr Sr S r S r S rS rS rSrSrSrSrSrg)erf2ia Two-argument error function. Explanation =========== This function is defined as: .. math :: \mathrm{erf2}(x, y) = \frac{2}{\sqrt{\pi}} \int_x^y e^{-t^2} \mathrm{d}t Examples ======== >>> from sympy import oo, erf2 >>> from sympy.abc import x, y Several special values are known: >>> erf2(0, 0) 0 >>> erf2(x, x) 0 >>> erf2(x, oo) 1 - erf(x) >>> erf2(x, -oo) -erf(x) - 1 >>> erf2(oo, y) erf(y) - 1 >>> erf2(-oo, y) erf(y) + 1 In general one can pull out factors of -1: >>> erf2(-x, -y) -erf2(x, y) The error function obeys the mirror symmetry: >>> from sympy import conjugate >>> conjugate(erf2(x, y)) erf2(conjugate(x), conjugate(y)) Differentiation with respect to $x$, $y$ is supported: >>> from sympy import diff >>> diff(erf2(x, y), x) -2*exp(-x**2)/sqrt(pi) >>> diff(erf2(x, y), y) 2*exp(-y**2)/sqrt(pi) See Also ======== erf: Gaussian error function. erfc: Complementary error function. erfi: Imaginary error function. erfinv: Inverse error function. erfcinv: Inverse Complementary error function. erf2inv: Inverse two-argument error function. References ========== .. [1] https://functions.wolfram.com/GammaBetaErf/Erf2/ cURup#US:Xa"S[US-*5-[[5- $US:Xa"S[US-*5-[[5- $[ X5e)Nr?rar*)r+rrr rr1rBr4r5s r7rC erf2.fdiffs`yy q=c1a4%j=b) ) ]S!Q$Z<R( ($T4 4r9c`[R[R[R4nU[RLdU[RLa[R$X:Xa[R$X;dX#;a[ U5[ U5- $[ U[5(a"URSU:XaURS$UR(dUUR(dDUR(aUR(d"UR(a(UR(a[ U5[ U5- $UR5nUR5nU(aU(a U"U*U*5*$U(dU(a[ U5[ U5- $grL) rrOrQr.rNr;rTrVr+rSr, is_infiniterX)rZr4r5chksign_xsign_ys r7r] erf2.evalszz1--qvv6 :aee55L V66M Xq6CF? " a ! !affQi1n66!9  99 Q%7%7AMM""q}}q6CF? "++-++- vQBK< q6#a&= r9cURURSR5URSR55$rLrkrms r7rnerf2._eval_conjugate s5yy1//1499Q<3I3I3KLLr9ctURSR=(a URSR$rLrsrms r7r5erf2._eval_is_extended_real s)yy|,,N11N1NNr9c 0[U5[U5- $rr;r1r4r5rs r7rerf2._eval_rewrite_as_erfs1vAr9c 0[U5[U5- $rrrfs r7rerf2._eval_rewrite_as_erfcsAwa  r9c Z[[[U-5[[U-5- -$rrrfs r7rerf2._eval_rewrite_as_erfis"$qs)D1I%&&r9c |[U5R[5[U5R[5- $r)r;r rrfs r7rerf2._eval_rewrite_as_fresnels'1v~~h'#a&..*BBBr9c |[U5R[5[U5R[5- $r)r;r rrfs r7rerf2._eval_rewrite_as_fresnelcror9c |[U5R[5[U5R[5- $r)r;r r'rfs r7rerf2._eval_rewrite_as_meijergs'1v~~g&Q)@@@r9c |[U5R[5[U5R[5- $r)r;r r&rfs r7rerf2._eval_rewrite_as_hyper"s'1v~~e$s1v~~e'<<tBx&GGH AJqL!%%*QVVQT":48"CC DE Fr9c |[U5R[5[U5R[5- $r)r;r rrfs r7rerf2._eval_rewrite_as_expint*s'1v~~f%Av(>>>r9c ,UR[5$rr r!s r7r"erf2._eval_expand_func-r$r9c&[UR6$r)rr+rms r7rerf2._eval_is_zero0sdii  r9rN)rrrrrrCrr]rnr5rrrrrrrrrr"rrrr9r7rVrVsiBJ5!!0MO!'CCA=F ?!!r9rVcH\rSrSrSrS SjrS Sjr\S5rSr Sr Sr g ) rHi3a Inverse Error Function. The erfinv function is defined as: .. math :: \mathrm{erf}(x) = y \quad \Rightarrow \quad \mathrm{erfinv}(y) = x Examples ======== >>> from sympy import erfinv >>> from sympy.abc import x Several special values are known: >>> erfinv(0) 0 >>> erfinv(1) oo Differentiation with respect to $x$ is supported: >>> from sympy import diff >>> diff(erfinv(x), x) sqrt(pi)*exp(erfinv(x)**2)/2 We can numerically evaluate the inverse error function to arbitrary precision on [-1, 1]: >>> erfinv(0.2).evalf(30) 0.179143454621291692285822705344 See Also ======== erf: Gaussian error function. erfc: Complementary error function. erfi: Imaginary error function. erf2: Two-argument error function. erfcinv: Inverse Complementary error function. erf2inv: Inverse two-argument error function. References ========== .. [1] https://en.wikipedia.org/wiki/Error_function#Inverse_functions .. [2] https://functions.wolfram.com/GammaBetaErf/InverseErf/ cUS:XaK[[5[URURS5S-5-[ R -$[X5eNr?rr*rr rr0r+rrrrAs r7rC erfinv.fdifffsG q=8C $))A, 7 :;;AFFB B$T4 4r9c[$rFrerAs r7rIerfinv.inversels  r9cU[RLa[R$U[RLa[R$UR(a[R $U[R La[R$[U[5(a-URSR(aURS$UR(a[R $URS5nUbE[U[5(a/URSR(aURS*$gggNrr) rrNrRrQrSr.rPrOrTr;r+r,rWr-s r7r] erfinv.evalss :55L !-- %% % YY66M !%%Z::  a  !&&)"<"<66!9  9966M ' ' + >z"c22 7T7TGGAJ; 8U2>r9c [SU- 5$Nr?rrs r7_eval_rewrite_as_erfcinverfinv._eval_rewrite_as_erfcinvsqs|r9c4URSR$rjr9rms r7rerfinv._eval_is_zeror;r9rNr) rrrrrrCrIrr]rrrrr9r7rHrH3s0/d5 *$r9rHcN\rSrSrSrS SjrS Sjr\S5rSr Sr Sr S r g ) rUia@ Inverse Complementary Error Function. The erfcinv function is defined as: .. math :: \mathrm{erfc}(x) = y \quad \Rightarrow \quad \mathrm{erfcinv}(y) = x Examples ======== >>> from sympy import erfcinv >>> from sympy.abc import x Several special values are known: >>> erfcinv(1) 0 >>> erfcinv(0) oo Differentiation with respect to $x$ is supported: >>> from sympy import diff >>> diff(erfcinv(x), x) -sqrt(pi)*exp(erfcinv(x)**2)/2 See Also ======== erf: Gaussian error function. erfc: Complementary error function. erfi: Imaginary error function. erf2: Two-argument error function. erfinv: Inverse error function. erf2inv: Inverse two-argument error function. References ========== .. [1] https://en.wikipedia.org/wiki/Error_function#Inverse_functions .. [2] https://functions.wolfram.com/GammaBetaErf/InverseErfc/ cUS:XaL[[5*[URURS5S-5-[ R -$[X5errrAs r7rC erfcinv.fdiffsI q=H9S499Q>> from sympy import erf2inv, oo >>> from sympy.abc import x, y Several special values are known: >>> erf2inv(0, 0) 0 >>> erf2inv(1, 0) 1 >>> erf2inv(0, 1) oo >>> erf2inv(0, y) erfinv(y) >>> erf2inv(oo, y) erfcinv(-y) Differentiation with respect to $x$ and $y$ is supported: >>> from sympy import diff >>> diff(erf2inv(x, y), x) exp(-x**2 + erf2inv(x, y)**2) >>> diff(erf2inv(x, y), y) sqrt(pi)*exp(erf2inv(x, y)**2)/2 See Also ======== erf: Gaussian error function. erfc: Complementary error function. erfi: Imaginary error function. erf2: Two-argument error function. erfinv: Inverse error function. erfcinv: Inverse complementary error function. References ========== .. [1] https://functions.wolfram.com/GammaBetaErf/InverseErf2/ cURup#US:Xa#[URX#5S-US-- 5$US:Xa>[[5[ R -[URX#5S-5-$[X5eNr?r*)r+rr0rr rrrrXs r7rC erf2inv.fdiffspyy q=tyy~q(A-. . ]8AFF?3tyy~q'8#99 9$T4 4r9cU[RLdU[RLa[R$UR(a!UR(a[R$UR(a#U[RLa[R $U[RLa!UR(a[R$UR(a [ U5$U[R La [U*5$UR(aU$U[R La [ U5$UR(a,UR(a[R$[ U5$UR(aU$gr)rrNrSr.rPrOrHrU)rZr4r5s r7r] erf2inv.eval s :aee55L YY19966M YY1:::  !%%ZAII55L YY!9  !**_A2;  YYH !**_!9  99yyvv ay 99H r9chURupUR(aUR(aggg)NTr9)r1r4r5s r7rerf2inv._eval_is_zero;s$yy 99#9r9rN) rrrrrrCrr]rrrr9r7rVrVs&0f54r9rVc^\rSrSrSr\S5rSSjrU4SjrSr Sr Sr S r \ r \ r\ rSS jrS rU4S jrSU4S jjrU4SjrSrU=r$)EiiDa  The classical exponential integral. Explanation =========== For use in SymPy, this function is defined as .. math:: \operatorname{Ei}(x) = \sum_{n=1}^\infty \frac{x^n}{n\, n!} + \log(x) + \gamma, where $\gamma$ is the Euler-Mascheroni constant. If $x$ is a polar number, this defines an analytic function on the Riemann surface of the logarithm. Otherwise this defines an analytic function in the cut plane $\mathbb{C} \setminus (-\infty, 0]$. **Background** The name exponential integral comes from the following statement: .. math:: \operatorname{Ei}(x) = \int_{-\infty}^x \frac{e^t}{t} \mathrm{d}t If the integral is interpreted as a Cauchy principal value, this statement holds for $x > 0$ and $\operatorname{Ei}(x)$ as defined above. Examples ======== >>> from sympy import Ei, polar_lift, exp_polar, I, pi >>> from sympy.abc import x >>> Ei(-1) Ei(-1) This yields a real value: >>> Ei(-1).n(chop=True) -0.219383934395520 On the other hand the analytic continuation is not real: >>> Ei(polar_lift(-1)).n(chop=True) -0.21938393439552 + 3.14159265358979*I The exponential integral has a logarithmic branch point at the origin: >>> Ei(x*exp_polar(2*I*pi)) Ei(x) + 2*I*pi Differentiation is supported: >>> Ei(x).diff(x) exp(x)/x The exponential integral is related to many other special functions. For example: >>> from sympy import expint, Shi >>> Ei(x).rewrite(expint) -expint(1, x*exp_polar(I*pi)) - I*pi >>> Ei(x).rewrite(Shi) Chi(x) + Shi(x) See Also ======== expint: Generalised exponential integral. E1: Special case of the generalised exponential integral. li: Logarithmic integral. Li: Offset logarithmic integral. Si: Sine integral. Ci: Cosine integral. Shi: Hyperbolic sine integral. Chi: Hyperbolic cosine integral. uppergamma: Upper incomplete gamma function. References ========== .. [1] https://dlmf.nist.gov/6.6 .. [2] https://en.wikipedia.org/wiki/Exponential_integral .. [3] Abramowitz & Stegun, section 5: https://web.archive.org/web/20201128173312/http://people.math.sfu.ca/~cbm/aands/page_228.htm cUR(a[R$U[RLa[R$U[RLa[R$UR(a[R$UR 5up#U(a[ U5S[-[-U--$gr) rSrrQrOr.extract_branch_factorrr r rZr}r.rds r7r]Ei.evals 99%% % !**_::  !$$ $66M 99%% %'') b6AaCF1H$ $ r9cp[URS5nUS:Xa[U5U- $[X5erL)rr+rrr1rBr[s r7rCEi.fdiffs41& q=s8C< $T4 4r9c>URS[S5- R(a,[TU]U5[ [ -R U5-$[TU]U5$r)r+rrr _eval_evalfr r )r1precrs r7rEi._eval_evalfsQ IIaLB ' 4 47&t,"/A/A$/GG Gw"4((r9c VSSKJn U"S[S5U-5*[[-- $)Nrrr)rrrr r rs r7rEi._eval_rewrite_as_uppergammas)F1jnQ.//!B$66r9c P[S[S5U-5*[[-- $)Nr?r)rrr r rs r7rEi._eval_rewrite_as_expints$q*R.*++ad22r9c [U[5(a[URS5$[[ U55$rj)rTrlir+rrs r7_eval_rewrite_as_liEi._eval_rewrite_as_lis1 a  affQi= #a&zr9c UR(a%[U5[U5-[[-- $[U5[U5-$r) is_negativeShiChir r rs r7_eval_rewrite_as_SiEi._eval_rewrite_as_Sis6 ==q6CF?QrT) )q6CF? "r9c 0[U5[U5-$r)r_eisr s r7rEi._eval_rewrite_as_tractables1vQr9c SSKJn [[SU/5R5nU"[ R U-U- U[ RU45$NrIntegralr\)sympy.integrals.integralsrrrnamerExp1rQr1r}rrr\s r7_eval_rewrite_as_IntegralEi._eval_rewrite_as_IntegralsE6 'aS166 7 ! a););Q%?@@r9c>SSKJn URSRUS5nURSR XS9nUR X5nUR (aqURU5upxUc [U5OUn[U5X--[-U"U5R(a[[-- $[R- $[T U]AXUS9$)Nr)r)rr)sympyrr+rrrrSas_coeff_exponentrrrr r rr.rr) r1r4rrrx0r[cers r7rEi._eval_as_leading_terms YYq\  1 %iil**1*8wwq ::((+DA!\3q6tDq6AF?Z/4,,": :23&&: :w,Q,EEr9c>URSRUS5nUR(a+UR"UR6nUR XU5$[ TU]XU5$rj)r+rrSr _eval_nseriesrr1r4rdrrrfrs r7rEi._eval_nseriess[ YYq\  1 % ::(($))4A??1. .w$Q400r9cT>SSKJn USnU[R[R4;a`UR Sn[ U5Vs/sHn[U5Xx-- PM snU"SXq-- U5/-n [U5U- [U 6-$[[U]3XX45$s snfNrrr?) rrrrOrQr+rrrrrrrrSs r7rEi._eval_aseriess,a QZZ!3!34 4 ! A05a911&91QT61%&'AF1HQ' 'R,Qq?? :s B%rrrr)rrrrrrr]rCrrrrr_eval_rewrite_as_Ci_eval_rewrite_as_Chi_eval_rewrite_as_Shirrrrrrrrs@r7rrDsrTn % %5) 7 3# ... A F1 @ @r9rc|^\rSrSrSr\S5rSrSrSr Sr Sr \ r \ r \ rS U4S jjrU4S jrS rS rU=r$)ria Generalized exponential integral. Explanation =========== This function is defined as .. math:: \operatorname{E}_\nu(z) = z^{\nu - 1} \Gamma(1 - \nu, z), where $\Gamma(1 - \nu, z)$ is the upper incomplete gamma function (``uppergamma``). Hence for $z$ with positive real part we have .. math:: \operatorname{E}_\nu(z) = \int_1^\infty \frac{e^{-zt}}{t^\nu} \mathrm{d}t, which explains the name. The representation as an incomplete gamma function provides an analytic continuation for $\operatorname{E}_\nu(z)$. If $\nu$ is a non-positive integer, the exponential integral is thus an unbranched function of $z$, otherwise there is a branch point at the origin. Refer to the incomplete gamma function documentation for details of the branching behavior. Examples ======== >>> from sympy import expint, S >>> from sympy.abc import nu, z Differentiation is supported. Differentiation with respect to $z$ further explains the name: for integral orders, the exponential integral is an iterated integral of the exponential function. >>> expint(nu, z).diff(z) -expint(nu - 1, z) Differentiation with respect to $\nu$ has no classical expression: >>> expint(nu, z).diff(nu) -z**(nu - 1)*meijerg(((), (1, 1)), ((0, 0, 1 - nu), ()), z) At non-postive integer orders, the exponential integral reduces to the exponential function: >>> expint(0, z) exp(-z)/z >>> expint(-1, z) exp(-z)/z + exp(-z)/z**2 At half-integers it reduces to error functions: >>> expint(S(1)/2, z) sqrt(pi)*erfc(sqrt(z))/sqrt(z) At positive integer orders it can be rewritten in terms of exponentials and ``expint(1, z)``. Use ``expand_func()`` to do this: >>> from sympy import expand_func >>> expand_func(expint(5, z)) z**4*expint(1, z)/24 + (-z**3 + z**2 - 2*z + 6)*exp(-z)/24 The generalised exponential integral is essentially equivalent to the incomplete gamma function: >>> from sympy import uppergamma >>> expint(nu, z).rewrite(uppergamma) z**(nu - 1)*uppergamma(1 - nu, z) As such it is branched at the origin: >>> from sympy import exp_polar, pi, I >>> expint(4, z*exp_polar(2*pi*I)) I*pi*z**3/3 + expint(4, z) >>> expint(nu, z*exp_polar(2*pi*I)) z**(nu - 1)*(exp(2*I*pi*nu) - 1)*gamma(1 - nu) + expint(nu, z) See Also ======== Ei: Another related function called exponential integral. E1: The classical case, returns expint(1, z). li: Logarithmic integral. Li: Offset logarithmic integral. Si: Sine integral. Ci: Cosine integral. Shi: Hyperbolic sine integral. Chi: Hyperbolic cosine integral. uppergamma References ========== .. [1] https://dlmf.nist.gov/8.19 .. [2] https://functions.wolfram.com/GammaBetaErf/ExpIntegralE/ .. [3] https://en.wikipedia.org/wiki/Exponential_integral c SSKJnJn [U5nX:wa [ XR5$UR (aUS::d%UR (d:SU-R (a&[[ X!S- -U"SU- U5-55$UR5up&U[RLagUR(a^US:dg[ X5S[-[-U-[RUS- --[US- 5- [U5US- --- $[S[-[-U-U-5S- X!S- --U"SU- 5-[ X5-$)Nr)gammarr*r?)rrrrr is_Integerrrrr. is_integerr r rRrr)rZnur}rrnu2rds r7r] expint.evalcs.On 9#> ! ==R1WR]]"?P?Pj!VZB5J)JKL L&&( ;  ==6"=B$q&(1==26229R!V3DDZPQ]UWZ[U[E\\] ]!Br ! $q(!1f+5eAFmCfRmS Sr9c URup#US:Xa X2S- -*[/SS/SSSU- //U5-$US:Xa[US- U5*$[X5er)r+r'rr)r1rBrr}s r7rC expint.fdiffysh  q=QK<QFQ1r6NB JJ J ]261%% %$T4 4r9c 6SSKJn X!S- -U"SU- U5-$)Nrrr?)rr)r1rr}rrs r7r"expint._eval_rewrite_as_uppergammas!F6{:a"fa000r9c US:Xa2[U[[*[-5-5*[[-- $UR(aUS:a[ U5*nXAS- -[ US- 5- [U5R[5-[U5[ US- 5- [[US- 5Vs/sHn[ X- S- 5XE--PM sn6--$U$s snfr) rr r r rrrE1r rrr)r1rr}rr4rfs r7_eval_rewrite_as_Eiexpint._eval_rewrite_as_Eis 7qA2b5))**QrT1 1 ]]rAvAAAv;ya00Ar1BBAya((%Q-H-Qi +AD0-HIJJ JKIs>C' c VUR[5R"[40UD6$r)r rrr!s r7r"expint._eval_expand_funcs!||B''8%88r9c @US:waU$[U5[U5- $r)rr)r1rr}rs r7rexpint._eval_rewrite_as_Sis 7K1vAr9cn>URSRU5(dURSnUS:Xa+UR"UR6nURXU5$UR(a1US:a+UR "UR6nURXU5$[ TU] XU5$rL)r+hasrrrrr)r1r4rdrrrrrs r7rexpint._eval_nseriessyy|""1BQw,,dii8qT2226,,dii8qT22w$Q400r9cz>SSKJn USnURSnU[RLauURSn[ U5V s/sH'n [R U -[Xy5-X-- PM) sn U"SX-- U5/-n [U*5U- [U 6-$[[U]3XX45$s sn fr) rrr+rrOrrRrrrrrr) r1rdrr4rrrrr}rfrrs r7rexpint._eval_aseriess,a YYq\ AJJ  ! AKPQR8T8a!OB$::QTA8TX]^_`a`d^dfgXhWiiAGAIa( (VT01CCUs .B8cSSKJn URupE[[ SU5R 5nU"Xd*-[ U*U-5-US[R45$Nrrr\r?) rrr+rrrrrrO)r1r+rrrdr4r\s r7r expint._eval_rewrite_as_IntegralsS6yy 'T277 82QBqD )Aq!**+=>>r9rr)rrrrrrr]rCrrr"rrrrrrrrrrs@r7rrsbdNTT*51 9... 1 D??r9rc[SU5$)a Classical case of the generalized exponential integral. Explanation =========== This is equivalent to ``expint(1, z)``. Examples ======== >>> from sympy import E1 >>> E1(0) expint(1, 0) >>> E1(5) expint(1, 5) See Also ======== Ei: Exponential integral. expint: Generalised exponential integral. li: Logarithmic integral. Li: Offset logarithmic integral. Si: Sine integral. Ci: Cosine integral. Shi: Hyperbolic sine integral. Chi: Hyperbolic cosine integral. r?)r)r}s r7rrs@ !Q<r9c\rSrSrSr\S5rSSjrSrSr Sr Sr S r \ r S r\rS rS rSSjrSSjrSrSrg )ria The classical logarithmic integral. Explanation =========== For use in SymPy, this function is defined as .. math:: \operatorname{li}(x) = \int_0^x \frac{1}{\log(t)} \mathrm{d}t \,. Examples ======== >>> from sympy import I, oo, li >>> from sympy.abc import z Several special values are known: >>> li(0) 0 >>> li(1) -oo >>> li(oo) oo Differentiation with respect to $z$ is supported: >>> from sympy import diff >>> diff(li(z), z) 1/log(z) Defining the ``li`` function via an integral: >>> from sympy import integrate >>> integrate(li(z)) z*li(z) - Ei(2*log(z)) >>> integrate(li(z),z) z*li(z) - Ei(2*log(z)) The logarithmic integral can also be defined in terms of ``Ei``: >>> from sympy import Ei >>> li(z).rewrite(Ei) Ei(log(z)) >>> diff(li(z).rewrite(Ei), z) 1/log(z) We can numerically evaluate the logarithmic integral to arbitrary precision on the whole complex plane (except the singular points): >>> li(2).evalf(30) 1.04516378011749278484458888919 >>> li(2*I).evalf(30) 1.0652795784357498247001125598 + 3.08346052231061726610939702133*I We can even compute Soldner's constant by the help of mpmath: >>> from mpmath import findroot >>> findroot(li, 2) 1.45136923488338 Further transformations include rewriting ``li`` in terms of the trigonometric integrals ``Si``, ``Ci``, ``Shi`` and ``Chi``: >>> from sympy import Si, Ci, Shi, Chi >>> li(z).rewrite(Si) -log(I*log(z)) - log(1/log(z))/2 + log(log(z))/2 + Ci(I*log(z)) + Shi(log(z)) >>> li(z).rewrite(Ci) -log(I*log(z)) - log(1/log(z))/2 + log(log(z))/2 + Ci(I*log(z)) + Shi(log(z)) >>> li(z).rewrite(Shi) -log(1/log(z))/2 + log(log(z))/2 + Chi(log(z)) - Shi(log(z)) >>> li(z).rewrite(Chi) -log(1/log(z))/2 + log(log(z))/2 + Chi(log(z)) - Shi(log(z)) See Also ======== Li: Offset logarithmic integral. Ei: Exponential integral. expint: Generalised exponential integral. E1: Special case of the generalised exponential integral. Si: Sine integral. Ci: Cosine integral. Shi: Hyperbolic sine integral. Chi: Hyperbolic cosine integral. References ========== .. [1] https://en.wikipedia.org/wiki/Logarithmic_integral .. [2] https://mathworld.wolfram.com/LogarithmicIntegral.html .. [3] https://dlmf.nist.gov/6 .. [4] https://mathworld.wolfram.com/SoldnersConstant.html cUR(a[R$U[RLa[R$U[R La[R $UR(a[R$gr)rSrr.rPrQrOrs r7r]li.evalBsS 9966M !%%Z%% % !**_::  9966M r9czURSnUS:Xa[R[U5- $[ X5erLr+rrPrrrs r7rCli.fdiffM4iil q=553s8# #$T4 4r9cURSnUR(dURUR55$grj)r+rr0rlr|s r7rnli._eval_conjugateTs2 IIaL%%99Q[[]+ +&r9c 0[U5[S5-$r)Lirrs r7_eval_rewrite_as_Lili._eval_rewrite_as_LiZ!ur!u}r9c *[[U55$r)rrrs r7rli._eval_rewrite_as_Ei]s#a&zr9c SSKJn U"S[U5*5*[R[[U55[[R [U5- 5- --[[U5*5- $Nrr)rrrrrrPrs r7rli._eval_rewrite_as_uppergamma`s`FAAw''CF c!%%A,&7789;>Aw<H Ir9c N[[[U5-5[[[[U5-5-- [R [[R [U5- 5[[U55- -- [[[U5-5- $r)Cir rSirrrPrs r7rli._eval_rewrite_as_Siesp1SV8 qAc!fH~-AEE#a&L)CAK789;>qQx=I Jr9c [[U55[[U55- [R[[R [U5- 5[[U55- -- $r)rrrrrrPrs r7rli._eval_rewrite_as_ShiksJCF c#a&k)AFFCc!f 4ECPQF 4S,TTUr9c [U5[SS[U55-[R[[U55[[R[U5- 5- --[ -$)N)r?r?)r*r*)rr&rrrPrrs r7rli._eval_rewrite_as_hyperpsYAuVVSV44CF c!%%A,&7789;EF Gr9c [[U5*5*[R[[R[U5- 5[[U55- -- [ SS[U5*5- $)N)rr))rrr)rrrrPr'rs r7rli._eval_rewrite_as_meijergtsYc!fW AEE#a&L(9CAK(G HH*lSVG<= >r9Nc 0U[[U55-$r)rrr s r7rli._eval_rewrite_as_tractablexs4A<r9cURSn[SU5Vs/sH n[U5U-[U5U-- PM" nn[[[U55-[ U6-$s snfrL)r+rrrrr)r1r4rdrrr}rfrs r7rli._eval_nseries{sa IIaL7 2 r9rcR\rSrSrSr\S5rS SjrSrSr S Sjr S S jr S r g)ria The offset logarithmic integral. Explanation =========== For use in SymPy, this function is defined as .. math:: \operatorname{Li}(x) = \operatorname{li}(x) - \operatorname{li}(2) Examples ======== >>> from sympy import Li >>> from sympy.abc import z The following special value is known: >>> Li(2) 0 Differentiation with respect to $z$ is supported: >>> from sympy import diff >>> diff(Li(z), z) 1/log(z) The shifted logarithmic integral can be written in terms of $li(z)$: >>> from sympy import li >>> Li(z).rewrite(li) li(z) - li(2) We can numerically evaluate the logarithmic integral to arbitrary precision on the whole complex plane (except the singular points): >>> Li(2).evalf(30) 0 >>> Li(4).evalf(30) 1.92242131492155809316615998938 See Also ======== li: Logarithmic integral. Ei: Exponential integral. expint: Generalised exponential integral. E1: Special case of the generalised exponential integral. Si: Sine integral. Ci: Cosine integral. Shi: Hyperbolic sine integral. Chi: Hyperbolic cosine integral. References ========== .. [1] https://en.wikipedia.org/wiki/Logarithmic_integral .. [2] https://mathworld.wolfram.com/LogarithmicIntegral.html .. [3] https://dlmf.nist.gov/6 cU[RLa[R$U[S5:Xa[R$gr)rrOr.rs r7r]Li.evals0  ?::  !A$Y66Mr9czURSnUS:Xa[R[U5- $[ X5erLrrs r7rCLi.fdiffrr9cJUR[5RU5$r)r revalf)r1rs r7rLi._eval_evalfs||B%%d++r9c 0[U5[S5- $r)rrs r7rLi._eval_rewrite_as_lir r9Nc HUR[5RSSS9$)Nr T)r2)r rr s r7rLi._eval_rewrite_as_tractables!||B'' $'??r9cXUR"UR6nURXU5$r)rr+r)r1r4rdrrrs r7rLi._eval_nseriess'  $ $dii 0qT**r9rrrr) rrrrrrr]rCrrrrrrr9r7rrs6=@ 5,@+r9rcV^\rSrSrSr\S5rS SjrSrSr S U4Sjjr Sr U=r $) TrigonometricIntegraliz(Base class for trigonometric integrals. ctU[RLa UR$U[RLaUR 5$U[R LaUR 5$UR(a UR$UR[[55nUc*URS5S:XaUR[5nUbURUS5$UR[[*55nUbURUS5$UR[S55nUc&URS5S:XaURS5nUbURU5$UR5up#US:XaX!:XagS[-[-U-URS5-U"U5-$)Nrr?rr*)rr._atzerorO_atinfrQ _atneginfrSrWrr _trigfunc_Ifactor _minusfactorrr rs r7r]TrigonometricIntegral.evalsj ;;;  !**_::<  !$$ $==? " 99;;   ' ' 1 6 :#--*a/++A.B ><<A& &  ' ' A2 7 ><<B' '  ' ' 2 7 :#--*a/++B/B >##B' ''') 6bg tAvax a((3r722r9c|[URS5nUS:XaURU5U- $[X5erL)rr+r6rrs r7rCTrigonometricIntegral.fdiff s:1& q=>>#&s* *$T4 4r9c JURU5R[5$r)rr rrs r7r)TrigonometricIntegral._eval_rewrite_as_Eis++A.66r::r9c NSSKJn URU5RU5$r)rrrr rs r7r1TrigonometricIntegral._eval_rewrite_as_uppergammas!F++A.66zBBr9c>URSRUS5S:wa[TU] XU5$UR U5RXU5nUR S5S:waUS-nUR [ SSS9nUR S5S:waU[[U5-- nURXRS5RXU5$)Nrr?cX-U- $rr)r\rds r75TrigonometricIntegral._eval_nseries..s !$q&r9F) simultaneous) r+rrrr6replacerrr)r1r4rdrr baseseriesrs r7r#TrigonometricIntegral._eval_nseriess 99Q<  Q "a '7(t4 4^^A&44Q4@ >>!  ! !OJ''-@u'U >>!  ! *s1v- -Jq))A,/==aDIIr9rrr) rrrrrrr]rCrrrrrrs@r7r1r1s6333>5;C J Jr9r1c^\rSrSrSr\r\Rr \ S5r \ S5r \ S5r \ S5rSrSr\rS rU4S jrS rS rU=r$) ri$a< Sine integral. Explanation =========== This function is defined by .. math:: \operatorname{Si}(z) = \int_0^z \frac{\sin{t}}{t} \mathrm{d}t. It is an entire function. Examples ======== >>> from sympy import Si >>> from sympy.abc import z The sine integral is an antiderivative of $sin(z)/z$: >>> Si(z).diff(z) sin(z)/z It is unbranched: >>> from sympy import exp_polar, I, pi >>> Si(z*exp_polar(2*I*pi)) Si(z) Sine integral behaves much like ordinary sine under multiplication by ``I``: >>> Si(I*z) I*Shi(z) >>> Si(-z) -Si(z) It can also be expressed in terms of exponential integrals, but beware that the latter is branched: >>> from sympy import expint >>> Si(z).rewrite(expint) -I*(-expint(1, z*exp_polar(-I*pi/2))/2 + expint(1, z*exp_polar(I*pi/2))/2) + pi/2 It can be rewritten in the form of sinc function (by definition): >>> from sympy import sinc >>> Si(z).rewrite(sinc) Integral(sinc(_t), (_t, 0, z)) See Also ======== Ci: Cosine integral. Shi: Hyperbolic sine integral. Chi: Hyperbolic cosine integral. Ei: Exponential integral. expint: Generalised exponential integral. sinc: unnormalized sinc function E1: Special case of the generalised exponential integral. li: Logarithmic integral. Li: Offset logarithmic integral. References ========== .. [1] https://en.wikipedia.org/wiki/Trigonometric_integral c0[[R-$rr rrrZs r7r4 Si._atinfns!&&yr9c2[*[R-$rrJrKs r7r5 Si._atneginfrss166zr9c[U5*$r)rrs r7r8Si._minusfactorvs 1v r9c,[[U5-U-$r)r rrZr}signs r7r7 Si._IfactorzsQx}r9c [S- [[[5U-5[[[*5U-5- S- [- -$r)r rrr rs r7rSi._eval_rewrite_as_expint~s=!tr*Q-/*R A2q0@-AA1DQFFFr9c xSSKJn [[SU/5R5nU"[ U5USU45$r)rrrrrr%rs r7rSi._eval_rewrite_as_Integrals66 'aS166 7Q!Q++r9cNURSRXUS9nURUS5nU[RLa-UR US[ U5R(aSOSS9nUR(aU$UR(dURU5$U$Nrrrrr r+rrrrNrrrrSr[r0rs r7rSi._eval_as_leading_termsiil**1d*Cxx1~ 155=99Qbh.B.Bs9LD <<J!!99T? "Kr9cx>SSKJn USnU[RLaURSn[ US-S-5Vs/sH1n[R U-[SU-5-USU-S--- PM3 snU"SXq-- U5/-n [ US-5Vs/sH4n[R U-[SU-S-5-USUS---- PM6 snU"SXq-- U5/-n [S- [U5[U 6-- [U5[U 6-- $[[U];XX45$s snfs snfr)rrrrOr+rrRrr r#rr$rrr) r1rdrr4rrrr}rfpqrs r7rSi._eval_aseriessJ,a AJJ  ! A"1a4!8_.,!IacN2Q1q\A,.16qvq1A0BCA#1a4[*(!IacAg$66QAYG(*-21QT61-=,>?Aa4#a&a.(3q6#q'>9 9R,Qq??.*s 8D2;D7cFURSnUR(aggrrr9r|s r7rSi._eval_is_zeror!r9r)rrrrrr$r6rr.r3rr4r5r8r7rr_eval_rewrite_as_sincrrrrrrs@r7rr$sDLIffGG, 7 @ r9rc^\rSrSrSr\r\Rr \ S5r \ S5r \ S5r \ S5rSrSrS rU4S jrS rU=r$) ria  Cosine integral. Explanation =========== This function is defined for positive $x$ by .. math:: \operatorname{Ci}(x) = \gamma + \log{x} + \int_0^x \frac{\cos{t} - 1}{t} \mathrm{d}t = -\int_x^\infty \frac{\cos{t}}{t} \mathrm{d}t, where $\gamma$ is the Euler-Mascheroni constant. We have .. math:: \operatorname{Ci}(z) = -\frac{\operatorname{E}_1\left(e^{i\pi/2} z\right) + \operatorname{E}_1\left(e^{-i \pi/2} z\right)}{2} which holds for all polar $z$ and thus provides an analytic continuation to the Riemann surface of the logarithm. The formula also holds as stated for $z \in \mathbb{C}$ with $\Re(z) > 0$. By lifting to the principal branch, we obtain an analytic function on the cut complex plane. Examples ======== >>> from sympy import Ci >>> from sympy.abc import z The cosine integral is a primitive of $\cos(z)/z$: >>> Ci(z).diff(z) cos(z)/z It has a logarithmic branch point at the origin: >>> from sympy import exp_polar, I, pi >>> Ci(z*exp_polar(2*I*pi)) Ci(z) + 2*I*pi The cosine integral behaves somewhat like ordinary $\cos$ under multiplication by $i$: >>> from sympy import polar_lift >>> Ci(polar_lift(I)*z) Chi(z) + I*pi/2 >>> Ci(polar_lift(-1)*z) Ci(z) + I*pi It can also be expressed in terms of exponential integrals: >>> from sympy import expint >>> Ci(z).rewrite(expint) -expint(1, z*exp_polar(-I*pi/2))/2 - expint(1, z*exp_polar(I*pi/2))/2 See Also ======== Si: Sine integral. Shi: Hyperbolic sine integral. Chi: Hyperbolic cosine integral. Ei: Exponential integral. expint: Generalised exponential integral. E1: Special case of the generalised exponential integral. li: Logarithmic integral. Li: Offset logarithmic integral. References ========== .. [1] https://en.wikipedia.org/wiki/Trigonometric_integral c"[R$r)rr.rKs r7r4 Ci._atinfs vv r9c[[-$r)r r rKs r7r5 Ci._atneginfs t r9c4[U5[[--$rrr r rs r7r8Ci._minusfactors!uqt|r9c@[U5[[-S- U--$rrr r rRs r7r7 Ci._Ifactor s1v"Qt ##r9c z[[[5U-5[[[*5U-5-*S- $r)rrr rs r7rCi._eval_rewrite_as_expints2JqM!O$r*aR.*:';;>r9c SSKJn [[SU/5R5n[ R [U5-U"S[U5- U- USU45- $r) rrrrrrrrr#rs r7rCi._eval_rewrite_as_IntegralsP6 'aS166 7||c!f$x3q61 q!Qi'HHHr9cURSRXUS9nURUS5nU[RLa-UR US[ U5R(aSOSS9nUR(a:URU5upgUc [U5OUn[U5Xr--[-$UR(aURU5$U$rZr+rrrrNrrrrSrrrr{r0r1r4rrr[rrrs r7rCi._eval_as_leading_termiil**1d*Cxx1~ 155=99Qbh.B.Bs9LD <<((+DA!\3q6tDq6AF?Z/ / ^^99T? "Kr9c>SSKJn USnU[R[R4;GaUR Sn[ US-S-5Vs/sH1n[RU-[SU-5-USU-S--- PM3 snU"SXq-- U5/-n [ US-5Vs/sH4n[RU-[SU-S-5-USUS---- PM6 snU"SXq-- U5/-n [U5[U 6-[U5[U 6-- n U[RLaU [[-- n U $[[U]CXX45$s snfs snfr)rrrrOrQr+rrRrr$rr#r r rrr) r1rdrr4rrrr}rfr^r_resultrs r7rCi._eval_aseries&si,a QZZ!3!34 4 ! A"1a4!8_.,!IacN2Q1q\A,.16qvq1A0BCA#1a4[*(!IacAg$66QAYG(*-21QT61-=,>?AVS!W%AQ(88F***!B$MR,Qq??.*s 8E/;E$r)rrrrrr#r6rrr3rr4r5r8r7rrrrrrrs@r7rrsM^IG$$?I @@r9rc\rSrSrSr\r\Rr \ S5r \ S5r \ S5r \ S5rSrSrS rS rg ) ri8ai Sinh integral. Explanation =========== This function is defined by .. math:: \operatorname{Shi}(z) = \int_0^z \frac{\sinh{t}}{t} \mathrm{d}t. It is an entire function. Examples ======== >>> from sympy import Shi >>> from sympy.abc import z The Sinh integral is a primitive of $\sinh(z)/z$: >>> Shi(z).diff(z) sinh(z)/z It is unbranched: >>> from sympy import exp_polar, I, pi >>> Shi(z*exp_polar(2*I*pi)) Shi(z) The $\sinh$ integral behaves much like ordinary $\sinh$ under multiplication by $i$: >>> Shi(I*z) I*Si(z) >>> Shi(-z) -Shi(z) It can also be expressed in terms of exponential integrals, but beware that the latter is branched: >>> from sympy import expint >>> Shi(z).rewrite(expint) expint(1, z)/2 - expint(1, z*exp_polar(I*pi))/2 - I*pi/2 See Also ======== Si: Sine integral. Ci: Cosine integral. Chi: Hyperbolic cosine integral. Ei: Exponential integral. expint: Generalised exponential integral. E1: Special case of the generalised exponential integral. li: Logarithmic integral. Li: Offset logarithmic integral. References ========== .. [1] https://en.wikipedia.org/wiki/Trigonometric_integral c"[R$rrrOrKs r7r4 Shi._atinf{ zzr9c"[R$r)rrQrKs r7r5 Shi._atneginfs!!!r9c[U5*$r)rrs r7r8Shi._minusfactors Awr9c,[[U5-U-$r)r rrRs r7r7 Shi._IfactorsAwt|r9c [U5[[[[-5U-5- S- [[-S- - $r)rr r r rs r7rShi._eval_rewrite_as_expints519QrT?1,--q01R4699r9cFURSnUR(aggrrr9r|s r7rShi._eval_is_zeror!r9cPURSRU5nURUS5nU[RLa-UR US[ U5R(aSOSS9nUR(aU$UR(dURU5$U$)Nrrrrr[rs r7rShi._eval_as_leading_terms~iil**1-xx1~ 155=99Qbh.B.Bs9LD <<J!!99T? "Kr9rN)rrrrrr"r6rr.r3rr4r5r8r7rrrrrr9r7rr8su=~IffG"": r9rc\rSrSrSr\r\Rr \ S5r \ S5r \ S5r \ S5rSrSrS rg ) riaX Cosh integral. Explanation =========== This function is defined for positive $x$ by .. math:: \operatorname{Chi}(x) = \gamma + \log{x} + \int_0^x \frac{\cosh{t} - 1}{t} \mathrm{d}t, where $\gamma$ is the Euler-Mascheroni constant. We have .. math:: \operatorname{Chi}(z) = \operatorname{Ci}\left(e^{i \pi/2}z\right) - i\frac{\pi}{2}, which holds for all polar $z$ and thus provides an analytic continuation to the Riemann surface of the logarithm. By lifting to the principal branch we obtain an analytic function on the cut complex plane. Examples ======== >>> from sympy import Chi >>> from sympy.abc import z The $\cosh$ integral is a primitive of $\cosh(z)/z$: >>> Chi(z).diff(z) cosh(z)/z It has a logarithmic branch point at the origin: >>> from sympy import exp_polar, I, pi >>> Chi(z*exp_polar(2*I*pi)) Chi(z) + 2*I*pi The $\cosh$ integral behaves somewhat like ordinary $\cosh$ under multiplication by $i$: >>> from sympy import polar_lift >>> Chi(polar_lift(I)*z) Ci(z) + I*pi/2 >>> Chi(polar_lift(-1)*z) Chi(z) + I*pi It can also be expressed in terms of exponential integrals: >>> from sympy import expint >>> Chi(z).rewrite(expint) -expint(1, z)/2 - expint(1, z*exp_polar(I*pi))/2 - I*pi/2 See Also ======== Si: Sine integral. Ci: Cosine integral. Shi: Hyperbolic sine integral. Ei: Exponential integral. expint: Generalised exponential integral. E1: Special case of the generalised exponential integral. li: Logarithmic integral. Li: Offset logarithmic integral. References ========== .. [1] https://en.wikipedia.org/wiki/Trigonometric_integral c"[R$rr}rKs r7r4 Chi._atinfrr9c"[R$rr}rKs r7r5 Chi._atneginfrr9c4[U5[[--$rrmrs r7r8Chi._minusfactors1v"}r9c@[U5[[-S- U--$rrjrRs r7r7 Chi._Ifactors!uqtAvd{""r9c [*[-S- [U5[[[[-5U-5-S- - $r)r r rr rs r7rChi._eval_rewrite_as_expint s7r"uQw"Q%"Yqt_Q%6"77:::r9cURSRXUS9nURUS5nU[RLa-UR US[ U5R(aSOSS9nUR(a:URU5upgUc [U5OUn[U5Xr--[-$UR(aURU5$U$rZrtrus r7rChi._eval_as_leading_term rwr9rN)rrrrrr!r6rrr3rr4r5r8r7rrrrr9r7rrssHTIG##; r9rcP\rSrSrSrSr\S5rS SjrSr \ r Sr Sr \ rS rg ) FresnelIntegrali z%Base class for the Fresnel integrals.TcpU[RLa[R$UR(a[R$[R nUnSnUR S5nUbU*nUnSnUR [5nUbUR[-U-nUnSnU(a X "U5-$g)NFrT) rrOrrSr.rPrWr _sign)rZr}prefactnewargchangedr.s r7r]FresnelIntegral.eval s  ?66M 9966M%%  , ,R 0 >hGFG  , ,Q / >iik')GFG 3v;& & r9cUS:Xa9UR[R[-URSS--5$[ X5er)r6rrr r+rrAs r7rCFresnelIntegral.fdiff; s< q=>>!&&)DIIaL!O";< <$T4 4r9c4URSR$rjrsrms r7r5&FresnelIntegral._eval_is_extended_realA r7r9c4URSR$rjr9rms r7rFresnelIntegral._eval_is_zeroF r;r9cZURURSR55$rjrkrms r7rnFresnelIntegral._eval_conjugateI rpr9rNr)rrrrrrrr]rCr5r~rrnr8r/rrr9r7rr s>0J''<5 --O$3-Lr9rc^\rSrSrSr\r\R*r \ \ S55r Sr SrSrSrSrU4S jrS rU=r$) riO a Fresnel integral S. Explanation =========== This function is defined by .. math:: \operatorname{S}(z) = \int_0^z \sin{\frac{\pi}{2} t^2} \mathrm{d}t. It is an entire function. Examples ======== >>> from sympy import I, oo, fresnels >>> from sympy.abc import z Several special values are known: >>> fresnels(0) 0 >>> fresnels(oo) 1/2 >>> fresnels(-oo) -1/2 >>> fresnels(I*oo) -I/2 >>> fresnels(-I*oo) I/2 In general one can pull out factors of -1 and $i$ from the argument: >>> fresnels(-z) -fresnels(z) >>> fresnels(I*z) -I*fresnels(z) The Fresnel S integral obeys the mirror symmetry $\overline{S(z)} = S(\bar{z})$: >>> from sympy import conjugate >>> conjugate(fresnels(z)) fresnels(conjugate(z)) Differentiation with respect to $z$ is supported: >>> from sympy import diff >>> diff(fresnels(z), z) sin(pi*z**2/2) Defining the Fresnel functions via an integral: >>> from sympy import integrate, pi, sin, expand_func >>> integrate(sin(pi*z**2/2), z) 3*fresnels(z)*gamma(3/4)/(4*gamma(7/4)) >>> expand_func(integrate(sin(pi*z**2/2), z)) fresnels(z) We can numerically evaluate the Fresnel integral to arbitrary precision on the whole complex plane: >>> fresnels(2).evalf(30) 0.343415678363698242195300815958 >>> fresnels(-2*I).evalf(30) 0.343415678363698242195300815958*I See Also ======== fresnelc: Fresnel cosine integral. References ========== .. [1] https://en.wikipedia.org/wiki/Fresnel_integral .. [2] https://dlmf.nist.gov/7 .. [3] https://mathworld.wolfram.com/FresnelIntegrals.html .. [4] https://functions.wolfram.com/GammaBetaErf/FresnelS .. [5] The converging factors for the fresnel integrals by John W. Wrench Jr. and Vicki Alley crUS:a[R$[U5n[U5S:a9USn[S-*US--SU-S- -SU-SU-S--SU-S--- U-$US-US-*U--[S5SU-S- -[SU-S----SU-S-[ SU-S-5-- $) Nrr?rr*rrarr.rrbr rrdr4rer^s r7rgfresnels.taylor_term s q566M A>"Q&"2&Qq!t QqS1W-qsAaC!G}acAg/FG1LL!t1uqj(AaD2a4!8,SSKJn USnU[R[R*4;GaURSn[ SU5Vs/sHfnSU-S-U:dM[R U-[SU-S-5-SSU-S--USU-S---SSU---[SU-5-- PMh n nSSU-- /[ SU5Vs/sHlnSU-S-U:dM[R U-[SU-S- 5-SSU-S--USU-S---SSU-S- --[SU-S- 5-- PMn sn-n U V s/sHn [S[- 5*U -PM n n U V s/sHn [S[- 5*U -PM n n U[RLaSOSn U [R-[US-5[U 6-[US-5[U 6--RU[S[- 5U-5-U"SXq-- U5-$[T U]AXX45$s snfs snfs sn fs sn fNrrrrr?r*r)rrrrOr+rrRrrr rr$rr#rrrr1rdrr4rrrr}rfr^r_r\rrs r7rfresnels._eval_aseries s`,a QZZ!**- - ! A  1+6%Q1q1I!IacAg$66acAg,QqS1W-AaC81Q3GI% 6AaC 1+6%Q1q1QQ]]A- !A#'0BBacAg,QqS1W-AaC!G >> from sympy import I, oo, fresnelc >>> from sympy.abc import z Several special values are known: >>> fresnelc(0) 0 >>> fresnelc(oo) 1/2 >>> fresnelc(-oo) -1/2 >>> fresnelc(I*oo) I/2 >>> fresnelc(-I*oo) -I/2 In general one can pull out factors of -1 and $i$ from the argument: >>> fresnelc(-z) -fresnelc(z) >>> fresnelc(I*z) I*fresnelc(z) The Fresnel C integral obeys the mirror symmetry $\overline{C(z)} = C(\bar{z})$: >>> from sympy import conjugate >>> conjugate(fresnelc(z)) fresnelc(conjugate(z)) Differentiation with respect to $z$ is supported: >>> from sympy import diff >>> diff(fresnelc(z), z) cos(pi*z**2/2) Defining the Fresnel functions via an integral: >>> from sympy import integrate, pi, cos, expand_func >>> integrate(cos(pi*z**2/2), z) fresnelc(z)*gamma(1/4)/(4*gamma(5/4)) >>> expand_func(integrate(cos(pi*z**2/2), z)) fresnelc(z) We can numerically evaluate the Fresnel integral to arbitrary precision on the whole complex plane: >>> fresnelc(2).evalf(30) 0.488253406075340754500223503357 >>> fresnelc(-2*I).evalf(30) -0.488253406075340754500223503357*I See Also ======== fresnels: Fresnel sine integral. References ========== .. [1] https://en.wikipedia.org/wiki/Fresnel_integral .. [2] https://dlmf.nist.gov/7 .. [3] https://mathworld.wolfram.com/FresnelIntegrals.html .. [4] https://functions.wolfram.com/GammaBetaErf/FresnelC .. [5] The converging factors for the fresnel integrals by John W. Wrench Jr. and Vicki Alley cXUS:a[R$[U5n[U5S:a9USn[S-*US--SU-S- -SU-SU-S- -SU-S--- U-$XS-*U--[S5SU--[SU----SU-S-[ SU-5-- $) Nrr?rr*rrrrarrs r7rgfresnelc.taylor_termH s q566M A>"Q&"2&Qq!t QqS1W-qsAaC!G}acAg/FG1LLTEA:~11b1Q3i)?@QqS1WiXYZ[X[nD\]]r9c [R[- S- [[R[-S- [ [ 5-U-5[[[R[- S- [ [ 5-U-5---$rrrs r7rfresnelc._eval_rewrite_as_erfU rr9c U[[SS5/[R[SS5/[S-*US--S- 5-$)Nr?rr*r)r&r rrr rs r7rfresnelc._eval_rewrite_as_hyperX sA5(1a.)AFFHQN+Cb!eVAqD[QS^TTTr9c [U[SS5--[S5[US-S5-[U*S5-- [ /S/[SS5/[SS5S/[S-*US--S- 5-$)Nrrr*r?rr)r r rrr'rs r7r!fresnelc._eval_rewrite_as_meijerg[ s1hq!n$$QQT1 (=dA2qk(IJ"qcHQN#3hq!na5H2q5&QRTUQU+VX.YZ [r9c SSKJn [[SU/5R5nU"[ [ US--S- 5USU45$r)rrrrrr#r rs r7r"fresnelc._eval_rewrite_as_Integral_ rr9cSSKJn URSRXUS9nUR US5nU[ R La-URUS[U5R(aSOSS9nUR(aU$U[ R[ R4;a3U[ RLaSOSnU[ R-U"X5-$URU5$) Nrrrrrrr?r)rrr+rrrrrrrrSrOrQrr0rs r7rfresnelc._eval_as_leading_termd s,iil**1d*Cxx1~ 1$$ $99Qbh.B.Bs9LD <<J ajj!"4"45 5QZZ'RAQVV8eAk) )99T? "r9c >SSKJn USnU[R[R*4;GaURSn[ U5Vs/sHfnSU-S-U:dM[R U-[SU-S-5-SSU-S--USU-S---SSU---[SU-5-- PMh n nSSU-- /[ SU5Vs/sHlnSU-S-U:dM[R U-[SU-S- 5-SSU-S--USU-S---SSU-S- --[SU-S- 5-- PMn sn-n U V s/sHn [S[- 5*U -PM n n U V s/sHn [S[- 5U -PM n n U[RLaSOSn U [R-[US-5[U 6-[US-5[U 6--RU[S[- 5U-5-U"SXq-- U5-$[T U]AXX45$s snfs snfs sn fs sn fr)rrrrOr+rrRrrr rr#rr$rrrrs r7rfresnelc._eval_aseriess s\,a QZZ!**- - ! A  (3"QacAgkI!IacAg$66acAg,QqS1W-AaC81Q3GI" 3AaC 1+6%Q1q1QQ]]A- !A#'0BBacAg,QqS1W-AaC!G SSKJn USnU[RLaURSn[ U5Vs/sHPnS[ [5- [SU-5-[S5*U*--[U5- SU- SU-S---PMR n nU"SUSU-S--- U5n [U 6RX1U5U -$UR[5n U [RLaURSn[ U5Vs/sHPnS[ [5- [SU-5-[S5*U*--[U5- SU- SU-S---PMR n nU"SUSU-S--- U5n [U 6RX1U5U -$[T U]9XX45$s snfs snf)Nrrr?r*r)rrrrOr+rrr rrrrWr rr) r1rdrr4rrrr}rflor\rs r7r_erfs._eval_aseries s,a AJJ  ! ADI!HNDLq48i!n,q0/r.$Q<(+,Q3!A#'*:;DL NaAaC!G na(AG**16: :  * *1 -  ? ! AEJ!HNDLq48i!n,q0/r.$Q<(+,Q3!A#'*:;DL NaAaC!G na(AG**16: :w$Qq77%NNs AF" AF'cUS:Xa3URSnS[[5- SU-[U5--$[ X5e)Nr?rrar*)r+rr rrr1rBr}s r7rC _erfs.fdiff s@ q= ! Ad2h;1U1X- -$T4 4r9c X[R[U5- [US-5-$r)rrPr;rrs r7_eval_rewrite_as_intractable"_erfs._eval_rewrite_as_intractable s!AAqD ))r9rr) rrrrrrr]rrCrrrrs@r7rr s0 865**r9rcX^\rSrSrSrU4SjrS SjrSrU4SjrS U4Sjjr Sr U=r $) ri zr Helper function to make the $\mathrm{Ei}(z)$ and $\mathrm{li}(z)$ functions tractable for the Gruntz algorithm. c^>SSKJn US[R[R4;a[ T U]XX45$URSn[U5Vs/sHn[U5SU- US---PM nnU"SXaS--- U5n [U6RX1U5U -$s snfr) rrrrOrQrrr+rrrr) r1rdrr4rrr}rfrrrs r7r_eis._eval_aseries s, 8AJJ(:(:; ;7(1; ; IIaL49!H =HqYq\QqSAEN *H = !AAJ, "Q&&qT2Q66 >s!B*cUS:Xa.URSn[RU- [U5- $[ X5e)Nr?r)r+rrPrrrs r7rC _eis.fdiff s8 q= ! A5519tAw& &$T4 4r9c 2[U*5[U5-$r)rrrs r7r!_eis._eval_rewrite_as_intractable sA2wr!u}r9c>URSRUS5nUR(a)UR"UR6nUR XUS9$[ TU]XUS9$)Nrr)r+rrSrrr)r1r4rrrrrs r7r_eis._eval_as_leading_term sc YYq\  1 % ::11499=A**1d*C Cw,Q,EEr9c>URSRUS5nUR(a+UR"UR6nUR XU5$[ TU]XU5$rj)r+rrSrrrrs r7r_eis._eval_nseries s[ YYq\  1 % ::11499=A??1. .w$Q400r9rrr) rrrrrrrCrrrrrrs@r7rr s' 75F11r9rN)T)Pr sympy.corersympy.core.addrsympy.core.cachersympy.core.functionrrrsympy.core.logicr sympy.core.numbersr r r r sympy.core.relationalrsympy.core.powerrsympy.core.singletonrsympy.core.symbolrrsympy.core.sympifyr(sympy.functions.combinatorial.factorialsrrr$sympy.functions.elementary.complexesrrr#sympy.functions.elementary.integersrr(sympy.functions.elementary.miscellaneousrr&sympy.functions.elementary.exponentialrrr %sympy.functions.elementary.hyperbolicr!r"(sympy.functions.elementary.trigonometricr#r$r%sympy.functions.special.hyperr&r'r8r;rrrVrHrUrVrrrrrr1rrrrrrrrrrr9r7rsF"$OO%77' ":&[[LL>?FF<CC8 *l-/l-^H?HD}B?}B@O!?O!bZ$_Z$zP;P;f[o[Bt@t@n}?_}?@ FeeNZ+Z+@<JO<J~E EPJ@ J@Xg gTn nj5-o5-p^8^8B^8^8L.*O.*b*1?*1r9