Otj6dZddlmZddlmZddlmZddlmZm Z m Z ddl m Z ddl mZmZmZmZddlmZdd lmZdd lmZdd lmZmZdd lmZdd lmZmZm Z ddl!m"Z"m#Z#m$Z$ddl%m&Z&m'Z'ddl(m)Z)m*Z*ddl+m,Z,m-Z-m.Z.ddl/m0Z0m1Z1ddl2m3Z3m4Z4m5Z5ddl6m7Z7m8Z8dCdZ9GddeZ:GddeZ;GddeZ<GddeZ=Gdd eZ>Gd!d"eZ?Gd#d$eZ@Gd%d&eZAGd'd(eZBd)ZCGd*d+eZDGd,d-eZEGd.d/eZFGd0d1eFZGGd2d3eFZHGd4d5eFZIGd6d7eFZJGd8d9eZKGd:d;eKZLGd<d=eKZMGd>d?eZNGd@dAeZOdBS)Dz 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)hypermeijergTc R|jdjr/|rd|d<|j|fi|tjfS|tjfS|r/|jdj|fi|\}}n"|jd\}}||t|zz||t|zz zdz }||t|zz||t|zz z dtzz }||fS)NrFcomplex)argsis_extended_realexpandrZero as_real_imagfuncr )selfdeephintsxyrims e/var/www/html/CA-Chatbot/venv/lib/python3.11/site-packages/sympy/functions/special/error_functions.pyreal_to_real_as_real_imagr8s" y|$"  "$E) DK....7 7!&> ! +"ty|"411511>>@@11y|((**1 ))A!G  tyyQqS11 11 4B ))A!G  tyyQqS11 1AaC 8B 8OceZdZdZdZddZddZedZe e dZ dZ d Z d Zd Zd Zd ZdZdZdZdZdZdZdZddZdZdZdZfdZeZxZ S)erfa. 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 Tc|dkr6dt|jddz zttz St ||Nr<r*rrr+rr rr1argindexs r7fdiffz erf.fdiffsF q==S$)A,/)***4883 3$T844 4r9ctSz8 Returns the inverse of this function. erfinvr@s r7inversez erf.inverses  r9c|jra|tjur tjS|tjur tjS|tjur tjS|jr tjSt|tr |j dSt|trtj|j dz S|jr tjSt|tr|j djr |j dS|t}|tjtjfvr|S|r ||  SdSNrr<) is_NumberrNaNInfinityOneNegativeInfinity NegativeOneis_zeror. isinstancerFr+erfcinverf2invextract_multiplicativelyr could_extract_minus_signclsargts r7evalzerf.evals> = ae||u  ""u ***}$ v c6 " " 8A;  c7 # # '538A;& & ; 6M c7 # #  (; 8A;   ( ( + + Q/0 0 0J  ' ' ) ) CII:   r9cv|dks |dzdkr tjSt|}t|dz tdz }t |dkr|d |dzz|dz z||zz Sdtj|zz||zz|t |zttzz SNrr*r<) rr.rrlenrOrrr nr4previous_termsks r7 taylor_termzerf.taylor_terms q55AEQJJ6M Aq1uaddl##A>""Q&&&r**QT1QU;QqSAA))AqD0!IaLL.b2IJJr9cf||jdSNrr0r+ conjugater1s r7_eval_conjugatezerf._eval_conjugate&yy1//11222r9c2|jdjdurdSdSNrTr+r,rhs r7 _eval_is_realzerf._eval_is_reals# 9Q< (D 0 04 1 0r9c2|jdjdurdSdSrl)r+ is_imaginaryrhs r7_eval_is_imaginaryzerf._eval_is_imaginarys# 9Q< $ , ,4 - ,r9cR|jd}t|j|jgSre)r+r is_finiter,r1zs r7_eval_is_finitezerf._eval_is_finites% IaLa&89:::r9cR|jdjdur|jdjSdSrl)r+r,rPrhs r7 _eval_is_zerozerf._eval_is_zeros- 9Q< (D 0 09Q<' ' 1 0r9cR|jdjdur|jdjSdSrl)r+r,is_extended_positiverhs r7_eval_is_positivezerf._eval_is_positive- 9Q< (D 0 09Q<4 4 1 0r9cR|jdjdur|jdjSdSrl)r+r,is_extended_negativerhs r7_eval_is_negativezerf._eval_is_negativer|r9c ddlm}t|dz|z tj|tj|dztt z z zSNr uppergammar*'sympy.functions.special.gamma_functionsrrrrMHalfr r1rukwargsrs r7_eval_rewrite_as_uppergammazerf._eval_rewrite_as_uppergammasQFFFFFFAqDzz!|QUZZ1%=%=d2hh%FFGGr9c tjtz |zttz }tjtzt |tt |zz zSNrrMr rr fresnelcfresnelsr1rurrXs r7_eval_rewrite_as_fresnelszerf._eval_rewrite_as_fresnelsBuqy!mDHH$ HSMMAhsmmO;<>>>r9c t|dz|z |ttj|dzzttz z SNr*rexpintrrr rs r7_eval_rewrite_as_expintzerf._eval_rewrite_as_expints;AqDzz!|aqvq!t 4 44T"XX===r9Nc ddlm}|rV|||tj}|tjur1tjt | t|dz zzStjt |t|dz zz S)Nr)limitr*) sympy.series.limitsrrrLrNrO_erfsrrM)r1rulimitvarrrlims r7_eval_rewrite_as_tractablezerf._eval_rewrite_as_tractables------  <%8QZ00Ca(((}uaRyyadU';;;uuQxxQTE ***r9c :tjt|z Sr)rrMerfcrs r7_eval_rewrite_as_erfczerf._eval_rewrite_as_erfcsutAwwr9c Bt tt|zzSrr erfirs r7_eval_rewrite_as_erfizerf._eval_rewrite_as_erfisr$qs))|r9cN|jd|||}||d}|tjur ||d|dkrdnd}||jvr!|jrd|zttz S| |S)Nrlogxcdirr-+dirr*) r+as_leading_termsubsrComplexInfinityr free_symbolsrPrr r0r1r4rrrXarg0s r7_eval_as_leading_termzerf._eval_as_leading_termsil**14d*CCxx1~~ 1$ $ $99Qdbjjssc9BBD  T\ S5b> !99T?? "r9c@ ddlm}|d}|tjtjfvr|jd |\}}n#ttf$r|cYSwxYw| }|j rt||z } fdt| D|d | zz |gz} tj t dz ttz t!| zz St#t$|||||S)NrOrderc~g|]9}tj|ztd|zdz zd|zdzzd|zzz :Sr*r<)rrOr.0rbrus r7 z%erf._eval_aseries..s`+++]A% 1Q37(;(;;q1Q37|aQRd?RS+++r9r<r*)sympy.series.orderrrrLrNr+leadterm ValueErrorNotImplementedError is_positiverrangerMrrr rsuperr; _eval_aseries) r1r`args0r4rrpoint_exnewnsru __class__s @r7rzerf._eval_aseriessD,,,,,,a QZ!34 4 4 ! A  1 22 34     B~ ?qt}}++++#Dkk+++.3eAagIq.A.A-BCuQTE 488 3sAw>>>S$--a4@@@sAA('A(r<r)!__name__ __module__ __qualname____doc__ unbranchedrBrG classmethodrZ staticmethodrrcrirnrqrvrxr{rrrrrrrrrrrrr8r/ __classcell__rs@r7r;r;1sJJXJ5555[B  K K W\ K333 ;;;(((555555HHH======NNN???>>>++++ # # #AAAAA*-LLLLLr9r;ceZdZdZdZddZddZedZe e dZ dZ d Z dd Zd Zd ZdZdZdZdZdZdZdZdZeZdZd S)ra& 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 Tr<c|dkr6dt|jddz zttz St ||)Nr<r]rr*r?r@s r7rBz erfc.fdiffosF q==c49Q<?*+++DHH4 4$T844 4r9ctSrDrRr@s r7rGz erfc.inverseus r9c|jrG|tjur tjS|tjur tjS|jr tjSt|trtj|j dz St|tr |j dS|jr tjS| t}|tjtj fvr| S|rd|| z SdSNrr*)rJrrKrLr.rPrMrQrFr+rRrTr rNrUrVs r7rZz erfc.eval|s = ae||u  ""v  u c6 " " '538A;& & c7 # # 8A;  ; 5L  ( ( + + Q/0 0 04K  ' ' ) ) !ssC4yy=  ! !r9c|dkr tjS|dks |dzdkr tjSt|}t |dz tdz }t |dkr|d |dzz|dz z||zz Sdtj|zz||zz|t|zttzz Sr\) rrMr.rrr^rOrrr r_s r7rczerfc.taylor_terms 665L UUa!eqjj6M Aq1uaddl##A>""Q&&&r**QT1QU;QqSAA!-**QT11Yq\\>$r((3JKKr9cf||jdSrerfrhs r7rizerfc._eval_conjugaterjr9c^|jdjdurdS|jdjdurdSdS)NrTF)r+r,rprhs r7rnzerfc._eval_is_reals< 9Q< (D 0 04 9Q< $ , ,5 - ,r9Nc b|tdd|SN tractableT)r2rrewriter;r1rurrs r7rzerfc._eval_rewrite_as_tractable)||C  ((4((SSSr9c :tjt|z Sr)rrMr;rs r7_eval_rewrite_as_erfzerfc._eval_rewrite_as_erfsus1vv~r9c Ztjttt|zzzSr)rrMr rrs r7rzerfc._eval_rewrite_as_erfisuqac{""r9c tjtz |zttz }tjtjtzt |tt |zz zz Srrrs r7rzerfc._eval_rewrite_as_fresnelssIuqy!mDHH$u HSMMAhsmmO$CDDDr9c tjtz |zttz }tjtjtzt |tt |zz zz Srrrs r7rzerfc._eval_rewrite_as_fresnelcsIuQwk$r(("u HSMMAhsmmO$CDDDr9c tj|ttz t tjggdgt ddg|dzzz Sr)rrMrr r'rr rs r7rzerfc._eval_rewrite_as_meijergsFuqbz'16(Bhr1oo=NPQSTPT"U"UUUUr9c tjd|zttz t tjgdtjzg|dz zz Sr)rrMrr r&rrs r7rzerfc._eval_rewrite_as_hypers@uqs488|E16(QqvXJA$F$FFFFr9c ddlm}tjt |dz|z tj|tj|dzt t z z zz Sr)rrrrMrrr rs r7rz erfc._eval_rewrite_as_uppergammasXFFFFFFutAqDzz!|QUZZ1-E-Ed2hh-N%NOOOr9c tjt|dz|z z |ttj|dzztt z zSr)rrMrrrr rs r7rzerfc._eval_rewrite_as_expintsButAqDzz!|#aqvq!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 Tr<c|dkr5dt|jddzzttz St ||r>r?r@s r7rBz erfi.fdiff.sC q==S1q)))$r((2 2$T844 4r9c|jrG|tjur tjS|jr tjS|tjur tjS|jr tjS|r ||  S|t}||tjurtSt|trt|j dzSt|tr"ttj |j dz zSt|tr)|j djrt|j dzSdSdSdSrI)rJrrKrPr.rLrUrTr rQrFr+rRrMrSrWrunzs r7rZz erfi.eval4s? ; "AEzzu  "v ajz! 9 6M % % ' ' CGG8O ' ' * * >QZ"f%% $|#"g&& .!%"'!*,--"g&& $271:+= $|# > $ $ $ $r9cT|dks |dzdkr tjSt|}t|dz tdz }t |dkr|d|dzz|dz z||zz Sd||zz|t |zt tzz Sr\)rr.rrr^rrr r_s r7rczerfi.taylor_termRs q55AEQJJ6M Aq1uaddl##A>""Q&&%b)AqD0AE:AaC@@1a4x9Q<<R!899r9cf||jdSrerfrhs r7rizerfi._eval_conjugate_rjr9c&|jdjSrermrhs r7_eval_is_extended_realzerfi._eval_is_extended_realby|,,r9c&|jdjSrer+rPrhs r7rxzerfi._eval_is_zeroey|##r9Nc b|tdd|Srrrs r7rzerfi._eval_rewrite_as_tractablehrr9c Bt tt|zzSr)r r;rs r7rzerfi._eval_rewrite_as_erfksr#ac(({r9c Pttt|zztz Sr)r rrs r7rzerfi._eval_rewrite_as_erfcnsac{Qr9c tjtz|zttz }tjtz t |tt |zz zSrrrs r7rzerfi._eval_rewrite_as_fresnelsqrr9c tjtz|zttz }tjtz t |tt |zz zSrrrs r7rzerfi._eval_rewrite_as_fresnelcurr9c |ttz ttjggdgt ddg|dz zSrrrs r7rzerfi._eval_rewrite_as_meijergys>bz'16(Bhr1oo5FANNNNr9c d|zttz ttjgdtjzg|dzzSrrrs r7rzerfi._eval_rewrite_as_hyper|s7s488|E16(QqvXJ1====r9c ddlm}t|dz |z |tj|dz tt z tjz zSr)rrrrrr rMrs r7rz erfi._eval_rewrite_as_uppergammasUFFFFFFQTE{{1}jj!Q$77R@15HIIr9c t|dz |z |ttj|dz zttz z Srrrs r7rzerfi._eval_rewrite_as_expints?QTE{{1}qA!6!66tBxx???r9c 6|tSrrrs r7rzerfi._eval_expand_funcrr9c*|jd|||}||d}||jvr!|jrd|zt t z S|jr||S||S)Nrrr*) r+rrrrPrr rsr0rs r7rzerfi._eval_as_leading_termsil**14d*CCxx1~~  T\ S5b> ! ^ #99T?? "yy~~r9cddlm}|d}|tjurv|jdfdt |D|d|zz |gz}t tdzttz t|zzStt| ||||S)Nrrc^g|])}td|zdz d|zd|zdzzzz *Sr)rrs r7rz&erfi._eval_aseries..sS'''AaC!G$$1q1Q37|(;<'''r9r<r*)rrrrLr+rr rrr rrrr r1r`rr4rrrrrurs @r7rzerfi._eval_aseriess,,,,,,a AJ   ! A''''"1XX'''*/%!Q$*:*:);>>JJJ@@@!!!-L B B B B B B B B Br9rc|eZdZdZdZedZdZdZdZ dZ dZ d Z d Z d Zd Zd ZdZdZdZdS)erf2a? 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/ c|j\}}|dkr+dt|dz zttz S|dkr+dt|dz zttz St ||)Nr<r]r*)r+rrr rr1rAr4r5s r7rBz erf2.fdiffsqy1 q==c1a4%jj=b) ) ]]S!Q$ZZ<R( ($T844 4r9ctjtjtjf}|tjus|tjur tjS||kr tjS||vs||vrt |t |z St |tr|jd|kr |jdS|j s#|j s|j r|j s|j r&|j rt |t |z S| }| }|r|r|| |  S|s|rt |t |z SdSrI) rrLrNr.rKr;rQrSr+rPr, is_infiniterU)rWr4r5chksign_xsign_ys r7rZz erf2.evals]z1-qv6 ::ae5L !VV6M #XXcq66CFF? " a ! ! afQi1nn6!9  9 #  #Q%7 #AM #" #'(} #q66CFF? "++--++--  !v !CQBKK<  ! !q66#a&&=  ! !r9c||jd|jdSrIrfrhs r7rizerf2._eval_conjugate s:yy1//1149Q<3I3I3K3KLLLr9cJ|jdjo|jdjSrIrmrhs r7r zerf2._eval_is_extended_real sy|,N11NNr9c @t|t|z Srr;r1r4r5rs r7rzerf2._eval_rewrite_as_erfs1vvAr9c @t|t|z Srrr*s r7rzerf2._eval_rewrite_as_erfcsAwwa  r9c pttt|ztt|zz zSrrr*s r7rzerf2._eval_rewrite_as_erfis&$qs))D1II%&&r9c t|tt|tz Sr)r;rrr*s r7rzerf2._eval_rewrite_as_fresnels11vv~~h''#a&&..*B*BBBr9c t|tt|tz Sr)r;rrr*s r7rzerf2._eval_rewrite_as_fresnelcr/r9c t|tt|tz Sr)r;rr'r*s r7rzerf2._eval_rewrite_as_meijergs11vv~~g&&Q)@)@@@r9c t|tt|tz Sr)r;rr&r*s r7rzerf2._eval_rewrite_as_hyper"s11vv~~e$$s1vv~~e'<'<<&>tBxx&GGH AJJqL!%**QVQT":":488"CC DE Fr9c t|tt|tz Sr)r;rrr*s r7rzerf2._eval_rewrite_as_expint*s11vv~~f%%Av(>(>>>r9c 6|tSrrrs r7rzerf2._eval_expand_func-rr9ct|jSr)rr+rhs r7rxzerf2._eval_is_zero0sdi  r9N)rrrrrBrrZrir rrrrrrrrrrrxrr9r7rrs BBJ555!![!0MMMOOO!!!'''CCCCCCAAA===FFF ???!!!!!!!!r9rcDeZdZdZd dZd dZedZdZdZ dS) rFaR 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/ r<c|dkrRttt||jddzzt jzSt||Nr<rr*rr rr0r+rrrr@s r7rBz erfinv.fdifffsR q==88C $)A, 7 7 :;;;AFB B$T844 4r9ctSrDr)r@s r7rGzerfinv.inversels  r9c|tjur tjS|tjur tjS|jr tjS|tjur tjSt|tr|j dj r |j dS|jr tjS| d}|5t|tr"|j dj r|j d SdSdSdSNrr) rrKrOrNrPr.rMrLrQr;r+r,rTrs r7rZz erfinv.evalss ::5L !-  % % Y 6M !%ZZ:  a   !&)"< 6!9  9 6M ' ' + + >z"c22> 7T>GAJ;  >>>>>r9c &td|z SNr<rrs r7_eval_rewrite_as_erfcinvzerfinv._eval_rewrite_as_erfcinvsqs||r9c&|jdjSrer rhs r7rxzerfinv._eval_is_zerorr9Nr) rrrrrBrGrrZr@rxrr9r7rFrF3s//d5555 [*$$$$$r9rFcJeZdZdZd dZd dZedZdZdZ dZ d S) rRa 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/ r<c|dkrStt t||jddzzt jzSt||r9r:r@s r7rBz erfcinv.fdiffsT q==HH9S49Q>> 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/ c6|j\}}|dkr,t|||dz|dzz S|dkrHttt jzt|||dzzSt||)Nr<r*)r+rr0rr rrrr s r7rBz erf2inv.fdiffsy1 q==tyy1~~q(A-.. . ]]88AF?3tyy1~~q'8#9#99 9$T844 4r9c&|tjus|tjur tjS|jr|jr tjS|jr|tjur tjS|tjur|jr tjS|jrt |S|tjurt| S|jr|S|tjurt |S|jr"|jr tjSt |S|jr|SdSr)rrKrPr.rMrLrFrR)rWr4r5s r7rZz erf2inv.eval s ::ae5L Y 19 6M Y 1:::  !%ZZAIZ5L Y !99  !*__A2;;  Y H !*__!99  9 !y !v ayy 9 H  r9c>|j\}}|jr |jrdSdSdS)NTr )r1r4r5s r7rxzerf2inv._eval_is_zero;s9y1 9  4    r9N)rrrrrBrrZrxrr9r7rSrSsX00f555[4r9rSceZdZdZedZddZfdZdZdZ dZ d Z e Z e Z e Zdd Zd Zfd Zdfd ZfdZxZS)Eia 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 c6|jr tjS|tjur tjS|tjur tjS|jr tjS|\}}|r%t |dtztz|zzSdSr) rPrrNrLr.extract_branch_factorrQr r rWrurr`s r7rZzEi.evals 9 % % !*__:  !$ $ $6M 9 &% %''))A  %b66AaCF1H$ $ % %r9r<ct|jd}|dkrt||z St||rI)rr+rrr1rArXs r7rBzEi.fdiffs>1&& q==s88C< $T844 4r9c|jdtdz jrDt|t t z|zSt|Sr=)r+rrr _eval_evalfr r )r1precrs r7rXzEi._eval_evalfsg IaLB ' 4 H77&&t,,"/A/A$/G/GG Gww""4(((r9c hddlm}|dtd|z ttzz S)Nrrr)rrrr r rs r7rzEi._eval_rewrite_as_uppergammas?FFFFFF 1jnnQ.///!B$66r9c dtdtd|z ttzz S)Nr<r)rrr r rs r7rzEi._eval_rewrite_as_expints)q*R..*+++ad22r9c t|trt|jdStt |Sre)rQrlir+rrs r7_eval_rewrite_as_lizEi._eval_rewrite_as_lis: a   !afQi== #a&&zzr9c |jr/t|t|zttzz St|t|zSr) is_negativeShiChir r rs r7_eval_rewrite_as_SizEi._eval_rewrite_as_SisA = #q66CFF?QrT) )q66CFF? "r9Nc @t|t|zSr)r_eisrs r7rzEi._eval_rewrite_as_tractables1vvQr9c ddlm}ttd|gj}|t j|z|z |t j|fSNrIntegralrY)sympy.integrals.integralsrirrnamerExp1rNr1rurrirYs r7_eval_rewrite_as_IntegralzEi._eval_rewrite_as_IntegralsW666666 'aS116 7 7x ! a);Q%?@@@r9cddlm}|jd|d}|jd||}|||}|jrt||\}}|t|n|}t|||zztz||j rttzn tjz St|||S)Nr)r)rr)sympyrr+rrrrPas_coeff_exponentrrr`r r rr.rr) r1r4rrrx0rXcers r7rzEi._eval_as_leading_terms Yq\  1 % %il**14*88wwq$ : :((++DAq!\3q666tDq66AdF?Z/4,8"!&: :ww,,QT,EEEr9rc|jd|d}|jr&|j|j}||||St |||Sre)r+rrPrc _eval_nseriesrr1r4r`rrrrfrs r7rvzEi._eval_nseriessk Yq\  1 % % : /(($)4A??1a.. .ww$$Q4000r9c\ddlm}|d}|tjtjfvrX|jdfdt |D|d|zz |gz}tz t|zStt| ||||S)Nrrc:g|]}t||zz Srrrs r7rz$Ei._eval_aseries..s(99911Q&999r9r<) rrrrLrNr+rrrrrQrrs @r7rzEi._eval_aseriess,,,,,,a QZ!34 4 4 ! A9999a999U1QT61%%&'AFF1HQ' 'R,,Qq$???r9rrr)rrrrrrZrBrXrrr^rc_eval_rewrite_as_Ci_eval_rewrite_as_Chi_eval_rewrite_as_Shirrnrrvrrrs@r7rQrQDsHTTn % %[ %5555))))) 777 333### ...    AAA F F F F F111111 @ @ @ @ @ @ @ @ @r9rQcteZdZdZedZdZdZdZdZ dZ e Z e Z e Z d fd Zfd Zd ZxZS) ra 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 ddlm}m}t|}||krt ||S|jr|dks|js> ! = MR1WWR]W"?PWjR!VZZB5J5J)JKKLL L&&((1 ;; F = T66"a==B$q&(1=26229R!V3D3DDZPQ]]UWZ[U[E\\] ]!Br ! $$q(!b1f+5eeAFmmCfRQRmmS Sr9c |j\}}|dkr%||dz z tgddgddd|z gg|zS|dkrt|dz | St||r9)r+r'rr)r1rArrus r7rBz expint.fdiffysy A q==QK<QFQ1r6NB J JJ J ]]261%%% %$T844 4r9c >ddlm}||dz z|d|z |zS)Nrrr<)rr)r1rrurrs r7rz"expint._eval_rewrite_as_uppergammas7FFFFFF26{::a"fa0000r9c dkr>t|tt tzz ttzz Sjrdkrt | dz zt dz z t|tztt dz z tfdtdz DzzS|S)Nr<cFg|]}t|z dz |zzSr*r{)rrbrr4s r7rz.expint._eval_rewrite_as_Ei..s2HHHQiQ ++AqD0HHHr9) rQr r r rrrE1rrrr)r1rrurr4s ` @r7_eval_rewrite_as_Eizexpint._eval_rewrite_as_Eis 77qA2b5)))***QrT1 1 ] rAvvAArAv;ya000Ar1B1BBAya(((HHHHH%Q--HHHIJJ JKr9c X|tjtfi|Sr)rrQrrs r7rzexpint._eval_expand_funcs)'t||B'88%888r9c P|dkr|St|t|z Sr?)rarb)r1rrurs r7rczexpint._eval_rewrite_as_Sis& 77K1vvAr9rcb|jd|sl|jd}|dkr&|j|j}||||S|jr,|dkr&|j|j}||||St |||SrI)r+hasrcrvrrr)r1r4r`rrrrxrs r7rvzexpint._eval_nseriessy|"" 31BQww,D,di8q!T222 3266,D,di8q!T222ww$$Q4000r9cd ddlm}|d}|jd|tjurZ|jd  fdt |D|d |zz |gz}t  z t|zStt| ||||S)Nrrr<c\g|](}tj|zt|z|zz )Sr)rrOr)rrbrrus r7rz(expint._eval_aseries..s8TTTa!OB$:$::QTATTTr9) rrr+rrLrrrrrr) r1r`rr4rrrrrrurs @@r7rzexpint._eval_aseriess,,,,,,a Yq\ AJ   ! ATTTTT5QR88TTTX]X]^_`acd`d^dfgXhXhWiiAGGAIa( (VT""00E1dCCCr9cddlm}|j\}}tt d|j}||| zt | |zz|dtjfSNrrhrYr<) rjrir+rrrkrrrL)r1r+rrir`r4rYs r7rnz expint._eval_rewrite_as_Integralsl666666y1 'T227 8 8xA2QBqD )Aq!*+=>>>r9r|)rrrrrrZrBrrrrcr}r~rrvrrnrrs@r7rrsddNTT[T*555111   999... 1 1 1 1 1 1 D D D D D???????r9rc"td|S)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)rus r7rrs@ !Q<<r9c~eZdZdZedZddZdZdZdZ dZ d Z e Z d Z e Zd Zd ZddZddZdZd S)r]a 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 c|jr tjS|tjur tjS|tjur tjS|jr tjSdSr)rPrr.rMrNrLrFs r7rZzli.evalBsU 9 6M !%ZZ% % !*__:  9 6M  r9r<c|jd}|dkrtjt|z St ||rIr+rrMrrrVs r7rBzli.fdiffM9il q==53s88# #$T844 4r9c||jd}|js'||SdSre)r+r~r0rgrts r7rizli._eval_conjugateTs< IaL% ,99Q[[]]++ + , ,r9c @t|tdzSr)Lir]rs r7_eval_rewrite_as_Lizli._eval_rewrite_as_LiZ!uur!uu}r9c :tt|Sr)rQrrs r7rzli._eval_rewrite_as_Ei]s#a&&zzr9c (ddlm}|dt|  tjtt|ttjt|z z zztt| z SNrr)rrrrrrMrs r7rzli._eval_rewrite_as_uppergamma`sFFFFFFAAw'''CFF c!%A,&7&7789;>Aw<<H Ir9c ttt|ztttt|zzz tjttjt|z tt|z zz ttt|zz Sr)Cir rSirrrMrs r7rczli._eval_rewrite_as_Sies1SVV8 qAc!ffH~-AE#a&&L))CAKK789;>qQx==I Jr9c tt|tt|z tjttjt|z tt|z zz Sr)rbrrarrrMrs r7rzli._eval_rewrite_as_ShiksWCFF c#a&&kk)AFCc!ff 4E4ECPQFF 4S,TTUr9c t|tddt|ztjtt|ttjt|z z zzt zS)N)r<r<)r*r*)rr&rrrMrrs r7rzli._eval_rewrite_as_hyperpseAuVVSVV444CFF c!%A,&7&7789;EF Gr9c &tt|  tjttjt|z tt|z zz t ddt| z S)N)rr))rrr)rrrrMr'rs r7rzli._eval_rewrite_as_meijergtsgc!ffW AE#a&&L(9(9CAKK(G HH*lSVVG<<= >r9Nc @|tt|zSr)rerrs r7rzli._eval_rewrite_as_tractablexs4A<<r9rc|jdfdtd|D}tttzt |zS)NrcZg|]'}t|zt||zz (Sr)rrrs r7rz$li._eval_nseries..}s3 C C C!c!ffq[IaLL1, - C C Cr9r<)r+rrrr)r1r4r`rrrrus @r7rvzli._eval_nseries{sO IaL C C C CuQ{{ C C CCAKK'#q'11r9c2|jd}|jrdSdSrlr rts r7rxzli._eval_is_zero& IaL 9 4  r9rrr|)rrrrrrZrBrirrrrcr}rr~rrrrvrxrr9r7r]r]s``F[5555,,, III JJJ.VVV0GGG>>>    2222 r9r]cLeZdZdZedZd dZdZdZd dZ d d Z dS)rab 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 cx|tjur tjS|tdkr tjSdSr)rrLr.rFs r7rZzLi.evals2  ??:  !A$$YY6MYr9r<c|jd}|dkrtjt|z St ||rIrrVs r7rBzLi.fdiffrr9c\|t|Sr)rr]evalf)r1rYs r7rXzLi._eval_evalfs"||B%%d+++r9c @t|tdz Sr)r]rs r7r^zLi._eval_rewrite_as_lirr9Nc `|tddS)NrT)r2)rr]rs r7rzLi._eval_rewrite_as_tractables'||B'' $'???r9rcN|j|j}||||Sr)r^r+rv)r1r4r`rrrxs r7rvzLi._eval_nseriess) $D $di 0q!T***r9rrr|) rrrrrrZrBrXr^rrvrr9r7rrs==@[ 5555,,,@@@@++++++r9rcNeZdZdZedZd dZdZdZd fd Z xZ S) TrigonometricIntegralz) Base class for trigonometric integrals. c|tjur|jS|tjur|S|tjur|S|jr|jS|tt}|3| ddkr|t}|| |dS|tt }|| |dS|td}|.| ddkr|d}|| |S|\}}|dkr||krdSdtztz|z| dz||zS)Nrr<rr*)rr._atzerorL_atinfrN _atneginfrPrTrr _trigfunc_Ifactor _minusfactorrSr rTs r7rZzTrigonometricIntegral.evals ;;;  !*__::<<  !$ $ $==?? " 9 ;   ' ' 1 6 6 :#--**a//++A..B ><<A&& &  ' ' A2 7 7 ><<B'' '  ' ' 2 7 7 :#--**a//++B//B >##B'' '''))A 66bAgg FtAvax a(((33r7722r9r<ct|jd}|dkr|||z St||rI)rr+rrrVs r7rBzTrigonometricIntegral.fdiff sE1&& q==>>#&&s* *$T844 4r9c \||tSr)rrrQrs r7rz)TrigonometricIntegral._eval_rewrite_as_Eis$++A..66r:::r9c ^ddlm}|||Sr)rrrrrs r7rz1TrigonometricIntegral._eval_rewrite_as_uppergammas6FFFFFF++A..66zBBBr9rc2|jd|ddkr#t|||S|||||}|ddkr|dz}|t dd}|ddkr|tt|zz }|||jd|||S)Nrr<c||z|z Srr)rYr`s r7z5TrigonometricIntegral._eval_nseries..s!Q$q&r9F) simultaneous) r+rrrvrreplacerrr)r1r4r`rr baseseriesrs r7rvz#TrigonometricIntegral._eval_nseriess 9Q<  Q " "a ' '77((At44 4^^A&&44Q4@@ >>!   ! ! !OJ''-@-@u'UU >>!   ! ! *s1vv- -Jq$)A,//==aDIIIr9rr|) rrrrrrZrBrrrvrrs@r7rrs3333[3>5555;;;CCC J J J J J J J J J Jr9rceZdZdZeZejZe dZ e dZ e dZ e dZ dZdZeZdZfd Zd ZxZS) ra 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 c*ttjzSrr rrrWs r7rz Si._atinfns!&yr9c,t tjzSrrrs r7rz Si._atneginfrss16zr9c"t| Sr)rrFs r7rzSi._minusfactorvs1v r9c6tt|z|zSr)r rarWrusigns r7rz Si._IfactorzsQx}r9c tdz ttt|zttt |zz dz tz zSr)r rrr rs r7rzSi._eval_rewrite_as_expint~sG!tr*Q--/**R A2q0@-A-AA1DQFFFr9c ddlm}ttd|gj}|t ||d|fSrg)rjrirrrkr%rms r7rnzSi._eval_rewrite_as_IntegralsO666666 'aS116 7 7xQ!Q+++r9c:|jd|||}||d}|tjur.||dt |jrdnd}|jr|S|j s| |S|SNrrrrr r+rrrrKrrr`rPr"r0rs r7rzSi._eval_as_leading_termsil**14d*CCxx1~~ 15==99Qbhh.B'Kss9LLD < J! 99T?? "Kr9c ddlm}|d}|tjur|jd fdt |dzdzD|d |zz |gz} fdt |dzD|d |zz |gz}t dz t t|zz t t|zz Stt| ||||S)Nrrclg|]0}tj|ztd|zzd|zdzzz 1SrrrOrrs r7rz$Si._eval_aseries..P...!IacNN2Q1q\A...r9r*r<crg|]3}tj|ztd|zdzzd|dzzzz 4Srrrs r7rz$Si._eval_aseries..V***!IacAg$6$66QAYG***r9) rrrrLr+rr r#rr$rrr) r1r`rr4rrrpqrurs @r7rzSi._eval_aseriess4,,,,,,a AJ   ! A...."1a4!8__...16qAvq1A1A0BCA****"1a4[[***-2U1QT61-=-=,>?Aa4#a&&a.(3q66#q'>9 9R,,Qq$???r9c2|jd}|jrdSdSrlr rts r7rxzSi._eval_is_zerorr9)rrrrr$rrr.rrrrrrrrn_eval_rewrite_as_sincrrrxrrs@r7rr$s DDLIfG[[[[GGG,,, 7   @@@@@ r9rceZdZdZeZejZe dZ e dZ e dZ e dZ dZdZdZfd ZxZS) ra 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 ctjSr)rr.rs r7rz Ci._atinfs v r9c ttzSr)r r rs r7rz Ci._atneginfs t r9c@t|ttzzSrrr r rFs r7rzCi._minusfactors!uuqt|r9cLt|ttzdz |zzSrrbr r rs r7rz Ci._Ifactor s1vv"Qt ##r9c ttt|zttt |zz dz Sr)rrr rs r7rzCi._eval_rewrite_as_expints<JqMM!O$$r*aR..*:';';;>r9c ddlm}ttd|gj}t jt|z|dt|z |z |d|fz Sr) rjrirrrkrrrr#rms r7rnzCi._eval_rewrite_as_Integralsi666666 'aS116 7 7|c!ff$xx3q661 q!Qi'H'HHHr9c|jd|||}||d}|tjur.||dt |jrdnd}|jrH| |\}}|t|n|}t|||zztzS|j r| |S|Srr+rrrrKrrr`rPrqrrrsr0r1r4rrrXrrsrts r7rzCi._eval_as_leading_termil**14d*CCxx1~~ 15==99Qbhh.B'Kss9LLD < ((++DAq!\3q666tDq66AdF?Z/ / ^ 99T?? "Kr9c> ddlm}|d}|tjtjfvr|jd fdt |dzdzD|d |zz |gz} fdt |dzD|d |zz |gz}t t|zt t|zz } |tjur| ttzz } | Stt|||||S)Nrrclg|]0}tj|ztd|zzd|zdzzz 1Srrrs r7rz$Ci._eval_aseries..,rr9r*r<crg|]3}tj|ztd|zdzzd|dzzzz 4Srrrs r7rz$Ci._eval_aseries...rr9)rrrrLrNr+rr$rr#r r rrr) r1r`rr4rrrrrresultrurs @r7rzCi._eval_aseries&sP,,,,,,a QZ!34 4 4 ! A...."1a4!8__...16qAvq1A1A0BCA****"1a4[[***-2U1QT61-=-=,>?AVVS!W%AQ(88F***!B$MR,,Qq$???r9)rrrrr#rrrrrrrrrrrnrrrrs@r7rrsMM^IG[[[$$[$???III   @@@@@@@@@r9rceZdZdZeZejZe dZ e dZ e dZ e dZ dZdZdZd S) raa 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 ctjSrrrLrs r7rz Shi._atinf{ zr9ctjSr)rrNrs r7rz Shi._atneginfs !!r9c"t| Sr)rarFs r7rzShi._minusfactorsAwr9c6tt|z|zSr)r rrs r7rz Shi._IfactorsAwt|r9c t|ttttz|zz dz ttzdz z Sr)rr r r rs r7rzShi._eval_rewrite_as_expints<19QrT??1,---q01R4699r9c2|jd}|jrdSdSrlr rts r7rxzShi._eval_is_zerorr9c4|jd|}||d}|tjur.||dt |jrdnd}|jr|S|j s| |S|S)Nrrrrrrs r7rzShi._eval_as_leading_termsil**1--xx1~~ 15==99Qbhh.B'Kss9LLD < J! 99T?? "Kr9N)rrrrr"rrr.rrrrrrrrxrrr9r7rara8s==~IfG[""["[[:::     r9raceZdZdZeZejZe dZ e dZ e dZ e dZ dZdZdS) rba  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 ctjSrrrs r7rz Chi._atinfrr9ctjSrrrs r7rz Chi._atneginfrr9c@t|ttzzSrrrFs r7rzChi._minusfactors1vv"}r9cLt|ttzdz |zzSrrrs r7rz Chi._Ifactors!uuqtAvd{""r9c t tzdz t|ttttz|zzdz z Sr)r r rr rs r7rzChi._eval_rewrite_as_expint s>r"uQw"Q%%"Yqt__Q%6"7"77:::r9c|jd|||}||d}|tjur.||dt |jrdnd}|jrH| |\}}|t|n|}t|||zztzS|j r| |S|Srrrs r7rzChi._eval_as_leading_term rr9N)rrrrr!rrrrrrrrrrrrr9r7rbrbsHHTIG[[[##[#;;;     r9rbcNeZdZdZdZedZd dZdZeZ dZ dZ e Z d S) FresnelIntegralz& Base class for the Fresnel integrals.Tc@|tjur tjS|jr tjStj}|}d}|d}|| }|}d}|t}||jtz|z}|}d}|r|||zSdS)NFrT) rrLrrPr.rMrTr _sign)rWruprefactnewargchangedrs r7rZzFresnelIntegral.eval s  ??6M 9 6M%  , ,R 0 0 >hGFG  , ,Q / / >ik')GFG  '33v;;& & ' 'r9r<c|dkr8|tjtz|jddzzSt ||r9)rrrr r+rr@s r7rBzFresnelIntegral.fdiff; sB q==>>!&)DIaL!O";<< <$T844 4r9c&|jdjSrermrhs r7r z&FresnelIntegral._eval_is_extended_realA r r9c&|jdjSrer rhs r7rxzFresnelIntegral._eval_is_zeroF rr9cf||jdSrerfrhs r7rizFresnelIntegral._eval_conjugateI rjr9Nr)rrrrrrrZrBr rvrxrir8r/rr9r7rr s00J''['<5555 ----O$$$333-LLLr9rczeZdZdZeZej Ze e dZ dZ dZ dZdZdZfdZxZS) ray 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 c|dkr tjSt|}t|dkr=|d}tdz |dzzd|zdz zd|zd|zdzzd|zdzzz |zS|dz|dz |zztdd|zdz ztd|zdzzzzd|zdzt d|zdzzz S) Nrr<rr*rr]rr.rr^r rr`r4rars r7rczfresnels.taylor_term s q556M A>""Q&&"2&Qq!t QqS1W-qsAaC!G}acAg/FG1LL!t1uqj(AaDD2a4!8,. s6661q1!IacAg$6$66acAg,QqS1W-AaC81Q3GI)4r9r<r*cg|]m}d|zdzktj|ztd|zdz zdd|zdzzd|zdzzzdd|zdz zztd|zdz zz nSrr<r*rr&s r7rz*fresnels._eval_aseries.. 6661q1]A- !A#'0B0BBacAg,QqS1W-AaC!G . (***1$qt**Q***r9cDg|]}tdtz  |zSrr+r,s r7rz*fresnels._eval_aseries.. r-r9r)rrrrLr+rrr$rr#rrr rr r1r`rr4rrrrrrrurs ` @r7rzfresnels._eval_aseries s,,,,,,a QZ!*- - - ! A666661++666AAaC 666661++6666A+****A*****Aaj((bAQV8s1a4yya03q!t99S!W3DD$q$qt**Q,''(*/%!Q$*:*:; ;ww$$Qq$777r9)rrrrr$rrrMrrrrcrrrrnrrrrs@r7rrO sSShI UFE  m m W\ maaaeee[[[333 # # #888888888r9rcxeZdZdZeZejZe e dZ dZ dZ dZdZdZfdZxZS) rau Fresnel integral C. Explanation =========== This function is defined by .. math:: \operatorname{C}(z) = \int_0^z \cos{\frac{\pi}{2} t^2} \mathrm{d}t. It is an entire function. Examples ======== >>> 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 c||dkr tjSt|}t|dkr=|d}tdz |dzzd|zdz zd|zd|zdz zd|zdzzz |zS||dz |zztdd|zztd|zzzzd|zdzt d|zzz S) Nrr<rr*rrrr]rrs r7rczfresnelc.taylor_termH s q556M A>""Q&&"2&Qq!t QqS1W-qsAaC!G}acAg/FG1LLQTEA:~11b1Q3i)?@QqS1WiXYZ[X[nnD\]]r9c 4tjtz dz ttjtzdz t t z|ztttjtz dz t t z|zzzzSrrrs r7rzfresnelc._eval_rewrite_as_erfU rr9c |ttddgtjtddgtdz |dzzdz zS)Nr<rr*r)r&r rrr rs r7rzfresnelc._eval_rewrite_as_hyperX sF5(1a..)AFHQNN+Cb!eVAqD[QS^TTTTr9c 0t|tddzztdt|dzdzt| dzz t gdgtddgtdddgtdz |dzzdz zS)Nrrr*r<rr)r r rrr'rs r7rz!fresnelc._eval_rewrite_as_meijerg[ s1hq!nn$$QQT1 (=dA2qkk(IJ"qcHQNN#3hq!nna5H2q5&QRTUQU+VX.YYZ [r9c ddlm}ttd|gj}|t t |dzzdz |d|fSr)rjrirrrkr#r rms r7rnz"fresnelc._eval_rewrite_as_Integral_ r r9cddlm}|jd|||}||d}|t jur.||dt|j rdnd}|j r|S|t j t j fvr.|t j urdnd}|t j z|||zS||S) Nrrrrrrr<r)rrr+rrrrrrr`rPrLrNrr0r"s r7rzfresnelc._eval_as_leading_termd s,,,,,,il**14d*CCxx1~~ 1$ $ $99Qbhh.B'Kss9LLD < #J aj!"45 5 5QZ''RAQV8eeAqkk) )99T?? "r9c ddlm}|d}|tjtj fvr|jd  fdt D}dd zz g fdt dDz}d|D}d|D}|tjurdnd } | tjzt dzt|zt dzt|zz |tdtz |zz|d zz |zSt|||S) Nrrcg|]g}d|zdzktj|ztd|zdzzdd|zdzzd|zdzzzdd|zzztd|zzz hSr%rr&s r7rz*fresnelc._eval_aseries..} s333acAgkk!IacAg$6$66acAg,QqS1W-AaC81Q3GI&1kkr9r<r*cg|]m}d|zdzktj|ztd|zdz zdd|zdzzd|zdzzzdd|zdz zztd|zdz zz nSr(rr&s r7rz*fresnelc._eval_aseries.. r)r9cDg|]}tdtz  |zSrr+r,s r7rz*fresnelc._eval_aseries.. r-r9cBg|]}tdtz |zSrr+r,s r7rz*fresnelc._eval_aseries.. s&***1$qt**Q,***r9r)rrrrLr+rrr#rr$rrr rrr/s ` @r7rzfresnelc._eval_aseriess s,,,,,,a QZ!*- - - ! A33333((333AAaC 666661++6666A+****A*****Aaj((bAQV8s1a4yya03q!t99S!W3DD$q$qt**Q,''(*/%!Q$*:*:; ;ww$$Qq$777r9)rrrrr#rrrMrrrrcrrrrnrrrrs@r7rr sSShI EE  ^ ^ W\ ^aaaUUU[[[333 # # #888888888r9rcFeZdZdZedZfdZddZdZxZ S)rzi Helper function to make the $\mathrm{erf}(z)$ function tractable for the Gruntz algorithm. c,|jr tjSdSr)rPrrM)rWrXs r7rZz _erfs.eval s ; 5L  r9cZ ddlm}|d}|tjura|jd fdt |D}|d d|zdzzz |}t |||||zS|t} | tjura|jd fdt |D}|d d|zdzzz |}t |||||zSt ||||S)Nrrcg|]`}dttz td|zztd | zzt|z dz d|zdzzzaSr<r*rrr rrrs r7rz'_erfs._eval_aseries.. NNN?@488i!nn,q00/r.$Q<<(+,Q3!A#'*:;NNNr9r<r*cg|]`}dttz td|zztd | zzt|z dz d|zdzzzaSrArBrs r7rz'_erfs._eval_aseries.. rCr9) rrrrLr+rrrvrTr rr) r1r`rr4rrrlorYrurs @r7rz_erfs._eval_aseries sr,,,,,,a AJ   ! ANNNNDI!HHNNNAaAaC!G na((AG**1a66: :  * *1 - -  ?? ! ANNNNDI!HHNNNAaAaC!G na((AG**1a66: :ww$$Qq$777r9r<c|dkr:|jd}dttz d|zt|zzSt ||)Nr<rr]r*)r+rr rrr1rArus r7rBz _erfs.fdiff sI q== ! Ad2hh;1U1XX- -$T844 4r9c `tjt|z t|dzzSr)rrMr;rrs r7_eval_rewrite_as_intractablez"_erfs._eval_rewrite_as_intractable s#AAqD ))r9r) rrrrrrZrrBrJrrs@r7rr s [8888865555*******r9rcFeZdZdZfdZd dZdZfdZd fd ZxZ S) rez~ Helper function to make the $\mathrm{Ei}(z)$ and $\mathrm{li}(z)$ functions tractable for the Gruntz algorithm. cVddlm}|dtjtjfvr$t ||||S|jdfdt|D}|d|dzzz |}t| ||||zS)NrrcFg|]}t|dz |dzzzSrr{rs r7rz&_eis._eval_aseries.. s0 = = =qYq\\QqSAEN * = = =r9r<) rrrrLrNrrr+rrrv) r1r`rr4rrrErFrurs @r7rz_eis._eval_aseries s,,,,,, 8AJ(:; ; ;77((E1d;; ; IaL = = = =E!HH = = = E!AAJ, " "Q&&q!T22Q66r9r<c|dkr,|jd}tj|z t|z St ||)Nr<r)r+rrMrerrHs r7rBz _eis.fdiff s= q== ! A519tAww& &$T844 4r9c Bt| t|zSr)rrQrs r7rJz!_eis._eval_rewrite_as_intractable sA2wwr!uu}r9c|jd|d}|jr'|j|j}||||St |||S)Nrr)r+rrPrJrr)r1r4rrrrrxrs r7rz_eis._eval_as_leading_term st Yq\  1 % % : D1149=A**14d*CC Cww,,QT,EEEr9rc|jd|d}|jr&|j|j}||||St |||Sre)r+rrPrJrvrrws r7rvz_eis._eval_nseries sk Yq\  1 % % : /1149=A??1a.. .ww$$Q4000r9rr|) rrrrrrBrJrrvrrs@r7rere s 7 7 7 7 75555FFFFF1111111111r9reN)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;rrrrFrRrSrQrrr]rrrrrarbrrrrrerr9r7res=FF"!!!!!$$$$$$OOOOOOOOOO%%%%%%777777777777'''''' """"""::::::::&&&&&&[[[[[[[[[[LLLLLLLLLL>>>>>>>>????????FFFFFFFFFF<<<<<<<<CCCCCCCCCC88888888    *l-l-l-l-l-/l-l-l-^HHHHH?HHHD}B}B}B}B}B?}B}B}B@O!O!O!O!O!?O!O!O!bZ$Z$Z$Z$Z$_Z$Z$Z$zP;P;P;P;P;P;P;P;f[[[[[o[[[Bt@t@t@t@t@t@t@t@n}?}?}?}?}?_}?}?}?@   FeeeeeeeeNZ+Z+Z+Z+Z+Z+Z+Z+@<J<J<J<J<JO<J<J<J~EEEEE EEEPJ@J@J@J@J@ J@J@J@Xggggg gggTnnnnn nnnj5-5-5-5-5-o5-5-5-p^8^8^8^8^8^8^8^8B^8^8^8^8^8^8^8^8L.*.*.*.*.*O.*.*.*b*1*1*1*1*1?*1*1*1*1*1r9