pub fn diffusion_2d_explicit( u: &[Vec<f64>], dx: f64, dy: f64, dt: f64, alpha: f64, ) -> Vec<Vec<f64>>