function [ Xm, sX, aX, kX ] = estatistica( X )
% [ X_medio, std_X, assimetria_X, kurtose_X ] = estatistica( X )

Xm = mean( X );
if nargout>=2 || nargout==0
    sX = std( X );
end
if nargout>=3 || nargout==0
    aX = skewness( X );
end
if nargout>=4 || nargout==0
    kX = kurtosis( X );
end
if nargout==0
%     disp( Xm )
%     disp( sX )
    disp( sprintf( 'Xm=%.3g, sX=%.3g, assimetria=%.1f, kurtose=%.1f', ...
        Xm, sX, aX, kX ) )
    clear Xm
end