J'essaie de résoudre le système d'équations différentielles donné par les variables eqn1 et eqn2.Matlab: dsolve avec des conditions ne fonctionnant pas
lambda1 = 3;
lambda2 = 2;
gamma1 = 1;
gamma2 = 1;
delta1 = 1;
delta2 = 1;
syms n1(t) n2(t)
eqn1 = diff(n1) == (lambda1 - gamma1)*n1 - delta1*(n1 + n2)*n1;
eqn2 = diff(n2) == (lambda2 - gamma2)*n2 - delta1*(n1 + n2)*n2;
c1 = n1(0) == 10;
c2 = n2(0) == 10;
[a, b] = dsolve(eqn1, eqn2, c1, c2)
Si je ne précise pas les conditions c1 et c2 et exécutez simplement:
[a, b] = dsolve(eqn1, eqn2)
alors tout va bien et je reçois:
a =
(6*exp(2*t)*n1(0))/(3*exp(2*t)*n1(0) + 4*exp((3*t)/2)*n2(0) + 6)
b =
(6*exp(t/2)*exp(t)*n2(0))/(3*exp(2*t)*n1(0) + 4*exp((3*t)/2)*n2(0) + 6)
mais dès que je tente pour spécifier n1 (0) et n2 (0), comme vous pouvez le voir dans le code source, j'obtiens l'avertissement suivant:
a =
Warning: The result cannot be displayed due a previously interrupted computation or out of memory. Run 'reset(symengine)' and rerun the commands to regenerate the result.
> In sym.disp at 36
In sym.display at 37
In competingForResources at 15
b =
[]
Et je suis incapable de tracer le résultat. Qu'est-ce qui me manque là? Comment spécifier les conditions de la bonne manière, donc j'obtiens un résultat plottable?