matplotlib Code

Brought to you by: cjgohlke, dsdale, efiring, heeres, and 8 others

[r3313]: / trunk / py4science / book / examples / scipy / example4.2 Maximize Restore History

12 lines (9 with data), 321 Bytes

>>> result = integrate.quad(lambda x: special.jv(2.5,x), 0, 4.5)
>>> print result
(1.1178179380783249, 7.8663172481899801e-09)

>>> I = sqrt(2/pi)*(18.0/27*sqrt(2)*cos(4.5)-4.0/27*sqrt(2)*sin(4.5)+
    sqrt(2*pi)*special.fresnl(3/sqrt(pi))[0])
>>> print I
1.117817938088701

>>> print abs(result[0]-I)
1.03761443881e-11 