matlab二分法求零点

matlab二分法求零点 >> clear >> f=inline('x^2+exp(-2*x)-2'); >> fplot(f,[0,2]);grid; >> hold on; >> a=1;b=2;dlt=1.0e-4;k=1; >> while abs(b-a)>dlt c=(a+b)/2; if(c)==0 break; elseif f(c)*f(b)<0 a=c; else b=c; end plot(c,0,'*') c k=k+1; end c = 1.5000 c = 1.2500 c = 1.3750 c = 1.4375 c = 1.4063 c = 1.3906 c = 1.3984 c = 1.3945 c = 1.3926 c = 1.3916 c = 1.3921 c = 1.3923 c = 1.3922 c = 1.3922 >>>> clear f=inline('x^2+exp(-2*x)-2'); fplot(f,[-2,0]);grid; >> hold on; a=-1;b=0;dlt=1.0e-4;k=1; while abs(b-a)>dlt c=(a+b)/2; if(c)==0 break; elseif f(c)*f(b)<0 a=c; else b=c; end plot(c,0,'*') c k=k+1; end c = -0.5000 c = -0.2500 c = -0.3750 c = -0.3125 c = -0.3438 c = -0.3281 c = -0.3203 c = -0.3164 c = -0.3184 c = -0.3193 c = -0.3198 c = -0.3201 c = -0.3202 c = -0.3203