> 100!;

  9332621544394415268169923885626670049071596826438162146859296389\
        5217599993229915608941463976156518286253697920827223758251\
        185210916864000000000000000000000000

> ifactor(%);

      97     48     24     16      9      7      5      5      4
   (2)    (3)    (5)    (7)    (11)   (13)   (17)   (19)   (23)

             3      3      2      2      2      2
         (29)   (31)   (37)   (41)   (43)   (47)   (53)  (59)  (61)

         (67)  (71)  (73)  (79)  (83)  (89)  (97)

> with(LinearAlgebra);

  [Add, Adjoint, BackwardSubstitute, BandMatrix, Basis, BezoutMatrix,

        BidiagonalForm, BilinearForm, CharacteristicMatrix,

        CharacteristicPolynomial, Column, ColumnDimension,

        ColumnOperation, ColumnSpace, CompanionMatrix,

        ConditionNumber, ConstantMatrix, ConstantVector,

        CreatePermutation, CrossProduct, DeleteColumn, DeleteRow,

        Determinant, DiagonalMatrix, Dimension, Dimensions,

        DotProduct, Eigenvalues, Eigenvectors, Equal,

        ForwardSubstitute, FrobeniusForm, GaussianElimination,

        GenerateEquations, GenerateMatrix, GetResultDataType,

        GetResultShape, GivensRotationMatrix, GramSchmidt,

        HankelMatrix, HermiteForm, HermitianTranspose, HessenbergForm,

        HilbertMatrix, HouseholderMatrix, IdentityMatrix,

        IntersectionBasis, IsDefinite, IsOrthogonal, IsSimilar,

        IsUnitary, JordanBlockMatrix, JordanForm, LA_Main,

        LUDecomposition, LeastSquares, LinearSolve, Map, Map2,

        MatrixAdd, MatrixInverse, MatrixMatrixMultiply, MatrixNorm,

        MatrixScalarMultiply, MatrixVectorMultiply, MinimalPolynomial,

        Minor, Multiply, NoUserValue, Norm, Normalize, NullSpace,

        OuterProductMatrix, Permanent, Pivot, QRDecomposition,

        RandomMatrix, RandomVector, Rank, ReducedRowEchelonForm, Row,

        RowDimension, RowOperation, RowSpace, ScalarMatrix,

        ScalarMultiply, ScalarVector, SchurForm, SingularValues,

        SmithForm, SubMatrix, SubVector, SumBasis, SylvesterMatrix,

        ToeplitzMatrix, Trace, Transpose, TridiagonalForm, UnitVector,

        VandermondeMatrix, VectorAdd, VectorAngle,

        VectorMatrixMultiply, VectorNorm, VectorScalarMultiply,

        ZeroMatrix, ZeroVector, Zip]

> ?LinearSolve
> A := Matrix([[2,3]]);

                            A := [2    3]

> b := Vector([5]);

                               b := [5]

> LinearSolve(A,b);

                          [5/2 - 3/2 _t0[2]]
                          [                ]
                          [     _t0[2]     ]

> A2 := Matrix([[2]]);

                              A2 := [2]

> b2 := Vector([[2]]);

                              b2 := [2]

> LinearSolve(A2,b2);

                                 [1]

>