exm01033_4.m

来自「这本书是matlab经典书籍」· M 代码 · 共 15 行

a=-8;
r=a^(1/3) 

m=[0,1,2];						%为三个方根而设
R=abs(a)^(1/3);				%模的开三方
Theta=(angle(a)+2*pi*m)/3;	%-pi<Theta<=pi的三个相角
rrr=R*exp(i*Theta)

t=0:pi/20:2*pi;x=R*sin(t);y=R*cos(t);
plot(x,y,'b:'),grid									%画一个半径为R的圆
hold on
plot(rrr(1),'.','MarkerSize',50,'Color','r')		%画第一象限的方根
plot(rrr([2,3]),'o','MarkerSize',15,'Color','b')	%画另二个方根
axis([-3,3,-3,3]),axis square						%保证屏幕显示呈真圆
hold off