LHAPDF/lhapdf5.3.1/eksarp.f

   1       function eksarp(x,v2,x1,x116,ptm,indpi)
   2       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
   3       dimension ptm(10)
   4
   5       z=v2
   6       xx=x-x1
   7       if (x.le.x116) then
   8          eksarp=ptm(1)+ptm(2)*z+ptm(3)*z**2+ptm(4)*xx+ptm(5)*xx**2
   9          else
  10             z1= log(x116/x1)
  11             xx1=x116-x1
  12             ff0=ptm(1)+ptm(2)*z1+ptm(3)*z1**2+ptm(4)*xx1+ptm(5)*xx1**2
  13             xx16=x-x116
  14             z16= log(x/x116)
  15             qexp=19.d0+(indpi-1)*(indpi-2)*8.d0/2.d0
  16             eksarp=ff0+ptm(6)*xx16+ptm(7)*xx16**2+ptm(8)*xx16**3
  17      +             +ptm(9)*xx16**qexp+ptm(10)*z16
  18       endif
  19 11    return
  20       end
  21
  22