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	king	1.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