% modifs
% GENERATION DES CLASSES
% volume minimum
%%%%%%%Vmin=1/200;
% definition du nombre de classes
%%%%%%%nc=50;
% raison geometrique pour la descretisation des classes
%%%%%%%s=1.2;

%agrégation pure
%FONC= FONC + TERME1-VrupY;
% noyau
%freq=1;

% -----------------------------------------------

% Distribution initiale des classes
di=exp(-v);
N=zeros(nc,1);
for i=1:nc-1
N(i)=eval(int(di,v,X(i),X(i+1)));
end
N(nc)=eval(int(di,v,X(nc),s*X(nc)));
semilogx(X,N,'-')
hold on

% solution analytique t=1

% solution numerique t=1