Dovahkin schreef:Tomyris schreef:Dovahkin schreef:
jup, 3 dimensionaal dus dan heb je drie vergelijkingen
Wat moet je doen? Ik was niet het beste in dat vak maar misschien kan ik je toch nog een beetje helpen?
heb je ook verstand van mathematica?
ik heb dit voedselweb model:
eqgrass = D[N1[t], t] == r*N1[t]*(1 - (N1[t]/k)) - f*N1[t]*N2[t]
eqres = D[N2[t], t] == N2[t]*(b*N1[t] - d) - phi*N2[t]*N3[t]
eqherb = D[N3[t], t] == N3[t]*(beta*N2[t] - delta)
waarbij N1, N2 en N3 populatiegroottes zijn, k de carrying capacity van N2 en de rest zijn constanten
nu moet ik met een plot laten zien hoe de evenwichten veranderen als je k varieert, maar mathematica zegt dat ie geen oplossing kan bedenken voor de vergelijkingen
heb dit tot nu toe; voor alle constanten en k heb ik al een waarde gedefinieerd, wilde eerst een 'gewoon' plot maken
eqgrass = D[N1[t], t] == r*N1[t]*(1 - (N1[t]/k)) - f*N1[t]*N2[t]
eqres = D[N2[t], t] == N2[t]*(b*N1[t] - d) - phi*N2[t]*N3[t]
eqherb = D[N3[t], t] == N3[t]*(beta*N2[t] - delta)
sys2 = {eqgrass, eqres, eqherb}
initial = {N1[0] == 100, N2[0] == 200, N3[0] == 200}
soll = DSolve[{sys2, initial}, {N1[t], N2[t], N3[t, t]
Plot[soll, {t, 0, 50}]
maar ik krijg steeds ""DSolve:

svar: 0.0010214285714285716` cannot be used as a variable."" :/
Ik gebruik mathematica eigenlijk niet, heb python en C++ geleerd.