;
; PURPOSE: Calculated Associated Lagendre Polynomals 
; WRITTEN: Jan 1992	Xuepu Zhao
; 
pro pdp06,cth,Nmax,P,dP

    if N_params() LT 1 then begin
       print,'pro pdp06,cth,Nmax,P,dP'
       return
    endif
    jth=n_elements(cth)
    P = fltarr(Nmax+1,Nmax+1,jth)
    dP = fltarr(Nmax+1,Nmax+1,jth)
    for j=0,jth-1 do begin
        x = cth(j)
        pdpj06,x,Nmax,Pj,dPj
        P(*,*,j) = Pj(*,*)
        dP(*,*,j) = dPj(*,*)
    endfor

end