Program Approx
      Implicit NONE
c        real PI
c        parameter (PI=3.141593)
        real V0Gap2m(3) / 0.1, 1., 10./
        real x_min /0.0/, x_max /10.0/
        real step /0.1/
        real x, T, y(3),ka,kapa,hyp
        Integer i,j,imax,jmax
c function:
c local:
        imax=(x_max-x_min)/step
        x = 0
        Do i=1,imax
           x = x + step
           if(x .lt.1)then
             Do j=1,3
               kapa=sqrt(V0Gap2m(j)*(1.-x))
               y(j)=16.*x*(1-x)*exp(-2.*kapa)
             Enddo
             write(*,' (4f10.4)') x, y(1), y(2), y(3)
           else
             Do j=1,3
               ka=sqrt(V0Gap2m(j)*(x-1.))
               y(j)=1./(1+sin(ka)*sin(ka)/(4.*x*(x-1.)))
             Enddo
             write(*,' (4f10.4)') x, y(1), y(2), y(3)
           Endif
        Enddo
       end