pub fn max_dt_diffusion(alpha: f64, dx: f64) -> f64
Computes the maximum stable time step for explicit diffusion (FTCS).