clear all
more off

J=500;
N=100;
dx=1/J;
T=1;
dt=T/N;

for j=0:J-1
    x(j+1)=dx*j;
    y(j+1)=g(dx*j);
end
plot(x,y)
pause
for n=1:N
    tn=n*dt;
    y=coneq(y,tn-dt,tn,10);
    plot(x,y);
    pause
end
