Page 1 sur 1

solve bizarre

Publié : ven. janv. 30, 2009 9:44 am
par Guillaume
Bonjour,

Code : Tout sélectionner

solve(diff(1/4*exp(-x/2)*(x^2+2*x),x),x)
fonctionne très bien mais

Code : Tout sélectionner

solve(diff(0.25*exp(-x/2)*(x^2+2*x),x),x)
plante la session.

Si on précise le point de départ

Code : Tout sélectionner

fsolve(diff(0.25*exp(-x/2)*(x^2+2*x),x),x,2,newton_solver)
on obtient bien une réponse. Apparemment, fsolve prend par défaut xo=0 et donc ça plante car f''(0)=0.
Est-ce qu'on peut imaginer un test de départ pour le xo de fsolve ?

Publié : ven. janv. 30, 2009 11:49 am
par parisse
Je n'ai pas de plantage avec la derniere version, j'obtiens pour
solve(diff(0.25*exp(-x/2)*(x^2+2*x),x),x)
[3.2360679775+4.72048977526e-10*i]
Peut-etre qu'une modif recente a corrige ce bug!

Publié : ven. janv. 30, 2009 5:14 pm
par Guillaume
Effectivement... J'ai dû loupé une maj.
Merci !