ASKSAGE: Sage Q&A Forum - Latest question feedhttp://ask.sagemath.org/questions/Q&A Forum for SageenCopyright Sage, 2010. Some rights reserved under creative commons license.Mon, 25 Feb 2013 03:08:03 -0600Newton's identities in Sagehttp://ask.sagemath.org/question/9847/newtons-identities-in-sage/**EDIT**
I actually need:
$s_k=[c_1s_{k-1}+...+c_{k-1}s_1-kc_k]$ ? could somebody help me to change tobias welch's answer so that it computes $s_k$ instead of $c_k$?
**END EDIT**
I'm combining netwon's identities with le verrier's algorithm
I need some help coding the following on python.
$c_k=\frac{-1}{k}(s_k+c_1s_{k-1}+c_2s_{k-2}+\dots+c_{k-1}s_1)$
where $s_k=Tr(A^k)$, for some square matrix A, $\forall k=1,2,3,\dots,n$
So, i'd like to type in $c(k)$ and python spits out the value for $c_k$ as defined above.
jtaaMon, 25 Feb 2013 03:08:03 -0600http://ask.sagemath.org/question/9847/