function y=fun(t)
y=1.0 / (1.0+25*t.*t);
