share/linearalgebra/linearalgebra.dem

   1 load("linearalgebra");
   2
   3 m : matrix([1,2],[1,2]);
   4
   5 nullspace(m);
   6
   7 columnspace(m);
   8
   9 ptriangularize(m-z*ident(2),z);
  10
  11 m : matrix([1,2,3],[4,5,6],[7,8,9]) - z * ident(3);
  12
  13 mm : ptriangularize(m,z);
  14
  15 algebraic : true;
  16
  17 tellrat(mm[3,3]);
  18
  19 mm : ratsimp(mm);
  20
  21 nullspace(mm);
  22
  23 m : matrix([1,2,3,4],[5,6,7,8],[9,10,11,12],[13,14,15,16]);
  24
  25 column_space(m);
  26
  27 apply('orthogonal_complement, args(nullspace(transpose(m))));
  28
  29 untellrat(z);