[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
roots(poly(r)) != r
From: |
Torben Mikael Hansen |
Subject: |
roots(poly(r)) != r |
Date: |
Fri, 2 Mar 2007 23:38:02 +0100 |
User-agent: |
KMail/1.9.5 |
Dear All
I have been playing around with poly and roots, don't understands this
a=rand(100,1);
r=exp(I*.1*pi*2*[a;-a]);
e=max(sort(abs(roots(poly(r))))-sort(abs(r)));
Shouldn't this give an "e" that close to zero. I always get something in the
order of 10.
Thanks
Torben
- roots(poly(r)) != r,
Torben Mikael Hansen <=