求助：matlab里球面坐标系下的非圆齿轮渐开线命令

爱无止境

我想根据一篇文献里的关系式重现齿轮圆锥曲线，前面公式及命令如下：
9 H; u7 l, y( ~fi1=[0:0.1:2*pi];
" [( T. j8 B9 h/ f& t- I) |& a4 D; jc=0.16;
1 U4 z, D+ B, g+ RR=43.6;
fi2=0.75*fi1+(c/3)*0.75*cos(3*fi1+3*pi/2); %0.75=z1/z2=3/4,我直接换成0.75代入关系式了。
( I8 g) b7 O/ adelta1=atan(0.75*(1-c*sin(3*fi1+3*pi/2)));
$ ^# a/ U8 ^" q( Hdelta2=acot(0.75*(1-c*sin(3*fi1+3*pi/2)));

5 ?, n/ x8 l) y' } s! V球面方程如下：
2 u0 `7 T# Q7 {2 H, ]x1=R*(sin(delta1).*cos(fi1));
) m- }; I- n! z9 c( Qy1=R*sin(delta1).*sin(fi1);
# f9 ?$ h5 p( Cz1=R*cos(delta1);
0 ?# l) h6 \' [& p" M4 o/ X3 Zr1=x1*i+y1*j+z1*k%这里的i,j,k前的系数是节曲线上的点在三个坐标上的坐标。写到程序里面的不是这样的，下面也是一样。这样写出来大家理解直观一些。
x2=R*(sin(delta2).*cos(fi2));
y2=R*sin(delta2).*sin(fi2);
z2=R*cos(delta2);
r2=x2*i+y2*j+z2*k
2 L% K# Y+ p1 l
r1和r2的值验算出来的都是43.6，没有问题。但是我是不知道用什么样的命令得到下面的球面图形。由于本人对matlab不很了解，到这里用命令老是出错，各位英雄女侠帮帮忙。感激不尽！！！
6 C. p/ R# @$ c7 j3 D9 R
# q9 F& g# I3 z6 `
0 Q4 n7 [! }6 Z0 z4 ~% ^: Q