pro csmphj,Nmax,phrj,cmphj,smphj

    cmphj=fltarr(Nmax+1)
    smphj=cmphj
    for m=0,Nmax do begin
        mph=m*phrj
        cmphj(m)=COS(mph)
        smphj(m)=SIN(mph)
    endfor

end