第4题
theta = linspace(0,pi,100);
phi = linspace(0,2*pi,100);
[Theta,Phi] = meshgrid(theta,phi);
r = sqrt( 0.8 + 0.5*( cos(2*Phi)-1 ).*( cos(4*Theta)-1 ) );
x = r.*sin(Theta).*cos(Phi);
y = r.*sin(Theta).*sin(Phi);
z = r.*cos(Theta);
figure
surf(x,y,z)
第5题
x = 0:pi/100:2*pi;
y = sin(pi/3) .* ( 0<=x & x < 2*pi/3 ) + sin(x) .* ( 2*pi/3 <= x & x < pi ) + cos(2*pi/3) .* ( pi <= x & x <= 2*pi ) ;
figure
plot(x,y)