Function second_derivative 

pub fn second_derivative(u: &[f64], dx: f64) -> Vec<f64>

Computes the second derivative of u with the standard 3-point Laplacian stencil.

Interior: (u[i-1] - 2 u[i] + u[i+1]) / dx². Boundaries are set to zero.