steffenson_solver

[alb]

je cherche les solutions de x^(1/4)-(ln(x))^2=0
les methodes où l'on donne l'intervalle fonctionnent
en principe la doc dit que pour steffenson_solver on peut se contenter d'une valeur de depart
Or fsolve(x^(1/4)-(ln(x))^2,x,1,steffenson_solver) renvoie Valeur Argument Incorrecte. Pourquoi ?

Pour info:
fsolve(x^(1/4)-(ln(x))^2,x=0..1,bisection_solver) renvoie [0.408917902095]
fsolve(x^(1/4)-(ln(x))^2,x=3..4,bisection_solver) renvoie [3.17520064334]
fsolve(x^(1/4)-(ln(x))^2,x=10^11..10^12,bisection_solver) renvoie [214910065296]

[parisse]

Un petit bug, la valeur initiale n'etait pas passee correctement. Par contre je n'arrive pas a empecher l'exception generee par la Gnu Scientific Library si la fonction/derivee n'est pas definie reelle pour l'estimation de depart.