Calculate all normalized associated Legendre polynomials and their
first and second derivatives for 0 <= l <= lmax and 0 <= m <= l for |x| <= 1. The
norm parameter specifies which normalization is used. The array index
of P_l^m is given by l(l+1)/2+m. To include or exclude the
Condon-Shortley phase factor of (-1)^m, set the fourth parameter to
either -1 or 1, respectively.
