function [mu,sigma]=stat(x,n)
xx=x(2:end-1);
nn=n(2:end-1);
N=sum(nn);
mu=sum(xx.*nn)/N;
var=sum(((xx-mu).^2).*nn)/N;
sigma=sqrt(var);