function testva(fun,varargin)
%
% testva(variabile_aleatoria,parametri)
% 
n = 20; % numero di intervalli
[v,density] = fun(varargin{:});
N = length(v); % numero di lanci
a = min(v);
b = max(v);
h = (b-a)/n;
x = linspace(a,b,n+1);
for i = 1:n
  f(i) = length(find(v > x(i) & v <= x(i+1)));
end
f = f / (N - 1); % v=a non viene contato
hold on
bar((x(1:n)+x(2:n+1))/2,f/h)
x = linspace(a,b,100);
plot(x,density(x),'r')
