"matrix([matrix([" returned with mat_function (Was: Re: Symbolic matrix power)

This is a repost, so it won't get lost.

On 10/25/2013 08:26 PM, andre maute wrote:
> On 10/25/2013 07:50 PM, Aleksas Domarkas wrote:
>> (%i1) A:matrix([2, 0, 0], [0, 2, 0], [-1, 0, 3])$
>> (%i2) load("diag")$
>> (%i3) integer_pow(x) := block ([k], declare (k, integer), x^k)$
>> (%i4) mat_function (integer_pow, A);
>> (%o4) matrix([2^k,0,0],[0,2^k,0],[2^k-3^k,0,3^k])
> Hi list,
> I tried the following
>
> [user at host ~]$ rmaxima
> Maxima 5.30.0 http://maxima.sourceforge.net
> using Lisp SBCL 1.1.2-1.fc18
> Distributed under the GNU Public License. See the file COPYING.
> Dedicated to the memory of William Schelter.
> The function bug_report() provides bug reporting information.
> STYLE-WARNING: redefining MAXIMA::$FILE_TYPE in DEFUN
> (%i1) display2d : false;
>
> (%o1) false
> (%i2) A : matrix([a,b],[c,d]);
>
> (%o2) matrix([a,b],[c,d])
> (%i3) load("diag");
>
> (%o3) "/home/user/opt/maxima/share/maxima/5.30.0/share/contrib/diag.mac"
> (%i4) integer_pow(x) := block ([k], declare (k, integer), x^k);
>
> (%o4) integer_pow(x):=block([k],declare(k,integer),x^k)
> (%i5) mat_function (integer_pow, A);
>
> (%o5) matrix([matrix([(sqrt(d^2-2*a*d+4*b*c+a^2)-d-a)^k* ...
>
> Where %o5 consists of a long long expression
> But can you see it the two occurrences of "matrix"?
>
> Andre
>
> _______________________________________________
> Maxima mailing list
> Maxima at math.utexas.edu
> http://www.math.utexas.edu/mailman/listinfo/maxima
>