x = 0:2000;
ss1 = sin(0.1*x)+sin(0.105*x);
plot(ss1)
title('sum of two sin waves')
pause


x = -2000:2000;
ss1 = sin(0.1*x);
for i=1:50
	ss1 = ss1 +sin( (0.1+0.002*i)*x);
end
plot(ss1)
title('sum of 50 sin waves')
pause



moviein(12);
for t=1:13 
        plot(sin(10*t-0.1*x)+sin(10.5*t-0.105*x))
	axis([0,2000,-4,4])
	title('sum of two sins, dw/w=0.05, Vp=100, vg=100')
        m(:,t) = getframe;
end
movie(m,20)
pause

moviein(12);
for t=1:13 
        plot(sin(0.5*t-0.1*x) )
	axis([0,2000,-4,4])
	title('traveling sin wave')
        m(:,t) = getframe;
end
movie(m,20)
pause

moviein(12);
for t=1:13 
        plot(sin(0.5*t-0.1*x)+sin(1.0*t-0.125*x))
	axis([0,2000,-4,4])
	title('sum of two sins, dw/w=0.5, Vp=5, vg=20')
        m(:,t) = getframe;
end
movie(m,20)
pause

moviein(12);
for t=1:13 
        plot(sin(1.0*t-0.1*x)+sin(0.5*t-0.12*x))
	axis([0,2000,-4,4])
	title('sum of two sins, dw/w=0.0, Vp=10, vg=-25')
        m(:,t) = getframe;
end
movie(m,20)
end
