function rmsout=rms(v)
%rms Root Mean Square
%rms(v) returns the root mean square of the elements of
%a column vector v. If v is a matrix, then rms(v) returns
%a row vector such that each element is the root mean square
%of the elements in the corresponding column of v.
vs=v.^2;
s=size(v);
rmsout=sqrt(sum(vs,1)/s(1));