Function compute_gradient 

pub fn compute_gradient(
    p: &Matrix<f64>,
    dx: f64,
    dy: f64,
) -> (Matrix<f64>, Matrix<f64>)

Computes the gradient of a scalar field.

Returns (∂p/∂x, ∂p/∂y)