Recent changes to 561: matrix(...taylor...)^^-1 wronghttp://sourceforge.net/p/maxima/bugs/561/Recent changes to 561: matrix(...taylor...)^^-1 wrongenWed, 21 Apr 2004 02:37:57 -0000matrix(...taylor...)^^-1 wronghttp://sourceforge.net/p/maxima/bugs/561/Consider
m: matrix\(\[taylor\(1+a\*x,x,0,1\),0\],
\[0,taylor\(1+d\*x,x,0,1\)\]\);
Now,
m^^-1
=> matrix\(\[1/\(d\*x+1\),0\],\[0,1/\(d\*x+1\)\]\)
There are two problems with this. First, the answer is
incorrect. Compare:
matrixmap\(ratdisrep,m\)^^-1
=> matrix\(\[1/\(a\*x+1\),0\],\[0,1/\(d\*x+1\)\]\)
Second, the answer is not in terms of taylor series. A
silent ratdisrep was done in the middle, losing truncation
information.
Stavros MacrakisWed, 21 Apr 2004 02:37:57 -0000http://sourceforge.net10cf800cac6e461a27d232c2b43a80d47cfdb6c9