Function coronal_loop_density 

pub fn coronal_loop_density(
    heating_rate: f64,
    loop_length: f64,
    temperature_k: f64,
) -> f64