求教牛顿法求根哪错了

[复制链接]

字体大小: 正常放大

silence123

7 主题	8 听众	317 积分

升级 5.67%

TA的每日心情

	慵懒 2016-10-11 09:39

签到天数: 173 天

[LV.7]常住居民III

自我介绍: 想学习数学建模

群组: 数学建模培训课堂2

群组: 2011年第一期数学建模

群组: 数学建模培训课堂1

电梯直达

1^#

发表于 2012-12-13 23:14 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

%牛顿法求解方程的根;
a=[];
a(1)=0.5;
e=0.0002;
syms x
f=x*exp(x)-1;
N=15;
z=diff(f);
fprintf(  'f的导数为%s,z')
for k=2:N
% f的导数为exp(x)+x*exp(x);
if exp(a(1))+a(1)*exp(a(1)) ==0
   disp(fprintf('奇异标志'));
else
   a(k)=a(k-1)-(a(k-1)*exp(a(k-1))-1)/(exp(a(k-1))+a(k-1)*exp(a(k-1)));
   if abs(a(k)-a(k-1))<=e;
         disp(fprintf('该方程的零点为d=%4.6f,a(k)'))
   else
         if k==N
            fprintf('奇异标志');
         end
   end
end
end

zan