function v = var(varargin)
%VAR  Variance
%
%   CALL: v = var(X,dim);
%
%     v  = Sample variance (second central moment)
%     X  = data vector or matrix
%   dim = dimension to sum across. (default 1'st non-singleton dimension of X)
% 
% Example:
%    R = wgumbrnd(2,2,[],100,2);
%   var(R)
%
% See also: wskewness, wkurtosis, mean

v = std(varargin{:}).^2;
return