pro csmphi,Nmax,phri,cmphi,smphi

    cmphi=fltarr(Nmax+1)
    smphi=cmphi
    for m=0,Nmax do begin
        mph=m*phri
        cmphi(m)=COS(mph)
        smphi(m)=SIN(mph)
    endfor

end