set (0, "DefaultAxesFontName", "Arial")
x = linspace(0, 1, 10);
figure(1);
psi_0 = sin(2 * pi * x);
%plot(x, psi_0);
%print -deps psi_0.eps;
psi_0_06 = iterate(psi_0, 0.1, 0.02, 3)
%subplot(5, 1, 2);
xx = linspace(0, 1, 100);
u = exp(-4*pi*pi*(1/6)*0.06)*sin(2*pi*xx);
plot(x, psi_0_06, xx, u);
%print -deps psi_0_06.eps;
psi_0_1 = iterate(psi_0_06, 0.1, 0.02, 2)
u = exp(-4*pi*pi*(1/6)*0.1)*sin(2*pi*xx);
plot(x, psi_0_1, xx, u);
%print -deps psi_0_1.eps;
psi_0_9 = iterate(psi_0_1, 0.1, 0.02, 40)
u = exp(-4*pi*pi*(1/6)*0.9)*sin(2*pi*xx);
plot(x, psi_0_9, xx, u);
%print -deps psi_0_9.eps;
psi_50 = iterate(psi_0_9, 0.1, 0.02, 2455)
u = exp(-4*pi*pi*(1/6)*50)*sin(2*pi*xx);
plot(x, psi_50, xx, u);
%print -deps psi_50.eps;