3582284059

Plik Edycja Format Widok Pomoc

NEWTON JEDNOWYMIAROWY:

fiiąrneHontxk, fxk,Liczbattęracji3-ńewtonCf,fprim, x0, tol, Maxlt) Liczbalteracj m 0; xk m x0; fxk = feval(f,xk); fprimxk = feval(fprim, xk); while (abs(fxk) |> tol) & (Liezbalteracji < Maxlt)

Liczbalteracji m Liczbalteracji + 1; xk = xk - fxic/fprimxk; fxk « feval(f,xk); fprimxk = feval (fprim, xk);

end

|NEWTON JEDNOWYMIAROWY - skrypt testujący:

f=i ni i neC * x. A3-x+2')

fprim=inlineC^,3*x.A2«l’)

x0=-2

|tol=10A-16

IMaxlt=1000

newtonCf,fprim, x0, tol, Maxlt)

NEWTON WIELOWYMIAROWY:

funct iMaxIt) kontynuować! s= 1; kontynuować2 = 1

i on [xk,fxk, Liczbalteracj i ]=newton_nwyni(f_f newton_jiwym_jakobian, x0, tol, ) Liczbalteracji « 0; xk = x0; f_rozmiar = size(f, 1);

while (Liczbalteracji < Maxlt) && (Ckontynuowacl » X) fi (kontynuowac2 * 1»

jakobianWartosc = feval(newton_nwynL_jakobian, xk);

f xk = feval Cf _t xk);

Liczbalteracji = Liczbalteracji + 1; dxk =■ jakobi anWartosc \ (~fxkj; xk - xk 4- dxk;

if max(abs(dxk)) < tol kontynuować! = 0; end    T

if max(abs(fxk)) < tol ^x kontynuowaćz = 0;

end