import numpy as np
from scipy import integrate
from scipy import special

I_quad = integrate.quad(np.sin, 0, np.pi)

print(I_quad[0])

J = special.jv(0,2.4)

print(J)