实验八.方程根的MATLAB求解

实验内容

实 验 步 骤、过 程 

1.(1)二分法

第一步：在 MATLAB 软件，建立一个实现二分法的 MATLAB 函数文件su1.m如下：

function x=su1(fname,a,b,e)
%fname 为函数名，a,b 为区间端点，e 为精度
fa=feval(fname,a); %把 a 端点代入函数，求 fa
fb=feval(fname,b); %把 b 端点代入函数，求 fb
if fa*fb>0 error('两端函数值为同号');
end
%如果 fa*fb>0，则输出两端函数值为同号
k=0
x=(a+b)/2
while(b-a)>(2*e) %循环条件的限制
fx=feval(fname,x);%把 x 代入代入函数，求 fx
if fa*fx<0%如果 fa 与 fx 同号，则把 x 赋给 b，把 fx 赋给 fb
b=x;
fb=fx;
else
%如果 fa 与 fx 异号&#