Function magnetic_field_energy_density 

pub fn magnetic_field_energy_density(magnetic_flux_density: f64) -> Option<f64>

Computes magnetic field energy density using u_B = B² / (2μ0).