x = 0:0.05:2;
t = 0:0.05:4;
for k=1:length(t)
	u = (intpsi1(x-t(k),x+t(k)))/2 + (phi1(x+t(k)) + phi1(x-t(k)))/2;
	plot(x, u);
	axis ([0, 4, -1, 1])
	M(k) = getframe;
end
movie (M,3)