On Mon, 19 Feb 2024 at 20:25, Mark “Essa King” Sukaiti < [email protected]> wrote:
> D(r)=-1/2*((3*r^2 - 1)*elliptic_ec((4*r/(r^2 + 2*r + 1)))^2 + (r^2 - 2*r + > 1)*elliptic_ec((4*r/(r^2 + 2*r + 1)))*elliptic_kc((4*r/(r^2 + 2*r + > 1))))/(pi^2*r^8 - 2*pi^2*r^7 - pi^2*r^6 + 4*pi^2*r^5 - pi^2*r^4 - > 2*pi^2*r^3 + pi^2*r^2) > > D(r).limit(r=0) > > The limit should be -0.125 (or -1/8) but it seems maxima doesnt know the > limit of elliptic_ec(x) for x->0 > > Sympy also fails giving -5/8. > To follow up on this the SymPy bug is now fixed in master: https://github.com/sympy/sympy/pull/26264 In [1]: e = ((1 - 3*x**2)*elliptic_e(4*x/(x**2 + 2*x + 1))**2/2 - (x**2 - 2*x + 1)*elliptic_e(4*x/(x* ...: *2 + 2*x + 1))*elliptic_k(4*x/(x**2 + 2*x + 1))/2)/(pi**2*x**8 - 2*pi**2*x**7 - pi**2*x**6 + ...: 4*pi**2*x**5 - pi**2*x**4 - 2*pi**2*x**3 + pi**2*x**2) In [2]: e.limit(x, 0) Out[2]: -1/8 Oscar -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/sage-support/CAHVvXxQst7%3DUUOYd%3DcRAGaMfahH_6LtvV73kwKPvzaC1Dy_Rtw%40mail.gmail.com.
