Bisher habe ich das Problem in Matlab gelöst, eben über eine Nullstellensuche bzw. der Matlab-Funktion fzero().

Matlab nutzt Brent's Method, was die von Rolf genannten Methoden kombiniert. Was besseres fiele mir sponatan auch nicht ein.