《4bb[理学]试验三 非线性方程求解》由会员分享,可在线阅读,更多相关《4bb[理学]试验三 非线性方程求解(5页珍藏版)》请在金锄头文库上搜索。
1、试验三 非线性方程求解1、 二分法function x_star, index, it=bisect(fun, a, b, ep)% 求解非线性方程的二分法,其中% fun(x) - 需要求根的函数% a,b - 初始区间的端点。% ep - 精度要求,当 (b-a)/20x_star=fa, fb; index=0; it=0;return;endk=1;while abs(b-a)/2=epx=(a+b)/2; fx=feval(fun,x);if fx*fa= it_max) 。% it - 迭代次数。if nargin = it_max) 。% it - 迭代次数。if nargin
2、= it_max) 。% it - 迭代次数。if nargin =tol) %循环直到 totol 为止j=j+1;to=(abs(simpson_integral(x,y,2j)-simpson_integral(x,y,1)/10;%精度测试式ends=s+simpson_integral(x,y,2j);endfunction y=f(x)% f-被积函数y=x5+sin(x)*x6*sqrt(x);function s=simpson_integral(a,b,m)%input: a-下限% b-上限% m-步数h=(b-a)/(2*m);s1=0;s2=0;for i=1:(m-1)x=a+2*i*h; s1=s1+f(f,x);endfor i=1:mx=a+(2*i-1)*h;s2=s2+f(f,x);ends=h/3*(f(f,a)+f(f,b)+2*s1+4*s2);