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

愛無止境

我想根據(jù)一篇文獻里的關系式重現(xiàn)齒輪圓錐曲線，前面公式及命令如下：
fi1=[0:0.1:2*pi];
c=0.16;
1 @2 K# `2 E6 X* b( n; v* yR=43.6;
fi2=0.75*fi1+(c/3)*0.75*cos(3*fi1+3*pi/2);      %0.75=z1/z2=3/4,我直接換成0.75代入關系式了。
3 F6 ?3 g; w# N% S/ Bdelta1=atan(0.75*(1-c*sin(3*fi1+3*pi/2)));
. Q3 P  E( D  B' g5 rdelta2=acot(0.75*(1-c*sin(3*fi1+3*pi/2)));
; ~) A& U. h- `( y
球面方程如下：
4 f4 ~4 V7 d$ n/ w$ sx1=R*(sin(delta1).*cos(fi1));
" q; b3 e: t) z2 t. c: s$ Ry1=R*sin(delta1).*sin(fi1);
z1=R*cos(delta1);
8 a: _/ x0 j/ ]# \r1=x1*i+y1*j+z1*k     %這里的i,j,k前的系數(shù)是節(jié)曲線上的點在三個坐標上的坐標。寫到程序里面的不是這樣的，下面也是一樣。這樣寫出來大家理解直觀一些。
- W" b* W. H7 v5 D5 Px2=R*(sin(delta2).*cos(fi2));
  G2 R0 ]5 s0 O" G2 {' L! Ly2=R*sin(delta2).*sin(fi2);
z2=R*cos(delta2);
9 v4 H2 |9 r9 a8 E5 ]$ q: Yr2=x2*i+y2*j+z2*k
3 G4 @# n, ?" ?# y* V$ t1 L
r1和r2的值驗算出來的都是43.6，沒有問題。但是我是不知道用什么樣的命令得到下面的球面圖形。由于本人對matlab不很了解，到這里用命令老是出錯，各位英雄女俠幫幫忙。感激不盡！！！
1 e2 `: M/ ]; W% F5 j+ D% o
2 @5 Y( E. y" M3 s