Test of Cp
> Cp(5,2,vec([0,0,1,0],[1,1,0,0]));
> CpMtest:=array(1..dimH,1..dimH);
>
for i from 1 to dimH do
for j from 1 to dimH do
CpMtest[i,j]:=int_product(Cp(1,2,Basis[i]),Cp(1,2,Basis[j]));
od:
od:
> print(CpMtest);
>