program debug
 implicit none
 double precision :: ao_overlap
 integer :: i,j
 do i=1,ao_num
  do j=1,ao_num
   print *, i,j, ao_overlap_matrix(i,j)
  enddo
 enddo
end