b=vector([1,0,5,3]);
A=matrix(QQ,[[1,4],[3,4],[5,1]]);
or equivalently.
A=matrix(QQ,3,2,[1,4,3,4,5,1]);
The first optional parameter is the field. We will need the following ones:
ZZ - integers,
QQ - rationals,
RR - real numbers,
CC - complex numbers,
GF(k) - the unique finite field with k elements.
The optional second and third parameters are the numbers of rows and columns.
A.rref() or A.echelon_form()
transpose(A)
A*B
A^{-1} or A.inverse()
A.augment(B) - horizontally
A.stack(B, subdivide=True) - vertically
The optional parameter subdivide=True is for aestheticeal purposes only - it draws a line between the parts.