<div class="gmail_quote">On Mon, Mar 16, 2009 at 9:50 PM, Henning Thielemann <span dir="ltr"><<a href="mailto:lemming@henning-thielemann.de">lemming@henning-thielemann.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im">Polynomial arithmetic is also contained in<br></div>
<a href="http://hackage.haskell.org/packages/archive/numeric-prelude/0.0.5/doc/html/MathObj-Polynomial.html" target="_blank">http://hackage.haskell.org/packages/archive/numeric-prelude/0.0.5/doc/html/MathObj-Polynomial.html</a><br>
<a href="http://darcs.haskell.org/htam/src/Polynomial.hs" target="_blank">http://darcs.haskell.org/htam/src/Polynomial.hs</a></blockquote><div><br></div><div>Nice. But the above code doen't seem to contain functions to find (or approximate) all complex roots, which is the function I need (maybe with special cases for linear, quadratic, cubic, since these can be done analytically). The DSP package does, albeit only<span class="Apple-style-span" style="font-size: small;"> </span><span class="Apple-style-span" style="font-family: '-webkit-sans-serif'; line-height: 34px; "><span class="Apple-style-span" style="font-size: small;">Laguerre's method, which might be good enough, although my old C# code that I wanted to port used<span class="Apple-style-span" style="font-family: arial, helvetica, sans-serif;"><span class="Apple-style-span" style="font-size: small;"> </span></span><span class="Apple-style-span" style="line-height: normal; "><span class="Apple-style-span" style="font-family: arial, helvetica, sans-serif;"><span class="Apple-style-span" style="font-size: small;">Weierstrass method</span></span><span class="Apple-style-span" style="font-family: -webkit-sans-serif; font-size: 13px; font-weight: normal; line-height: 34px; "></span></span></span></span></div>
<div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><a href="http://darcs.haskell.org/htam/src/Polynomial.hs" target="_blank"></a>vector-space seems to be a tough dependency, since it relies on type families. NumericPrelude's VectorSpace class is a multi-parameter type class.</blockquote>
<div><br></div><div>NumericPrelude is very impressive, but it scares me a little, since it is so big :-)</div><div><br></div></div>