Forum / Numerische Methoden der Ingenieurwissenschaften / Lösung programmier Bsp für den 2.Test

Bin selber nach langen hin und her drauf gekommen (zumindest mal die Bisectionsmethode), wen`s noch interessiert

    iter = 0
    c = (a+b)/2 
   
    while((iter < nrIter) & (abs(f(c)) > tol) & (abs((b-a)/2) > tol)):
        iter = iter+1 
        c = (a+b)/2
        if(numpy.sign(f(c)) == numpy.sign(f(a))):
            a = c
        else:
           b = c        
        print 'Iter=', iter, 'Wert =', c 
                             
    
    return c,iter\n~~~

Eine Alternative.

Hat wer das Sekantenverfahren programmiert?

Hat wer das Sekantenverfahren programmiert?

       N       (max. Anzahl von Iterationen)

i = 1

Solange i <= N, wiederhole:
    Falls |f(x1)| < EPS
        Beende Schleife

    x = x - f(x1)*(x1 - x0)/(f(x1)-f(x0));

    Falls |x-x1| < |x1| * EPS
        Beende Schleife

    i = i + 1
    x0 = x1
    x1 = x

Gib x aus
Gib f(x) aus\n~~~

Hier als Pseudocode. Im Grunde genau wie das Newtonverfahren nur wird ein zusätzlicher Startwert übergeben und die Formel ist statt Newton halt Sekantenformel.

Kommt einer dieser Themen fix bei der Prüfung zum programmieren?

Fix is nix ;) Runge Kutta könnte ich mir auch noch vorstellen, das kommt

Hat das Runge-Kutta-Verfahren vielleicht irgendjemand geschafft zu programmieren und könnte uns den Quellcode hier zu Verfügung stellen? Wäre überaus freundlich!