pro gfni06,nmax,cmphi,smphi,gg,hh,Gfi,dGfi

    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