matlab 凸轮程序,求助，关于凸轮的一段matlab程序。

在书上抄了一段matlab设计凸轮轮廓的一段程序，把它制成M文件，但运行不了，求各位高手帮帮忙啊，谢谢了！

程序如下：

function RollerCam

beta=60*pi/180;rb=3;h=0.5;rc=0.5;

rf=0.375;m=.375;n=23;

phi=linspace(0,beta,n);

ph=[phi beta+phi];

[Rx,Ry]=ContourRoller(0,rb,h,beta,rc,m,rf);

rb=sqrt(Rx(1)^2,Ry(1)^2);

delta=atan2(Ry(1),Rx(1));

[Rx,Ry,Cx,Cy]=ContourRoller(phi,rb,h,beta,rc,m,rf);

ang=linspace(2*beta+delta,2*pi+delta,40);

plot(Rx,Ry,'k',Rx(1)*cos(ang),Rx(1)*sin(ang),'k')

hold on

plot(0,0,'k+',Cx(1:5:2*n),Cy(1:5:2*n),'k+')

axis equal

phd=linspace(0,2*pi,50);

[x,phx]=meshgrid(Cx(1:5:2*n),phd);

y=meshgrid(Cy(1:5:2*n),phd);

hold on

plot(x+rc.*cos(phx),y+rc.*sin(phx),&#