c     given function  h small we compute the m-th order cummulant
c     for a function f, Kf_m  at h. 
c     h is given in X=(x^2-c)/r coordinates

      subroutine kum( n, m, h, c, r,alpha, beta, gamma, 
     *  dffalpm, workf, kummfh)
      integer n, m
      double precision f(n), h(n),kummfh(n)
      double precision c,r
      double precision workf(n), newf(n), newftemp(n)
      double precision aux
      double precision alpha, beta, gamma, dffalpm(n)
c     computes h(gamma*X + beta)* dffalpm(X)
      call mataffine(n,h,newftemp,beta,gamma)
      call prod(n,newftemp,dffalpm, newf)
c-----------------------------------------------
c     computes h(Gammaf(X))
      call compos(n,workf,h,newftemp)

c-------------------------------------
c     we compute K_m h (X)
      call matadd(n,newftemp,newf,kummfh)
      aux=alpha**m
      call matsca(n,kummfh,1.d0/aux)
      return
      end