pro gfnj,Nmax,cmphj,smphj,gg,hh,Gfj,dGfj

    Gfj = fltarr(Nmax+1,Nmax+1)
    dGfj = fltarr(Nmax+1,Nmax+1)
    for n=0, Nmax do begin
        for m=0, n do begin
        Gfj(n,m) = gg(n,m)*cmphj(m)+hh(n,m)*smphj(m)
        dGfj(n,m) = m*(-1.*gg(n,m)*smphj(m)+hh(n,m)*cmphj(m))
        endfor
    endfor

end