Path: news.mathworks.com!not-for-mail
From: "prabhakaran m" <prabha.gahon@gmail.com>
Newsgroups: comp.soft-sys.matlab
Subject: Re: Newton's Method
Date: Fri, 20 Apr 2012 05:16:19 +0000 (UTC)
Organization: The MathWorks, Inc.
Lines: 25
Message-ID: <jmqrf2$e61$1@newscl01ah.mathworks.com>
References: <jmqiva$bci$1@newscl01ah.mathworks.com> <jmqlt6$mg1$1@newscl01ah.mathworks.com>
Reply-To: "prabhakaran m" <prabha.gahon@gmail.com>
NNTP-Posting-Host: www-04-blr.mathworks.com
Content-Type: text/plain; charset=UTF-8; format=flowed
Content-Transfer-Encoding: 8bit
X-Trace: newscl01ah.mathworks.com 1334898979 14529 172.30.248.35 (20 Apr 2012 05:16:19 GMT)
X-Complaints-To: news@mathworks.com
NNTP-Posting-Date: Fri, 20 Apr 2012 05:16:19 +0000 (UTC)
X-Newsreader: MATLAB Central Newsreader 2438980
Xref: news.mathworks.com comp.soft-sys.matlab:765276

Roger, i am also doing the same newton method. here code is not converging. what may the problem. also i need explanation when to use error greater or lesser than tolerance.

noiter=10
tolerence=1e-6;
del_int=[x0; y0; zta; rho; Rint]; % initial parameters
del=[];
del(1,:)=del_int;
for j=2:noiter
for i=n
%here in the loop i have found the H and J matrix. i have excluded its calcuations.
end
Jinv=inv(J);
tempo=del(j-1,:);
del(j,:)=tempo'-Jinv*H
error=max(abs(del(j,:)-del(j-1,:)))
if error <= tolerence
    final=j;
    break;
end
end
final_val=del(final,:)

when i evaluate the error, it goes on increasing. what may be the problem. also if reduce the tolerance value to 1e-2 it converges in the 2 iteration itself.
here is the single page reference form where i took the procedure
http://www.sendspace.com/file/x19u7q