Sehr geehrter Herr Prof. Frauendorfer
Ihr Problem hat eine einfache Loesung.
Modifizieren sie Ihre Differentialgleichungen so,
dass die Abnahme einer Konzentration bei negativem Wert gestopt wird.
Dies kann mithilfe eines If[]-Befehle geschehen.
Zu Ihrem Beispiel:
vorher:
eqn ={ a'[t]== -k1*(a[t])^n1, a[0]==a0,
b'[t]== +k1*(a[t])^n1
-k2*(b[t])^n2, b[0]==b0,
c'[t]== k2*(b[t])^n2, c[0]==c0};
nachher:
eqn ={ a'[t]== If[a[t] > 0, -k1*(a[t])^n1, 0], a[0]==a0,
b'[t]== +k1*(a[t])^n1
-k2*(b[t])^n2, b[0]==b0,
c'[t]== k2*(b[t])^n2, c[0]==c0};
MfG, Markus van Almsick
|