function [x,yn]=xtgrid(gridsize,timesteps);
  dx=1/gridsize;
  for i=0:gridsize
    xi=i*dx;
    x(i+1)=xi;
    y(i+1)=f(xi);
  end
  yn=heateq(y,0.0,0.1,timesteps);
end