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

爱无止境

我想根据一篇文献里的关系式重现齿轮圆锥曲线，前面公式及命令如下：
fi1=[0:0.1:2*pi];
* I  h. u, V& Z' t) ^4 e6 _! p! Uc=0.16;
% e  O8 K% a# ~, l( n0 h2 sR=43.6;
! u' [. s; z4 ^3 n; Efi2=0.75*fi1+(c/3)*0.75*cos(3*fi1+3*pi/2);      %0.75=z1/z2=3/4,我直接换成0.75代入关系式了。
8 j" A, v! V. x3 |& f& j( A. @delta1=atan(0.75*(1-c*sin(3*fi1+3*pi/2)));
delta2=acot(0.75*(1-c*sin(3*fi1+3*pi/2)));
; J9 S: Z8 k2 s2 j
球面方程如下：
x1=R*(sin(delta1).*cos(fi1));
y1=R*sin(delta1).*sin(fi1);
z1=R*cos(delta1);
7 X3 K" y4 z2 j! ]! c- Vr1=x1*i+y1*j+z1*k     %这里的i,j,k前的系数是节曲线上的点在三个坐标上的坐标。写到程序里面的不是这样的，下面也是一样。这样写出来大家理解直观一些。
x2=R*(sin(delta2).*cos(fi2));
y2=R*sin(delta2).*sin(fi2);
+ R6 L$ e7 q0 ]z2=R*cos(delta2);
- G2 u* \- g& e! E) E: s+ vr2=x2*i+y2*j+z2*k

3 M( W  A: p% _3 _  Nr1和r2的值验算出来的都是43.6，没有问题。但是我是不知道用什么样的命令得到下面的球面图形。由于本人对matlab不很了解，到这里用命令老是出错，各位英雄女侠帮帮忙。感激不尽！！！
& E' I' c( K( ^" f/ F/ U) E  d5 W( g/ p
6 f% a4 e+ Q$ ^$ d% |% s  y- U
+ ?9 u2 @9 k! }; D  m: F% A