> v:=array(1..8):for i to 8 do v[i]:=i^2 od:print(v); 
A:=array(1..2,1..2):A[1,1]:=cos(x):A[1,2]:=sin(x):A[2,1]:=-A[1,2]:A[2,2]:=A[1,1]:print(A); 

> map(diff,A,x);

> with(linalg);
A:=matrix(3,3,[5,-1,0,0,1,9,4,6,3]); 
B:=matrix(3,3,[[9,3,-6],[3,2,0],[9,8,7]]); 

> evalm(A&*B);

> evalm(evalm(A&*B)-evalm(B&*A));

> evalm(A&*B+B&*A);
evalm(A^5+B^5);

> Id:=&*();

> restart;
> 
with(linalg):
> A:=matrix(3,4,[0,2,1,-2,3,5,-5,1,2,4,-2,2]);
b:=vector([-2,1,2]);
linsolve(A,b);

> eqns:={x+y+2*z=9,3*x-2*y=6,4*x+7*y=0};

> A:=genmatrix(eqns,[x,y,z],'flag'); 
A:=genmatrix(eqns,[x,y,z],'b');

> restart;
> p:=2*x*1-4*x*2+3*x*3-4*x*4-11*x*5=28;

> solve(p,x);

>