partielle differentialgleichungen

[conharpython]

Hallo liebe Pythonnutzer,

ich bin grade dabei zu lernen, wie Differentialgleichungen in Phython gelöst werden. Das habe ich auch soweit verstanden, glaub ich:

Code: Alles auswählen

def model(y,t,k):
    dydt = - y
    return dydt

# initial condition
y0 = 5

# time points
t = np.linspace(0,20)

# solve ODEs
y1 = odeint(model,y0,t)

Was nicht noch nicht verstanden habe ist, wie ich partielle Differentialgleichungen in Python implementiere. Zum Beispiel:

du[x,t]/dt=u[x,t]+v[x,t]
dv[x,t]/dx=u[x,t]+v[x,t]

Mit u[0,t]=1 und v[x,0]==2 oder so.

Danke schön

[ThomasL]

Code: Alles auswählen

def model(y,t,k):
    dydt = - y
    return dydt

Diese Funktion gibt also nur den negierten Parameter y zurück.
Warum übergibt man dann noch t und k ?