You can redefine the sym_matrixp function as
sym_matrixp(expr) :=
if atom(expr) then featurep(expr, sym_matrix)
else if member(part(expr, 0), ["+", "-", "."]) then
every(sym_matrixp, args(expr))
else if part(expr, 0)="^^" then sym_matrixp(first(expr))
else false$
Andrej
On Mon, Mar 2, 2009 at 2:17 PM, Paulo Grahl <pgrahl at gmail.com> wrote:
> Thanks a lot!
> This saved me probably from a couple of afternoons digging into maxima's help!
> (as I am still learning it)
>
> one additional question:
> If I
> declare([ x,y] sym_matrix);
> How do I make maxima recognize that x+y , x.y, x^^2, etc.
> also are sym_matrix ?
> thanks.
> Paulo