pub fn directivity(pattern: &[(f32, f32, f32)]) -> f32
Compute directivity from a radiation pattern. D = 4*pi * max(P) / integral(P * sin(theta) dtheta dphi)