%initialgjetning
x=[1;1];

%denne koden er basert paa at vi overskriver en enkelt variabel istedet for aa lagre hele rekken av iterasjoner. derfor trenger vi en ekstra variabel for aa holde styr paa avstanden mellom iterasjonene.
y=x-[100;100];

%hvor tett skal det vaere mellom iterasjonene foer vi er fornoyde?
tol=10^(-16);

%fikspunkt for x-cos x = 0
while norm(y-x) > tol % vi holder paa saa lenge avstanden mellom iterasjonene er storre en tol
    y=x; % y er naa forrige iterasjon

    A=[4*x(2),-2*x(2);-x(1)/2, 2*x(1)]; %matrisen
    h=[x(1)^2+x(2)^2-1;x(1)^2/4+2*x(2)^2-1]; %hoyresiden
    
    x=x-A*h/(7*x(1)*x(2)); %newton

    pause(.5) %en liten pause, slik at vi skal faa tid til aa nyte synet av den nye iterasjonen
    disp([x(1) x(2) norm(x-y)]) %som printes her
end