figure(1);
F=@(x) 54000*x+410000*(x-0.024).*(x>0.024);
t=linspace(0,0.03);dt=t(2)-t(1);
z(1)=0;z(2)=3*dt;
for i=2:99;
z(i+1)=2*z(i)-z(i-1)-F(z(i))/5*dt^2;
end;
for i=1:100;
w=z(i);
set(h(1),'ydata',yy{1}-w);
set(h(2),'ydata',yy{2}+linspace(0,-w));
set(h(3),'ydata',yy{3}+linspace(0,0.024-w).*(w>0.024));
set(h(5),'ydata',yy{5}-w);
pause(0.1);
end;