Example 2 function cycloid(r,n) 号w=1rad/s theta =0:pi/100:2*pi; center y =r; cyc_x =[] cyc_y [] jj=1; for t 0:pi/30:n*2*pi center x =r*t; x =r*cos(theta)+center x; y =r*sin(theta)+center y; current x center x-r*sin(t); current y center_y-r*cos(t); cyc x [cye x,current x]; cye_y [cyc_y,current y]; plot (x,y,cyc_x,cyc_y,current x,current y,'rd'); axis equal axis ([-r,2*pi*r*n+r,0,3*r]) mov(jj)=getframe; jj=jj+1; end movie (mov) returnExample 2 function cycloid(r,n) % w=1rad/s theta = 0:pi/100:2*pi; center_y = r; cyc_x = []; cyc_y = []; jj = 1; for t = 0:pi/30:n*2*pi center_x = r*t; x = r*cos(theta) + center_x; y = r*sin(theta) + center_y; current_x = center_x-r*sin(t); current_y = center_y-r*cos(t); cyc_x = [cyc_x, current_x]; cyc_y = [cyc_y, current_y]; plot(x,y,cyc_x,cyc_y,current_x,current_y,'rd'); axis equal axis([-r,2*pi*r*n+r,0,3*r]) mov(jj) = getframe; jj = jj+1; end movie(mov) return