Tach,
Die folgende Funktion stellt die Rosenbrock-Wanner Methode der Ordnung 2 dar:
def row_2_step(f, Jf, yi, h): yi = atleast_2d(yi) n = yi.shape[0] a = 1.0/(2.0 + sqrt(2.0)) I = identity(n) J = Jf(yi) A = I - a*h*J b1 = f(yi) k1 = solve(A, b1) b2 = f(yi+0.5*h*k1) - a*h*dot(J,k1) k2 = solve(A, b2) yip1 = yi + h*k2 return yip1
das ist nicht vollständig; welche Imports gehören noch zu der Methode, numpy?
... leider verstehe ich ein paar Python-Funktionen nicht bzw. wüsste keine 'Entsprechung' in Matlab.
Welche?
mfg
Woodfighter
P.S. Bitte weise auf Crosspostings (https://stackoverflow.com/questions/31170329/implementation-of-rosenbrock-wanner-method-python-vs-matlab, https://math.stackexchange.com/questions/1346162/python-vs-matlab) hin, damit Leute sich nicht, für bereits beantwortete Dinge, Arbeit machen.