program main
	interger :: i,j,a,b,nn,ll,nst
    real(8), allocatable :: mat(:,:)

    mat(:,:)=0.d0
    ll=2
    nn=ll*ll
    nst=2^nn !!2的n次方个态
    do a=0,nst-1
    do i=0,nn-1
    

end