solve bizarre

Guillaume · Message par **Guillaume** » ven. janv. 30, 2009 9:44 am

Bonjour,

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 ?

Message par **parisse** » ven. janv. 30, 2009 11:49 am

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!

Guillaume · Message par **Guillaume** » ven. janv. 30, 2009 5:14 pm

Effectivement... J'ai dû loupé une maj.
Merci !