% Dette skriptet løser et ikke-lineært ligningssystem med Newtons
% metode.
%
% Bruker eksempelet:
%
%   2x - cos(y) = 0
%   2y - sin(x) = 0
%

% Initial-verdier:
x = 0.7; y = 0.7;
% Antall iterasjoner:
MaxIter = 5;

for n=1:MaxIter
    F = [2*x - cos(y);        % NB: F er en kolonnevektor
	 2*y - sin(x)];      
    J = [2,       sin(y);     % Jacobimatrisa
	 -cos(x), 2      ];
    dx = -J\F                 % Løser lign.systemet J*dx = -F
    x = x + dx(1);            % Oppdaterer løsningen
    y = y + dx(2);
end
x, y