Function required_shield_thickness 

pub fn required_shield_thickness(
    linear_attenuation_coefficient: f64,
    initial_intensity: f64,
    target_intensity: f64,
) -> Option<f64>

Computes required shield thickness with x = ln(I0 / I) / mu.