% Bouncing Ball x = 0; y = 2; dx = 2; dy = 3; g = -9.8; t = 0; dt = 0.01; Bounce = 0; N = 0; while(Bounce < 2) ddx = 0; ddy = g; dx = dx + ddx*dt; dy = dy + ddy*dt; t = t + dt; x = x + dx*dt; y = y + dy*dt; if(y<0) dy = abs(dy); Bounce = Bounce + 1; end N = mod(N+1,5); if(N==0) plot(x,y,'ro'); hold on xlim([0,10]); ylim([0,5]); pause(0.01); end end