def ellipse_conf(fmean,fcov):
# ------------------------------------------------------------------------
# calcul et representation graphique d'une ellipse de confiance a 90#
# pour un v.a. gaussien bi-dimensionnel
#    de moyenne m et de matrice de covariance Q
# parametres :
#    fmean (vecteur colonne) : moyenne m
#    fcov : matrice de covariance Q
# bibliotheques requises :
#    import numpy as np
#    import matplotlib.pyplot as plt # fonctions d'affichage MATLAB-like
#    import math
# ------------------------------------------------------------------------
    a = 4.6052
    npoints = 100
    t = np.linspace(0,2*math.pi,npoints)
    U = np.linalg.cholesky(fcov)
    round = np.zeros((2,npoints))
    round[0,:] = np.cos(t)
    round[1,:] = np.sin(t)
    ellipse = math.sqrt(a)*np.dot(U,round)
    ellipse = ellipse+np.tile(fmean,(1,npoints))
    plt.plot(ellipse[0,:],ellipse[1,:],'-k')