不同方法求方程的根matlab代码.doc

《不同方法求方程的根matlab代码.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、1、%使用二分法找到方程600x^4-550x^3+200x^2-20x-1=0在区间[0.1,1]上的根，%误差限为e=10^-4disp('二分法')a=0.2;b=0.26;tol=0.0001;n0=10;fa=600*(a.^4)-550*(a.^3)+200*(a.^2)-20*a-1;fori=1:n0p=(a+b)/2;fp=600*(p.^4)-550*(p.^3)+200*(p.^2)-20*p-1;iffp==0

2、

3、(abs((b-a)/2)

4、迭代次数为：')disp(i)break;endiffa*fp>0a=p;elseb=p;endendifi==n0&&~(fp==0

5、

6、(abs((b-a)/2)

7、0.^3)-1650*p0.^2+400*p0-20);if(abs(p-p0)

8、4)-550*(p0.^3)+200*(p0.^2)-20*p0-1;q1=600*(p1.^4)-550*(p1.^3)+200*(p1.^2)-20*p1-1;fori=2:n0p=p1-q1*(p1-p0)/(q1-q0);ifabs(p-p1)

9、p-pp)

10、

11、0^-4disp('muller法')x0=0.1;x1=0.2;x2=0.25;h1=x1-x0;h2=x2-x1;d1=((600*(x1.^4)-550*(x1.^3)+200*(x1.^2)-20*x1-1)-(600*(x0.^4)-550*(x0.^3)+200*(x0.^2)-20*x0-1))/h1;d2=((600*(x2.^4)-550*(x2.^3)+200*(x2.^2)-20*x2-1)-(600*(x1.^4)-550*(x1.^3)+200*(x1.^2)-20*x1-1))/h2;d=(d2-d1)/(h2+h

12、1);fori=3:n0b=d2+h2*d;D=(b*b-4*(600*(x2.^4)-550*(x2.^3)+200*(x2.^2)-20*x2-1)*d)^0.