from scipy import *

""" Integrates bessel function j1(t)/t between 0 and x
    for x between 0 and 10.
"""

def func(x):
    return integrate.quad(lambda t: special.j1(t)/t, 0, x)

x = r_[0:30:0.1]

for tx in x:
    print tx, func(tx)[0]