@BlackJack:

Vielen herzlichen Dank!!!
Es lag daran, dass dieses eine b mit dem anderen nichts zu tun hatte!
Das Unterprogramm läuft nun, auch wenn die Werte nicht ganz stimmig sind.

Danke!

So sieht der gesamte Quelltext des Unterprogs aus:

from scipy.linalg import *
from scipy.linalg import gmres
from numpy import *


"Funktionen F"
def g1(x, x_1, x_2, h, a, d):
c = x[0]*((3./2.)-h+h*x[0]) + ((h*a*x[0]*x[1])/(x[0]+d)) - 2.*x_1[0] + \
0.5*x_2[0]
print "g1=", c
return c

def g2(x ...

Ja genau, bei der ersten Berechnung werden feste Werte zugewiesen
und bei den folgenden nimmt er den neuberechnetet Wert für die Funktionsauswertung.

Hier die Problem-Funktionen:



def g2u(x, h, b):
g = (((x[1]**2)*h*b*-1.)/(x[0]**2))
print "g2u=", g
return g
def g2v(x, h, b):
l = ((3./2.)-h ...

Ich habe mir zwischenwerte angeben lassen und habe jetzt den Punkt gefunden welches das Problem darstellt.
Das Prog liefert:

sol1= [ 0.5 0.5]
g1u= 1.90816326531
g1v= 1.42857142857
g2u= -0.1
g2v= 1.6
[[ 1.90816327 1.42857143]
[-0.1 1.6 ]]
g1= 0.855242038286
g2= 0.372650654
b= [ 0.85524204 0 ...

Ja ich habe numpy importiert.

Es ist nur sehr verwirrend, da das Prog mit anderen Funktionen ohne Probleme läuft.

Vielleicht erläuter ich mal genauer das Problem:

also
g1u(x, h, a, d), g1v(x, h, a, d) , g2u(x, h, b) , g2v(x, h, b)
sind Funktionen, die als Input a,b,d,h reelle Zahlen Floats haben ...

Hallo zusammen!
Ich bearbeite gerade ein projekt in Python.
Unter anderem beinhaltet es ein bdf-verfahren mit newton iterationen.

Da ich seit kurzem erst mit python programmiere, habe ich probleme;
eines davon ist

Q = array([[g1u(x, h, a, d),g1v(x, h, a, d)],[ g2u(x, h, b),g2v(x, h, b ...