Changes http://wiki.axiom-developer.org/339InterpreterGivesUpInLoop/diff
--
The following is a very surprising for a novice, and a bit surprising for me:
\begin{axiom}
c 0==1; c 1==1; c n==c(n-1)+s*(reduce(+,[q^k*c(k)*c(n-1-k) for k in 1..n-1]));
c 5
\end{axiom}
Note that the type of 'c 5' is 'FRAC POLY INT'. This is, because Axiom is not
smart enough to realize that $q^k$ is a polynomial if $k$ is nonnegative. The
following is more surprising, though:
\begin{axiom}
coefficient(univariate(c(5),s),3)
\end{axiom}
Axiom manages - after a *lot* of trying- to coerce 'univariate(c(5),s)' to a
'SUP FRAC POLY INT'. Wow.
\begin{axiom}
[coefficient(univariate(c(n),s),1) for n in 0..4]
\end{axiom}
But why does it give up here?
Martin
--
forwarded from http://wiki.axiom-developer.org/[EMAIL PROTECTED]
