Matlab 用二分法求方程先x^3-3*x+2-e^x=0的实根,要求精度e<10^-8
Matlab用二分法求方程先x^3-3*x+2-e^x=0的实根,要求精度e<10^-8,谢谢各位高手了!请求完整的Matlab程序,完整~~~...
Matlab 用二分法求方程先x^3-3*x+2-e^x=0的实根,要求精度e<10^-8,谢谢各位高手了!
请求完整的Matlab程序,完整~~~ 展开
请求完整的Matlab程序,完整~~~ 展开
2个回答
展开全部
源代码:
在matlab中保存为:demimethod.m
function [x,k]=demimethod(a,b,f,emg)
% a,b :求解区间的两个端点
% f :所求方程的函数名
% emg :精度指标
% x:所求近似解
% k: 循环次数
fa=feval(f,a);
fab=feval(f,(a+b)/2);
k=0;
while abs(b-a)>emg
if fab==0
x=(a+b)/2;
return;
elseif fa*fab<0
b=(a+b)/2;
else
a=(a+b)/2;
end
fa=feval(f,a);
fab=feval(f,(a+b)/2);
k=k+1;
end
x=(a+b)/2;
命令窗口输入以下内容:
>> f=@(x)x^3-3*x+2-exp(x);
>> demimethod(-5,5,f,10^(-8))
ans =
-1.9846
在matlab中保存为:demimethod.m
function [x,k]=demimethod(a,b,f,emg)
% a,b :求解区间的两个端点
% f :所求方程的函数名
% emg :精度指标
% x:所求近似解
% k: 循环次数
fa=feval(f,a);
fab=feval(f,(a+b)/2);
k=0;
while abs(b-a)>emg
if fab==0
x=(a+b)/2;
return;
elseif fa*fab<0
b=(a+b)/2;
else
a=(a+b)/2;
end
fa=feval(f,a);
fab=feval(f,(a+b)/2);
k=k+1;
end
x=(a+b)/2;
命令窗口输入以下内容:
>> f=@(x)x^3-3*x+2-exp(x);
>> demimethod(-5,5,f,10^(-8))
ans =
-1.9846
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询