c      ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c      cc                                                          ccc
c      cc      Fitting Formulae for Photoneutrino Process          ccc
c      cc                                  (C) Naoki Itoh          ccc
c      cc                                                          ccc
c      ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
        implicit real*8 ( a-h, o-z )
        real*8 lamb
        parameter ( nn = 10000 )
        dimension temp( nn ), rhom( nn ), lamb( nn )
     &           ,qp  ( nn ), fp  ( nn ), xi  ( nn )
     &           ,qpho( nn )
c
       call init (       temp, rhom                    )
       call clamb( lamb, temp                          )
       call cxi  ( xi  ,       rhom, lamb              ) 
       call cqp  ( qp  ,       rhom, lamb              )  
       call cfp  ( fp  , temp,       lamb, xi          )  
       call cqpho( qpho,       rhom, lamb, qp   , fp   )  
       call outpt( qpho, temp, rhom )
c
        stop
        end