x_QX:=PolynomialAlgebra(Q); x_f:=x_QX.1^2-5; x_O :=MaximalOrder(x_f); BasisMatrix(x_O);