FUNCTION pleg, x, n

; returns the l=0..(n-1) Pl(x) legendre polynomials (all evaluated at x).
; uses Bonnet's recursion equation for legendre polynomials

	P = dblarr(n_elements(x), n)

    P[*,0] = 1.0d0
    if n GT 1 then P[*,1] = x
    
    for i = 2, n-1 do p[*,i] = ((2.d*i-1.d)*x*p[*,i-1] - (i-1.d)*p[*,i-2])/(double(i))

    return, reform(p)

END