clear;clc;
R=10;
Rr=4;
thita0=pi/180*5;


[xyz]=sphere（10）;
axis（[-15?15?-15?15?-5?5]‘equal‘‘manual‘）;hold?on;grid?on;
surf（3*x3*y3*z）;%axis?off
xlabel（‘x‘）
ylabel（‘y‘）
zlabel（‘z‘）
h1=surf（xyz）;
h2=surf（xyz）;
thita1=pi/180*30/360;thita2=pi/180*30;


stop?=?uicontrol（‘style‘‘toggle‘‘units‘‘norm‘‘pos‘[.91?.94?.08?.05]?...
???‘backgr‘‘w‘‘fontw‘‘bold‘‘string‘‘stop‘）;
k=0;
while?get（stop‘value‘）?==?0
????x1=2*x+R*cos（k*thita1）;y1=2*y+R*sin（k*thita1）;z1=2*z+0;
x2=x+R*cos（k*thita1）+Rr*sin（k*th