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

Lösung

Matlab Code

function Bsp4( )
k = 1.06315;
pot = (1.5 * k) / (1 + k);


f = @(p)(p^(-1.5) * ((p-(1/(1+k))) / (1-(1/(1+k))))^pot - 0.15 );
a = 1;
b = 2;
maxerr = 0.0001;

c = (a*f(b) - b*f(a))/(f(b) - f(a));

while abs(f(c)) > maxerr
    a = b;
    b = c;
    c = (a*f(b) - b*f(a))/(f(b) - f(a));
end

    disp(c);
end