用5种颜色画出x^5-1=0的收敛域.
不怎么会matlab,老师布置的作业就硬着头皮做.....下面是我编的程,但是运行不出结果,求告手指点!!!
a=[1,0,0,0,0,1];x=roots(a)
A=x(1)
B=x(2)
C=x(3)
D=x(4)
E=x(5)
h=0.01
for a=-2:h:2
for b=-2:h:2
z=a+b*i;
y=newton(z);
if (abs(y-A)<1.0e-6)
plot(a,b,'r');
hold on
elseif (abs(y-B)<1.0e-6)
plot(a,b,'y');
hold om
elseif (abs(y-C)<1.0e-6)
plot(a,b,'g');
hold om
elseif (abs(y-D)<1.0e-6)
plot(a,b,'b');
hold om
elseif (abs(y-E)<1.0e-6)
plot(a,b,'p');
hold om
end
end
end
function y=newton(z)
if (z==0)
y=0;
return;
end
for l=1:1:100000
y=z-(z^5-1)/(5*z^4);
if(abs(y-z)<1.0e-7)
break;
end
z=y;
end
运行不出图来,跪求高手