@parisse : J'ai réutilisé Khicas aujourd'hui pour tester le solveur d'équa diff... et j'ai noté des erreurs flagrantes !
En fait, le solveur oublie des constantes d'intégration ou bien oublie de changer leur nom au fur et à mesure que de nouvelles constantes apparaissent, même dans des cas très simples.
Par exemple :
-
desolve(y''=y',x,y)
renvoie c_0 e^x
au lieu de c_0 e^x + c_1
-
desolve(y''=0,x,y)
renvoie x c_0 + c_0
au lieu de x c_0 + c_1
-
desolve(y'''=0,x,y)
renvoie (x^2 c_0 + 2x c_0 + c_0)/2
au lieu de (x^2 c_0 + 2x c_1 + c_2)/2
-
desolve(y'''-3y'+2y=0,x,y)
renvoie c_0 e^x
au lieu de (x c_0 +c_1) e^x + c_2 e^(-2*x)
On voit apparaître également très souvent des "Error : Bad Argument Value", par exemple avec
desolve(y''+y=0,x,y)
, l'erreur complète dans ce cas étant "Unable to divide, perhaps due to rounding error"
Enfin, dans des cas que je n'ai pas réussi à reproduire, j'ai obtenu une transformée inverse de Laplace... mais non calculée.
Bref, la fonction "desolve" semble avoir de gros ennuis...