function fiksgraf
% FIKSGRAF: Illustrerer fikspunktiterasjoner paa et intervall [0, 1].
%           FIKSGRAF bruker en funksjon y = g(x). 

a = 0; b = 1;      % Intervallet [0, 1]

% Plotter y=x og y=g(x)
h = (b-a)/20;
t = [a:h:b];
y = g(t);
plot(t,y,t,t)
axis([a,b,a,b]);
axis('square');
hold on

% Fikspunktiterasjonene starter med initialverdi x0:
x0 = input('Startverdi x0 = ');
x1 = g(x0)
plot([x0,x0],[0,x1],'--r');
pause
plot([x0,x1],[x1,x1],'--r');
pause
x0 = x1;
for n = 2:10
    x1 = g(x0)
    plot([x0,x0],[x0,x1],'--r');
    pause
    plot([x0,x1],[x1,x1],'--r');
    pause
    x0 = x1;
end
hold off



