Function electric_field_energy_density 

pub fn electric_field_energy_density(electric_field: f64) -> Option<f64>

Computes electric field energy density using u_E = 0.5 * ε0 * E².