////////////////////////////////////////////////////
function dens(X,c)
// demande le vecteur et la couleur du trace
n=length(X)
h=n^(-1/5);a=sqrt(5);b=3/(4*a*h*n);
xmin=min(X)-a*h;xmax=max(X)+a*h;
A=[xmin:0.1:xmax];
B=zeros(A);
for j=1:length(A)
  C=abs((A(j)-X)/h);
  B(j)=b*sum((1-(C.^2)/5).*(C<a));
end
plot2d(A,B,c);



////////////////////////////////////////////////////
function edexp(n)
xbasc();
dens(grand(1,n,"exp",2),5)
X=[0:.01:10];
plot2d(X,0.5*exp(-0.5.*X),3);
legends(['Densite','Estimation'],[3,5],1);
xtitle('Estimation de la densite exponentielle')



////////////////////////////////////////////////////
function edgau(n)
xbasc();
dens(grand(1,n,"nor",0,1),5)
X=[-4:.01:4];
plot2d(X,exp(-0.5.*X.^2)/(sqrt(2*%pi)),3);
legends(['Densite','Estimation'],[3,5],1);
xtitle('Estimation de la densite gaussienne')
legends(['Densite','Estimation'],[3,5],1);