% Forslag til loesning av oppgave 1b, oving 1.

% Oppgi g og dg.
g = inline('log(4+x-x^2)')
dg = inline('(1-2*x)/(4+x-x^2)')

r = 1.288677966823868; % Eksakt loesning, funnet fra 1a).

% Fikspunkt-iterasjonen
x = 1.5;
Nit = 20;
disp('Start fikspunktiterasjonene')
pause

for n=1:Nit
    x = g(x);
    disp(sprintf('n = %2d, x = %11.8f, |r-x| = %8.2e',n,x,abs(r-x)));
end

% Analyse
disp(' ');
disp('Finn et passende intervall for konvergensanalyse'); 
pause
disp('Plott dg(x)'); pause;

figure(1);
t = linspace(1,2,101);
fplot(dg,[1,2]); grid
xlabel('x'); ylabel('dg(x)');
axis('square');
pause

disp(' ');
disp('Plott g(x)'); pause
figure(2);
fplot(g,[1,2]);
hold on 
plot(t,t,'g');
xlabel('x'); ylabel('g(x)');
grid
axis('square');
hold off





    
    

