> with(linalg);A:=array([[1,2],[1,3]]); b:=array([1,-2]); : > linsolve(A,b); > A:=matrix(2,2,[a11,a12,a21,a22]);