% Resonanz bei Schwingungen 
clf; axis([0 1 -2 2]); hold on; 
t = [0:0.001:1];

for k=0:100
   w = 10-k/10; 
   p = sin(2*pi*(10-w)*t); 
   q = sin(2*pi*(10+w)*t);
   h = plot(t,p,t,q,t,p+q);
   M(k+1) = getframe;
   delete(h); 
end

hold off;
movie(M,1,10);