Function solid 

pub(crate) fn solid<const SECOND: bool>(
    l: i32,
    m: i32,
    position: [f64; 3],
) -> Solid

r^l P_l^m(z/r) exp(i m phi), zero for |m| > l. The second derivative is compiled out unless requested.