c----------------------------------------------------------------------
      subroutine EIGENF(nm,n,rho,wfctns,ev,omega,phi,phi1)
c constructs the eigenfunctions for a basis-set expansion on the 
c gridpoints rho with coefficents from the corresponding eigenvectors 
c
c author: Heiko Appel, 11.23.99
c     in: nm   - size of Hamiltonian and overlap matrices
c         n    - number of grid points for the eigenfunctions
c         rho  - scaled grid
c         ev   - eigenvectors (and therefore expansion coefficients for the
c                radial functions)
c         phi  - auxilary array which is filled with the basis set functions
c         phi1 - auxilary array used for normalization of the radial 
c                functions
c    out: wfctns  - array which is filled by the subroutine with the radial 
c                   functions
c  calls: --
      implicit real*8(a-h,o-z)
      real*8 rho(0:n),phi(0:nm,0:n)
      real*8 wfctns(0:nm,0:n)
      real*8 ev(0:nm,0:nm),phi1(0:n)
      real*8 omega,pi,norm,u0,r0