Le dimanche 11 mars 2012 à 01:31 -0600, Val Kalatsky a écrit :
Can you?
The question should be: Why sympy does not have Fresnel integrals?
Well, the development version has hypergeometric functions, which is
enough to integrate sin(x**2):
from sympy import *
x = Symbol('x')
integrate(sin(x**2))
x**3*gamma(3/4)*hyper((3/4,), (3/2, 7/4), -x**4/4)/(4*gamma(7/4))
integrate(sin(x**2), (x, 0, oo))
sqrt(2)*sqrt(pi)/4
It can also create a numerical function for the result, but it uses
mpmath, so it won't play well with numpy:
f = lambdify(x, integrate(sin(x**2)))
f(4567.)
mpf('0.62676518399179759')
f(np.linspace(-5., 5., 10))
Traceback (most recent call last):
File ipython-input-35-97e0996f77c3, line 1, in module
f(np.linspace(-5., 5., 10))
File string, line 1, in lambda
File /home/ronan/dev/sympy/sympy/mpmath/functions/hypergeometric.py,
line 199, in hyper
z = ctx.convert(z)
File /home/ronan/dev/sympy/sympy/mpmath/ctx_mp_python.py, line 662,
in convert
return ctx._convert_fallback(x, strings)
File /home/ronan/dev/sympy/sympy/mpmath/ctx_mp.py, line 556, in
_convert_fallback
raise TypeError(cannot create mpf from + repr(x))
TypeError: cannot create mpf from array([ -1.5625e+02, -5.71797363e
+01, -1.48843545e+01,
-1.92901235e+00, -2.38149672e-02, -2.38149672e-02,
-1.92901235e+00, -1.48843545e+01, -5.71797363e+01,
-1.5625e+02])
On Sun, Mar 11, 2012 at 1:06 AM, aa telukp...@gmail.com wrote:
why sympy cannot integrate sin(x**2)??
___
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion
___
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion
___
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion
