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