|  | LIB "multigrading.lib";
intmat P[2][4] =
1,3,4,6,
2,4,5,7;
// should be a matrix with two columns
// for example: [-2, 1, 0, 0], [3, -3, 0, 1]
intmat U = integralSection(P);
print(U);
==>     -1     1
==>     -2     1
==>      2    -1
==>      0     0
print(P * U);
==>      1     0
==>      0     1
kill U;
 |