牛頓法和割線法

1、點，精度要求 10-10-10。分法fun c tio nX=bisection(f x ,x a, xb , n ,d elta )%二分法解方程 %f x就是由方程轉(zhuǎn)化得關(guān)于x得函數(shù)，有fx=0。xa解區(qū)間上限 xb解區(qū)間下限 解區(qū)間人為判斷輸入 n最多循環(huán)步數(shù)，防止死循環(huán)。%delta為允許誤差x= x a;fa=eva 1( fx );x = xb; fb =eval (fx);for i= l:n作業(yè)（第五章）：1 1、在區(qū)間 （0,10,1 、5 5）別用二分法、頓法與割線法編程求下面得函數(shù)得零xc =( X a + xb)/2 ; x= X c ; fc=e va l(fx)X

2、= i , x c ,fcxb = xc ;syins sfx=cosfx=cos(2+K)茁2:Ma=OMa=O：xb=l.xb=l. 5:5:tT=L)00tT=L)00；dede LtLt a=a= lOi-ift;lOi-ift;SS =bis=bis ctian.(fxjctian.(fxj xbjxbj n ndede ItIt a)a)34.34. 000000000.51490.5149-0.-0. COOOCOOOif fc興fa end頓法functionX=n ewt on(fx,e,x0,m) x=x0;k=0;F=eval(fx);ifabs(F)=eX=xF;dis

3、p（X）； eturnendwhile kmk;returnx=x0；g=evaldiff(fx);fprintf（/牛頓法迭代M次沒有找到方程得根/ ）returnX0=x1;endfprintf(/n%s% 4ft%s%dX= ,X，k=,k)%俞出結(jié)果牛頓法結(jié)果：迭代5次 結(jié)果0、5149ATLAB0 0 CmfriandCmfriandWindowWindow 割線法：function X=gx9(fx, x0,x1,m,e) x=x0;f0=eval(fx);syitLSsyitLSX Xfyfy：=cjs=cjs (2*1)2(2*1)2-宴2;2;e=10-10e=10-10：1

4、1=100011=1000：KK =evton=evton(fij(fij e,e, xO,xO,- -0.0. 514951490 05.00005.0000endx=x1;f1=eva1(fx);ifabs(f0)二eendforX=x0,f0:;ifabs(f0)abs(f1)b=xO;xO=x1;x1=b;=fO;fO=f1;f1=b;endt=x1xO)/(f1fO);x0=x1；fO=f1x1=x1t*fx=x1;eva1(fx);bs(1)0y0=y1; c o nti n ueen dX 0=x k;x1 =x;X =g X 9( f X, X 0, X 1, m ,e) ; %割線法i f xbX=x ； disp (X)；r etur ny 0=y l；end%x=xO ; f 0 =ev a l(f x);= x1 ；f 1=eval(f XxO ,f 0:;ndf o r k=2f unc tion X=gx9(fX， x0 ， x1,m ， e)ifa bs(fO ) = eifabs(f0)a bsf1 )b=X 0 ;x0 =x1 ;x1 =b;= f0 ;fO=f 1;1=b;t=( x1 x0 ) ( f1 f0 ) ;x0= X 1;f0=f 1 ;x1=x 1 t*1;x