Function rbf_kernel 

pub fn rbf_kernel(
    x1: &[f64],
    x2: &[f64],
    amplitude: f64,
    length_scale: f64,
) -> f64

RBF (squared-exponential) kernel.

k(x, x’) = amplitude² · exp(-||x - x’||² / (2·length_scale²))