编辑代码

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次方个态

end