Forum / Numerische Methoden der Ingenieurwissenschaften / Lösung programmier Bsp für den 2.Test
Lösung programmier Bsp für den 2.Test
Hat schon irgendwer das Bisectionsverfahren, das Newton Cotes-Verfahren: Simpsonregel, oder das Runge Kutta Verfahren Programmiert, und könnte das auch hochladen?? Wäre sehr dankbar im Namen aller die sich in programmieren nicht so wirklich auskennen.
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?
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.
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!
Harald @Hari1Autor
Maschinenbau · Technische Universit...
Bin selber nach langen hin und her drauf gekommen (zumindest mal die Bisectionsmethode), wen`s noch interessiert