[求助]曲线的三维作图

[zjgydxcbt]

怎么把这个程序的月球轨道转化三维曲线？？？？（另外画,与此程序无关）（z坐标不变）
function solar_system
% 模拟太阳系运动
t=linspace(0,2*pi,100);
fill(cos(t),sin(t),'r');
hold on;
plot(4*cos(t),sin(t)*4,'k');
set(gca,'position',[0 0.11 0.775 0.815]) %太阳和地球的轨道
a=0.1;b=0;
xe=4*cos(a)+cos(t)*0.6;
ye=4*sin(a)+sin(t)*0.6;
He=fill(xe,ye,'b'); %地球
xm=4*cos(a)+cos(b);
ym=4*sin(a)+sin(b);
set(gcf,'doublebuffer','on');
Hm=plot(xm,ym,'c.','markersize',24); %月亮
aa=gca;
axis([-6,6,-6,6]);
axis square; %静态日地月
k=1;da=0.1;db=0.5;
xlabel('Please press "space" key and stop this program!',...
'fontsize',12,'color','r');
title('simulate solar system') %显示字
axes('position',[0.75,0.11,0.25,0.8]);
fill(0.2+cos(t)*0.18,0.75+sin(t)*0.08,'r');
ylim([0,1]);xlim([0,0.9]);
text(0.5,0.75,'Sun');hold on;
fill(0.2+cos(t)*0.11,0.5+sin(t)*0.05,'b');
text(0.5,0.5,'Earth');
plot(0.2,0.3,'c.','markersize',24);
text(0.5,0.3,'Moon');
axis off
axes(aa); %右边图像说明
while k;
s=get(gcf,'currentkey');
if strcmp(s,'space');
clc;k=0;
end %按SPACE停止
a=a+da;
b=b+db;
xe=4*cos(a)+cos(t)*0.6;
ye=4*sin(a)+sin(t)*0.6;
xm=4*cos(a)+cos(b);
ym=4*sin(a)+sin(b);
set(He,'xdata',xe,'ydata',ye);
set(Hm,'xdata',xm,'ydata',ym);
pause(0.1);
if a<80;
plot(xm,ym);
end
end
figure(gcf);