clear all
more off

gridsize=1000;
dx=1/gridsize;
for i=0:gridsize
    xi=i*dx;
    x(i+1)=xi;
    y(i+1)=f(xi);
end
