Résolution système d'équations différentielles

[Fanny Gallais]

Bonjour,

Je travaille sur un système de 2 équations différentielles modélisant la concentration d'un médicament dans deux compartiments (sang et tissu).
Ma condition initiale à t=0 est c(100,0). Je fais la résolution par la fonction lsoda du package deSolve.
Mon problème est que je voudrais modéliser une administration répétée de médicament. Autrement dit, je veux ajouter 100 non seulement à t=0 mais aussi à t=20, 40, ...
Pour cela j'ai essayé de rajouter un terme à ma première équation différentielle : il vaut 100 pour les temps d'intérêt, 0 sinon. Mais cela ne semble pas fonctionner.

Auriez-vous une idée pour régler ce problème ?

Merci

[Logez Maxime]

Bonjour,

est-ce qu'en rajoutant quelque chose comme ceci dans ta fonction ça peut t'aider :

Code : Tout sélectionner

function(x) ... +0^(x%%20)*100

Comme ça dès que tu obtiens un multiple de 20 tu rajoutes 100. Après je n'ai aucune expérience avec le package que tu utilises ni avec les résolutions d'équa diff.

Cordialement,
Maxime