| 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;
|