Function kernel_matrix 

pub fn kernel_matrix(
    x: &[Vec<f64>],
    amplitude: f64,
    length_scale: f64,
) -> Matrix

Compute the kernel matrix K[i, j] = k(X[i], X[j]).