%                           rms.m
%  Scope:   This MATLAB macro computes root mean square (RMS) of a sample.
%  Usage:   y = rms(x)
%  Description of parameters:
%           x  -  input, sample (row or column) vector with  m  elements
%           y  -  output, RMS value (normalized by  m-1); if the sample 
%                 vector has only one element then the result is 0.
%  Last update:  09/27/00
%  Copyright (C) 1996-00 by LL Consulting. All Rights Reserved.

function  y = rms(x)

[m,n] = size(x);
if  n ~= 1
   x = x';
   m = n;
end
if  m == 1
   y = 0.;
else
   avg = sum(x)/m;
   y = - avg * avg;
   for  i = 1:m
      y = x(i) * x(i) + y;
   end
   y = sqrt(y/(m-1));
end