function functn(x,i,a)
      implicit double precision (a-h,o-z)
      dimension x(50),a(10)
      functn=a(1)*dexp(a(2)*x(i))+a(3)*dexp(a(4)*x(i))
      return
      end