求助：matlab里球面坐標系下的非圓齒輪漸開線命令

愛無止境

我想根據一篇文獻里的關系式重現齒輪圓錐曲線，前面公式及命令如下：
fi1=[0:0.1:2*pi];
6 B& F; p) [; A4 z, a1 [c=0.16;
: x- G) F: G2 C( x% a/ G7 JR=43.6;
, o' Z' O, {' m4 I6 M7 X3 j8 zfi2=0.75*fi1+(c/3)*0.75*cos(3*fi1+3*pi/2);      %0.75=z1/z2=3/4,我直接換成0.75代入關系式了。
- ?3 W. S( W' [" v/ o* m6 n2 Ndelta1=atan(0.75*(1-c*sin(3*fi1+3*pi/2)));
+ |# g* ~  |; M& l5 _3 bdelta2=acot(0.75*(1-c*sin(3*fi1+3*pi/2)));

( g8 |& `7 z" k3 i" `球面方程如下：
x1=R*(sin(delta1).*cos(fi1));
y1=R*sin(delta1).*sin(fi1);
- C4 w  r! H) ^8 p; Dz1=R*cos(delta1);
' N+ e5 `9 `. vr1=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);
1 [. r5 |: J4 z! L# u- g0 Dr2=x2*i+y2*j+z2*k

1 R% k8 j  p7 Z. `4 kr1和r2的值驗算出來的都是43.6，沒有問題。但是我是不知道用什么樣的命令得到下面的球面圖形。由于本人對matlab不很了解，到這里用命令老是出錯，各位英雄女俠幫幫忙。感激不盡！！！
8 Z- W) m" P) w; |0 d  v+ _
+ u) C3 Z" g9 Q# e7 U
4 S( }8 \6 w+ p