708.sagews

from numpy import linspace
from scipy.special import sici
import pylab

r = 15
x = linspace(-r, r)
y = sici(x)[0]

pylab.plot(x,y)
pylab.grid(True)
pylab.savefig('sin-int.png')