pub fn heat_conduction_1d(
    slice: &[Float],
    delta_time: Float,
    depth_steps: &[Float],
    diffusivities: &[Float]
) -> Vec<Float>