> ?factor
> factor(6);

                                  6

> factor(x^2-y^2);

                           (x - y) (x + y)

> ifactor(6);

                               (2)  (3)

> ifactor(100!);

      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(linalg)
> ;
Warning, new definition for norm
Warning, new definition for trace

  [BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp,

        Wronskian, addcol, addrow, adj, adjoint, angle, augment,

        backsub, band, basis, bezout, blockmatrix, charmat, charpoly,

        cholesky, col, coldim, colspace, colspan, companion, concat,

        cond, copyinto, crossprod, curl, definite, delcols, delrows,

        det, diag, diverge, dotprod, eigenvals, eigenvalues,

        eigenvectors, eigenvects, entermatrix, equal, exponential,

        extend, ffgausselim, fibonacci, forwardsub, frobenius,

        gausselim, gaussjord, geneqns, genmatrix, grad, hadamard,

        hermite, hessian, hilbert, htranspose, ihermite, indexfunc,

        innerprod, intbasis, inverse, ismith, issimilar, iszero,

        jacobian, jordan, kernel, laplacian, leastsqrs, linsolve,

        matadd, matrix, minor, minpoly, mulcol, mulrow, multiply,

        norm, normalize, nullspace, orthog, permanent, pivot,

        potential, randmatrix, randvector, rank, ratform, row, rowdim,

        rowspace, rowspan, rref, scalarmul, singularvals, smith,

        stackmatrix, submatrix, subvector, sumbasis, swapcol, swaprow,

        sylvester, toeplitz, trace, transpose, vandermonde, vecpotent,

        vectdim, vector, wronskian]

> solve({2*x+y=3, x-y=0}, {x,y});

                            {y = 1, x = 1}

> solve({2*x+y^2=3, x-y=0}, {x,y});

                   {x = -3, y = -3}, {y = 1, x = 1}

> A := matrix(2,2,[1,2,3,4]);

                                 [1    2]
                            A := [      ]
                                 [3    4]

> b := vector([u,v]);

                             b := [u, v]

> ?linsolve
> linsolve(A,b);

                      [-2 u + v, 3/2 u - 1/2 v]

>