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

爱无止境

我想根据一篇文献里的关系式重现齿轮圆锥曲线，前面公式及命令如下：
- @6 w7 o3 M# J- C# sfi1=[0:0.1:2*pi];
c=0.16;
/ n" P9 h1 Y0 p2 {9 uR=43.6;
fi2=0.75*fi1+(c/3)*0.75*cos(3*fi1+3*pi/2);      %0.75=z1/z2=3/4,我直接换成0.75代入关系式了。
9 T0 C4 c9 i+ @# r+ U& gdelta1=atan(0.75*(1-c*sin(3*fi1+3*pi/2)));
; q4 I2 \, q) `6 c0 y1 E3 Idelta2=acot(0.75*(1-c*sin(3*fi1+3*pi/2)));
# J2 Y/ T) v9 a7 R
, _: S5 _. _0 b3 e  h球面方程如下：
x1=R*(sin(delta1).*cos(fi1));
y1=R*sin(delta1).*sin(fi1);
" X" }% |; x' |: c3 W) dz1=R*cos(delta1);
0 @- v0 V6 h7 j# s5 @0 E- `r1=x1*i+y1*j+z1*k     %这里的i,j,k前的系数是节曲线上的点在三个坐标上的坐标。写到程序里面的不是这样的，下面也是一样。这样写出来大家理解直观一些。
0 @1 k- {* K' i4 w8 L$ |7 Q+ Jx2=R*(sin(delta2).*cos(fi2));
y2=R*sin(delta2).*sin(fi2);
z2=R*cos(delta2);
r2=x2*i+y2*j+z2*k
$ o: Y7 `' n* F$ G
0 v( s* m/ l# {, Z1 S  G3 Ur1和r2的值验算出来的都是43.6，没有问题。但是我是不知道用什么样的命令得到下面的球面图形。由于本人对matlab不很了解，到这里用命令老是出错，各位英雄女侠帮帮忙。感激不尽！！！
. d1 }' o; I1 x% V
( \6 X8 s2 c( e% O" |# h, X# ~