j'ai un systeme equation diff du type suivant:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
dx(s_i,t)/dt=f(x(s_i,t), x(s_i,t)\sum_i x(s_i,t), y(s_i,t), y(s_i,t)\sum_i y(s_i,t), z(s_i,t), z(s_i,t)\sum_i z(s_i,t))
dy(s_i,t)/dt=g(x(s_i,t), x(s_i,t)\sum_i x(s_i,t), y(s_i,t), y(s_i,t)\sum_i y(s_i,t), z(s_i,t), z(s_i,t)\sum_i z(s_i,t))
dz(s_i,t)/dt=h(x(s_i,t), x(s_i,t)\sum_i x(s_i,t), y(s_i,t), y(s_i,t)\sum_i y(s_i,t), z(s_i,t), z(s_i,t)\sum_i z(s_i,t))
avec x, y, z: une fonction de s_i et de t et i=1:50 (s_i c'est une discretisation d'une deuxiemme variables).
f,g,h: des fonctions de type polynome
je peux vous envoyez les fonctions si cela est important

Je sais demontre mathematiquent que si la condition initiale dans [0,1] alors la solution reste dans [0,1].

Malheureusement qt je simule les solutions et pour certaine valeurs des parametres (intervalle asses petit), l'une des composante de la solution qui devrai tendre vers 0 commence par etre negative du type -0.8e-15 et ensuite cette composante devient clairement negative (~-1000).

Je me pose la question si c'est un probleme du solver ou un que j'ai une erreur dans l'ecriture de ma fonction

merci pour votre aide