th=(0:.25:179.75)*pi/180;th=th+(th==0)*eps;
fy=(0:.5:359.5)*pi/180;     %球坐标中的th为0--180度，fy为0--360度
kl=2*pi*0.65;
F=(cos(kl*cos(th))-cos(kl))./sin(th)/(1-cos(kl));
r=abs(F);    
x=r.*sin(th).*cos(fy);     %将球坐标中的r,th,fy转化为直角坐标
y=r.*sin(th).*sin(fy);
z=r.*cos(fy);
%[X,Y]=meshgrid(x,y);
plot3(x,y,z);
我希望出现较完整的三维图形，可是不行，能不能帮我改一改？谢谢！