Function inverse_square_intensity 

pub fn inverse_square_intensity(
    reference_intensity: f64,
    reference_distance: f64,
    target_distance: f64,
) -> Option<f64>

Computes inverse-square intensity from a reference value.

Formula: I2 = I1 * (r1 / r2)^2.