function ellisse1(xbar,c2,S)
%function ellisse1(xbar,c2,S)
% Grafico dell'ellisse 
%
%      (t - xbar)^T inv(S) (t-xbar) = c2   centro: xbar=[x1;x2] 
%
%                                     con c2 = \chi_2^2 = 1.39
%                                         (per test normalita')
%   x1   centro prima variabile (scalare)
%   x2   centro seconda variabile (scalare)
%   S    matrice di distanza statistica 
%        (per test normalita', e' la m. covarianza)

x1=xbar(1); x2=xbar(2);
%c2=1.39;
[v,e]=eig(S); 
theta = 0:.01:2*pi;
ntheta=length(theta);
T=[x1;x2]*ones(1,ntheta)+v*sqrt(c2*e)*[cos(theta);sin(theta)];

plot(T(1,1),T(2,1),'.');hold on
for k=2:ntheta,
plot(T(1,k),T(2,k),'.')
end
hold off

