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

Aus VoWi
Zur Navigation springen Zur Suche springen

Lösung[Bearbeiten | Quelltext bearbeiten]

Matlab Code

function Bsp1

x0 = 40;
fprintf('x0: %.10f\n',x0);
fprintf('ITERATIV:\n');
xn = x0;
for i = 1:15
    xn = 1/((xn^2) * 2); %T(xn+1) = T(xn)
    fprintf('x%d: %10.15f\n',i,xn);
end

xn = x0;
fprintf('NEWTON:\n');
for i = 1:20
    F = xn^3 - 1/2;
    Fx = 1/(3*xn^2);
    xn = xn - Fx * F;
    fprintf('x%d: %10.15f\n',i,xn);
end

end