global m;

x=-30:0.25:30;
m = 241;

T = 4;

u0 = KdVinit(x,1);

plot(x,u0);
pause

options=odeset('AbsTol',1e-10,'RelTol',1e-9);
[tt, yy]=ode15s('KdVFunc',[0 T],u0,options);

for i=1:size(yy,1)  
    plot(x,yy(i,:)');
    drawnow;
%    pause
end

mesh(x,tt,yy);
xlabel('rom');
ylabel('tid');