program debug implicit none double precision :: eplf_integral, ao_overlap PROVIDE ao_prim_num_max integer :: i,j read(*,*) i,j print *, eplf_integral(i,j,eplf_gamma,point) print *, ao_overlap(i,j) end