弦截法求解方程通用程序matlab,王能超 计算方法——算法设计及MATLAB实现课后代码...

为衡量数值解的精度，我们求出该方程的解析解为y?e?x?x.在此也以文件的形式表示如下：

function y=solvef3(x) y=exp(-x)+x;

令f=@f3; a=0; b=1; N=10; ya=1; 运行：A2=Adams2PC(f,a,b,N,ya)； A4=Adams4PC(f,a,b,N,ya); CA4=CAdams4PC(f,a,b,N,ya); y=solvef3(a:(b-a)/N:b); m=[A2,A4(:,2),CA4(:.2),y’]

其中m共有5列数据：左→右依次为离散节点值、二阶Adams预报校正系统所求解、四阶Adams预报校正系统所求解、改进四阶Adams预报校正系统所求解和准确解,精度依次递增.

第四章 方程求根 4.1 二分法

用二分法求解非线性方程f(x)?0在区间[a,b]内的根. MATLAB文件：(文件名：demimethod.m) function [x,k]=demimethod(a,b,f,emg) %a,b表示求解区间[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&

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！