用MATLAB画收敛域的程序,急!!!!牛顿法求收敛域,matlab编程遇到问题

用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

运行不出图来,跪求高手


本文来自互联网用户投稿,文章观点仅代表作者本人,不代表本站立场,不承担相关法律责任。如若转载,请注明出处。 如若内容造成侵权/违法违规/事实不符,请点击【内容举报】进行投诉反馈!

相关文章

立即
投稿

微信公众账号

微信扫一扫加关注

返回
顶部