bessel.py

import pylab
from numpy import *
from scipy import *
from scipy.special import *
from matplotlib import*

x = arange(0,100,0.1)
y = i1(x)
pyplot.plot(x, y)
pyplot.show()