exn517.m

来自「基天子matlab的双缝干涉从单缝(线光源)发出的单色光射到与之平行的双缝上」· M 代码 · 共 19 行

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第五章例5－1－7程序
% 用切线法求超越方程的近似解
% 西安电子科技大学出版社出版  陈怀琛编著 2007年5月
%
clear,clf
fplot('x.^3+10*x.^2-2*sin(x)-50',[-12, 5])
hold on,grid,e=1;
x0=input(' x0=给出要求的解的近似坐标(书上分别为-10,-2,2) ');
x10=x0;plot(x0,0,'*'),shg, pause(2)		% 画出新点在x轴上的位置
while e>0.0001
   f= x0.^3+10*x0.^2-2*sin(x0)-50;	 	% 求f(x0)
   g=3*x0^2+20*x0-2*cos(x0); 			% 求函数的导数g(x0)
   xx=x0-f/g; 		% 切线法求解公式
   e=abs(xx-x0);x0=xx 		% 精度控制,把新值赋予x0
   plot(x0,0,'*'),shg, pause(2)		% 画出新点在x轴上的位置
end  
pause
x1=fzero('exn517f',x10)               % 直接调用求过零点函数的解法