{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

hw2_finished

# hw2_finished - Boolean turned on for making a movie...

This preview shows pages 1–3. Sign up to view the full content.

% Boolean turned on for making a movie: Make_Movie=1; % Geometry of the closed container: Right_Wall = 1; Left_Wall = 0; Bottom_Wall = 0; Top_Wall = 1; IS_WALL_RIGHT = 1; IS_WALL_LEFT = 1; IS_WALL_BOTTOM = 1; IS_WALL_TOP = 1; % Physical parameters: g =.0981; alpha=.8 ; beta =.99 ; % Acceleration: acceleration_x= 0; acceleration_y= -g; % Parameters for time evolution: final_time=100; dt =.03; % Initial time: t=0; % Initial velocity vector: v_hor= .3; v_ver= 0; % Original position of the ball: Radius=.05 ; x=.5 ; y= 1-Radius; count=0; if Make_Movie fig=figure; aviobj = avifile( 'Bouncing_Balll.avi' ); end while t<final_time && count < 2000 if t+dt>final_time dt=final_time-t; end t=t+dt; v_hor = RK2(dt, v_hor, acceleration_x); v_ver = RK2(dt, v_ver, acceleration_y); % Treatment at the bottom wall:

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
if ( RK2(dt, y, v_ver ) ) < Bottom_Wall + Radius && IS_WALL_BOTTOM % would cross the wall. % Position the object on the wall and Draw: x = RK2( dt, x, v_hor ); y = Bottom_Wall + Radius; Draw_Disk(x,y,Radius); axis([0 1, 0 1]); pause(dt/10); if Make_Movie F = getframe(fig); aviobj = addframe(aviobj,F); end % Define the velocity to reflect how the ball is bouncing off:
This is the end of the preview. Sign up to access the rest of the document.

{[ snackBarMessage ]}