NB. Bring on the heavy weights! Find one root. NB. The other could be found by polynomial division NB. which I'm quite sure is demonstrated in the j phrases.
load'~addons/math/misc/amoeba.ijs' f=: _5 2 700&p. *:@:f amoeba (<30) 3,:0 NB. returns result and error estimate +---------+----------+ |0.0830989|1.8977e_14| +---------+----------+ f 0.0830989 _3.17315e_6 1e_6 1e6 1e_6 *:@:p. amoeba (<40) 1 ,: 4 +------------+-----------+ |_1.81899e_12|6.70744e_13| +------------+-----------+ f=: 1 0 1&p. NB. complex roots *:@:f amoeba (<30) 3,:1j1 |domain error | minindx=.tryerr i.minval=. <./tryerr=.uverb simptry ---------------------------------------------------------------------- For information about J forums see http://www.jsoftware.com/forums.htm
