def make_sq_exp_kernel(L0):

    def sq_exp_kernel(x1, x2):
        x1 = np.expand_dims(x1, 2)
        x2 = x2.T
        return np.exp(((- np.sum(((x1 - x2) ** 2), axis=1)) / (2 * (L0 ** 2))))
    return sq_exp_kernel
def make_sq_exp_kernel(L0):

    def sq_exp_kernel(x1, x2):
        x1 = np.expand_dims(x1, 2)
        x2 = x2.T
        return np.exp(((- np.sum(((x1 - x2) ** 2), axis=1)) / (2 * (L0 ** 2))))
    return sq_exp_kernel