bisection methods?

[木心]

谢谢!
我觉得现在的问题是怎么把 y 和 x 一一对应??

a=0;b=2;


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

if (x^3+x*c+c^3)*(x^3+x*c+c^3)== 0 , %c is a solution, done
break;
end
if (x^3+x*c+c^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

d=(a+b)/2; %new midpoint is best estimate

y=d;
disp(x);
disp('y=');
disp(y);
end

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