不同方法求方程的根matlab代碼

不同方法求方程的根matlab

代碼(總4頁)-CAL-FENGHAI.-(YICAI)-CompanyOne1-CAL-本頁僅作為文檔封面，使用請直接刪除#1、%使用二分法找到方程600x八4-550xA3+200x八2-20x-1=0在區(qū)間［,1］上的根，%誤差限為e=10八-4disp('二分法')a=;b=;tol=;n0=10;fa=600*(a.A4)-550*(a.A3)+200*(a.A2)-20*a-1;fori=1:n0p=(a+b)/2;fp=600*(p.A4)-550*(p.A3)+200*(p.A2)-20*p-1;iffp==0||(abs((b-a)/2)<tol)disp(用二分法求得方程的根p=')disp(p)disp('二分迭代次數(shù)為：’)disp(i)break;endiffa*fp>0a=p;elseb=p;endendifi==n0&&~(fp==0||(abs((b-a)/2)<tol))disp(n0)disp('次二分迭代后沒有求出方程的根?)end2、%使用牛頓法找到方程600xA4-550xA3+200xA2-20x-1=0在區(qū)間［,1］上的根，%誤差限為e=10A-4disp('牛頓法’)p0=;fori=1:n0p=p0-(600*(p0.A4)-550*(p0.A3)+200*(p0.A2)-20*p0-1)./(2400*(p0.A3)-1650*p0.A2+400*p0-20);if(abs(p-p0)<tol)dispd用牛頓法求得方程的根p=')disp(p)disp(‘牛頓迭代次數(shù)為：’)disp(i)break;endp0=p;endifi==n0&&~(abs(p-p0)<tol)disp(n0)disp('次牛頓迭代后沒有求出方程的根?)end3、%使用割線法找到方程600x八4-550x八3+200x八2-20x-1=0在區(qū)間［,1］上的根，%誤差限為e=10八-4disp('割線法’)p0=;p1=;q0=600*(p0.八4)-550*(p0.八3)+200*(p0.八2)-20*p0-1;q1=600*(p1.A4)-550*(p1.A3)+200*(p1.A2)-20*p1-1;fori=2:n0p=p1-q1*(p1-p0)/(q1-q0);ifabs(p-p1)<toldispd用割線法求得方程的根p=')disp(p)disp('割線法迭代次數(shù)為：’)disp(i)break;endp0=p1;q0=q1;pp=p1;p1=p;q1=600*(p.A4)-550*(p.A3)+200*(p.A2)-20*p-1;endifi==n0&&~(abs(p-pp)<tol)disp(n0)disp('次割線法迭代后沒有求出方程的根?)end4、%使用試位法找到方程600xA4-550xA3+200xA2-20x-1=0在區(qū)間［,1］上的根，%誤差限為e=10A-4disp('試位法')p0=;p1=;q0=600*(p0.A4)-550*(p0.A3)+200*(p0.A2)-20*p0-1;q1=600*(p1.A4)-550*(p1.A3)+200*(p1.A2)-20*p1-1;fori=2:n0p=p1-q1*(p1-p0)/(q1-q0);ifabs(p-p1)<toldisp('用試位法求得方程的根p=')disp(p)disp('試位法迭代次數(shù)為：’)disp(i)break;endq=600*(p.A4)-550*(p.A3)+200*(p.A2)-20*p-1;ifq*q1<0p0=p1;q0=q1;endpp=p1;p1=p;q1=q;endifi==n0&&~(abs(p-pp)<tol)disp(n0)disp('次試位法迭代后沒有求出方程的根?)end5、%使用muller方法找到方程600乂八4-550乂八3+200乂八2-20x-1=0在區(qū)間［,1］上的根，%誤差限為e=10A-4disp('muller法')x0=;x1=;x2=;h1=x1-x0;h2=x2-x1;d1=((600*(x1.A4)-550*(x1.A3)+200*(x1.A2)-20*x1-1)-(600*(x0.A4)-550*(x0.A3)+200*(x0.A2)-20*x0-1))/h1;d2=((600*(x2.A4)-550*(x2.A3)+200*(x2.A2)-20*x2-1)-(600*(x1.A4)-550*(x1.A3)+200*(x1.A2)-20*x1-1))/h2;d=(d2-d1)/(h2+h1);fori=3:n0b=d2+h2*d;D=(b*b-4*(600*(x2.A4)-550*(x2.A3)+200*(x2.A2)-20*x2-1)*d)A;if(abs(d-D)<abs(d+D))E=b+D;elseE=b-D;endh=-2*(600*(x2.A4)-550*(x2.A3)+200*(x2.A2)-20*x2-1)/E;p=x2+h;ifabs(h)<toldispd用muller方法求得方程的根p=')disp(p)disp（'muller方法迭代次數(shù)為：’）disp（i）break;endx0=x1;x1=x2;x2=p;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-