pro gfni,gg,hh,cmphi,smphi,Gfi,dGfi

    nmax=N_ELEMENTS(cmphi)-1
    Gfi = fltarr(nmax+1,nmax+1)
    dGfi = fltarr(nmax+1,nmax+1)
    for n=0, nmax do begin
        for m=0, n do begin
            Gfi(n,m) = gg(n,m)*cmphi(m)+hh(n,m)*smphi(m)
            dGfi(n,m) = m*(-1.*gg(n,m)*smphi(m)+hh(n,m)*cmphi(m))
        endfor
    endfor

end