Function lap1_backward 

fn lap1_backward(
    f: &[f64],
    mask: &[u8],
    l: usize,
    a: usize,
    n_axis: usize,
    coord: usize,
    _n_total: usize,
) -> f64

Backward one-sided second derivative (matching MATLAB’s bd/bf functions)