function gauss5(x)
       common/pawpar/p(8)

       x1=x-p(2)
       gauss5 = 
     *  5.0e-3*p(1)*exp(-0.5*(x1)**2/p(3)**2)/sqrt(2*3.1415926)/p(3)
     * + p(4) + p(5)*x1 + p(6)*x1**2 + p(7)*x1**3 + p(8)*x1**4
       return
       end