Test of Inner Product
> Mtest:=array(1..dimH,1..dimH);
> 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);
Test passed!!!!!