... 1f(x)2f(x)(2)Taxâydựnghàmhalley1()đểthựchiệnthuậttoántrên 241CHƯƠNG 5: CÁC PHƯƠNG TRÌNH PHI TUYẾN §1.KHÁINIỆMCHUNGNếuphương trình đạisốhaysiêuviệtkháphứctạpthìítkhitìmđượcnghiệmđúng.Bởivậyviệctìmnghiệmgầnđúngvàướclượngsaisốlàrấtcầnthiết. ... f0=x0‐(x1‐x0)*(x1‐x0)/(x2‐2.*x1+x0); x0=f0;endy=f0;Đểtìm nghiệmcủaphương trình x=(2‐ex+x2)/3=g(x)tadùng chương trình ctaitstef.mclearall,clcf=inline(ʹ(2.‐exp(x)+x.^2)/3ʹ);[x,y]=aitstef(f,0.,1e‐4,50)§9.PHƯƠNGPHÁPMUELLERTrongphươngphápdâycungkhitìmnghiệmtrongđoạn[a,b]taxấpxỉ ... 259P=Q;Y=feval(f,P);endifabs(p‐P(3))<abs(p‐P(2))R=[P(1)P(3)P(2)];P=R;Y=feval(f,P);end%cactritinhlanmoiP(3)=p;Y(3)=feval(f,P(3));y=Y(3);%dieukiendunglaperr=abs(z);relerr=err/(abs(p)+delta);if(err<delta)|(relerr<delta)|(abs(y)<eps)breakendendĐểgiảiphương trình sin(x)‐0.5*x=0tadùng chương trình ctmuller.mclearall,clcformatlongf=inline(ʹsin(x)‐0.5*xʹ);x=muller(f,1.8,2.2,50)§10.PHƯƠNGPHÁPHALLEYPhéplặpNewtontìmnghiệmcủahàmphương trình x=g(x)là:′f(x)g(x)=x‐f(x)(1)Tốcđộhộitụtăngđángkểkhihàmcónghiệm...