(Put your own values for rp and epsl which pressure angle result less than 20)

%Problem 8.29 % cycloidal displacement cam function %lift h=20 mm % beta = 60 degrees Rp=input('the minimum radius of the pitch rp='); h=input('the total lift (rise) of the segment h='); beta=input('the total angle of the segment beta='); epsl=input('the eccentricity epselon='); 'thea displ veloc accel radius(R) curve press.ang.' format short for t=0:3:60 theta=t*pi/180; s=h*(theta/beta)-(1/(2*pi))* sin(2*pi*theta/beta); v=(h/beta)*(1-cos(2* pi*theta/beta)); a=((2 *pi*h)/(beta^2))*sin(2*pi*theta/beta); R=Rp+s; rho=((R^2+v^2)^(3/2))/(R^2+2*v^2-a*R); phi (180/pi)*atan((v-epsl)/(s+sqrt(Rp^2-epsl^2))); disp ([t,s,v,a,R,rho,phi]); end