THETA=pi/4;
Xs1=rb.*sin(tan(u)+etab)-rb.*tan(u).*cos(tan(u)+etab);
Ys1=rb.*cos(tan(u)+etab)+rb.*tan(u).*sin(tan(u)+etab);
XT1=Xs1+ST1.*cos(tan(u)+etab);
YT1=Ys1-ST1.*sin(tan(u)+etab); %齿轮端面
%空间螺旋面
XTK1=XT1.*cos(THETA)-YT1.*sin(THETA);
YTK1=XT1.*sin(THETA)+YT1.*cos(THETA);
ZTK1=THETA.*19*5/2/sin(20*pi/180);
plot3(XTK1,YTK1,ZTK1)
出错 lilunkongjianchixing20160412 (line 23)
plot3(XTK1,YTK1,ZTK1)