Test of Inner Product

> Mtest:=array(1..dimH,1..dimH);

Mtest := array(1 .. 36,1 .. 36,[])

> for i from 1 to dimH do
for j from 1 to dimH do
Mtest[i,j]:=int_product(Basis[i]-Basis[1],Basis[j]);
od:
od:

> print(Mtest);

matrix([[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...

Test passed!!!!!