function [x,t]=rk(x,t,h)
  k1=h*f(x,t);
  k2=h*f(x+0.5*k1,t+0.5*h);
  k3=h*f(x+0.5*k2,t+0.5*h); 
  k4=h*f(x+k3,t+h);
  x=x+1/6*(k1+2*k2+2*k3+k4);
  t=t+h;