bisection methods?

[木心]

好像图做出来了...但结果好像有点乱...我再查一遍...
谢谢

y=zeros(1,200);
i=1;
for x=0:0.01:2
a=-2;b=2; % try 200 x, get 200 y
while (b-a)/2>10^(-12)
c=(a+b)/2;

if (x^3+x*a+a^3)*(x^3+x*c+c^3)== 0 , %c is a solution, done
break;
end
if (x^3+x*a+a^3)*(x^3+x*c+c^3)<0 %a and c make the new interval
b=c;
else %c and b make the new interval
a=c;
end
end

y(i)=(a+b)/2; %new midpoint is best estimate
i=i+1;
disp('x=');
disp (x);

disp('y=');
disp (y) ;
end
x=0:0.01:2
plot(x,y)