735a736,744
>         # mu* = mu2 + SIGMA_21 * SIGMA_11^-1(S-mu1)
>         # mu* = mu2 + SIGMA_21 * SIGMA_11^-1(résidus)
>         # mu* = mu2 + SIGMA_21 * invQt
>         # theta = [theta_x , theta_y , theta_. ]
>         # Ktest => SIGMA_12
>         # = exp(theta_.) * [ exp( -1/2 * ( exp(theta_x)*(x1_i - x2_j)^2 + exp(theta_y)*(y1_i - y2_j)^2 ) ) ]
>         # = sigma_2 * [ kernel(x1_i, x2_j) ]
>         # exp(theta_.) = sigma_2 
>         # exp(theta_x) = 1/l_x^2 => l_x = exp(-theta_x/2)
753a763,773
>         print("HYPERPARAMETERS !!!!!")
>         print("mean_func", self.theta[:switch])
>         print("kernel   ", self.theta[switch:-1])
>         # exp(theta_.) = sigma_2 
>         # exp(theta_x) = 1/l_x^2 => l_x = exp(-theta_x/2)
>         lx = np.exp(-self.theta[switch]/2)
>         ly = np.exp(-self.theta[switch+1]/2)
>         sig2 = np.exp(self.theta[-2])
>         print("l_x=%f \nl_y=%f"%( lx*3200, ly*3200) )
>         print("l_moy=%f"%(1/2.*(lx+ly)*3200))
>         print("sigma_2=%f"%sig2)