************************************************************************
*
*     ExternalEvolutionMatrixPh2Ph.f:
*
*     This function returns the "external" evolution matrix on the
*     interpolation grid. This is possible only if the internal grids are 
*     used or if one single external grid is provided.
*
************************************************************************
      function ExternalEvolutionMatrixPh2Ph(i,j,alpha,beta)
*
      implicit none
*
      include "../commons/grid.h"
      include "../commons/EvolutionOperator.h"
**
*     Input Variables
*
      integer i,j
      integer alpha,beta
**
*     Output Variables
*
      double precision ExternalEvolutionMatrixPh2Ph
*
      ExternalEvolutionMatrixPh2Ph = Ph2PhQCD(0,i,j,alpha,beta)
*
      return
      end