king/prof_prep_taopirata/rms.m

function [rms_value]=rms(Field_in)
%agel FORGET 22/01/03
a=size(Field_in);
ndim=length(a);
n=1;
Field_tmp=Field_in.*Field_in;
for i=1:ndim
        n=n*a(i);
        Field_tmp=sum(Field_tmp,ndim-i+1);
end
rms_value=sqrt(Field_tmp/n);

1	function [rms_value]=rms(Field_in)
2	%agel FORGET 22/01/03
3	a=size(Field_in);
4	ndim=length(a);
5	n=1;
6	Field_tmp=Field_in.*Field_in;
7	for i=1:ndim
8	n=n*a(i);
9	Field_tmp=sum(Field_tmp,ndim-i+1);
10	end
11	rms_value=sqrt(Field_tmp/n);
12