TU Wien:Computernumerik VU (Schranz-Kirlinger)/Übungen SS16/Beispiel 22

Aus VoWi
Zur Navigation springen Zur Suche springen

Lösung[Bearbeiten | Quelltext bearbeiten]

Matlab Code

function Bsp2

x0 = 0.0001;

fprintf('x0: %.10f\n', x0);
xn1 = x0;
xn2 = x0;
xn3 = x0;
for i = 1:10
    xn1 = -log(xn1);
    xn2 = exp(-xn2);
    xn3 = (xn3 + exp(-xn3))/2;
    fprintf('x%d -- I1: %.10f, I2: %.10f, I3: %.10f \n',i, xn1, xn2, xn3);
end

fprintf('I1: %.10f, I2: %.10f, I3: %.10f \n', xn1 + log(xn1), xn2 + log(xn2), xn3 + log(xn3));

end