anchor.diff:    !  recursive function involved in the bielectronic integral
     END_DOC
     implicit none
  +  include 'Utils/constants.include.F'
     integer, intent(in)            :: a,c,n_pt
  -  double precision, intent(in)   :: B_10(n_pt),B_01(n_pt),B_00(n_pt)
  -  double precision, intent(out)  :: res(n_pt)
  -  double precision               :: res2(n_pt)
  +  double precision, intent(in)   :: B_10(n_pt_max_integrals),B_01(n_pt_max_integrals),B_00(n_pt_max_integrals)
  +  double precision, intent(out)  :: res(n_pt_max_integrals)
  +  double precision               :: res2(n_pt_max_integrals)
     integer                        :: i
     
     if(c<0)then
  
positive.code:   BEGIN_DOC
    !  recursive function involved in the bielectronic integral
    END_DOC
    integer, intent(in)            :: c, n_pt
    double precision, intent(in)   :: B_10(n_pt),B_01(n_pt),B_00(n_pt)
    double precision, intent(out)  :: res(n_pt)
    integer                        :: i
    
    if(c==1)then
  
negative.code:   
    
    implicit none
    include 'Utils/constants.include.F'
    double precision               :: p,q
    integer                        :: a_x,b_x,c_x,d_x,a_y,b_y,c_y,d_y,a_z,b_z,c_z,d_z
    integer                        :: i, n_iter, n_pt, j
  