Test of Cp

> Cp(5,2,vec([0,0,1,0],[1,1,0,0]));

0

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

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

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

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

>