In [3]:
from scipy.integrate import quad
from scipy.special import airy
from matplotlib import pyplot as plt
from numpy import linspace
N=100000
def f(x):
return airy(x)[0]
x=linspace(-100,10,N+1)
y=f(x)
%matplotlib qt
plt.plot(x,y)
Out[3]:
In [ ]: