x=avsread('phtmebl.fld');
y=avsread('fbplchuber.fld');

r1=mean(mean(y(69:71,45:47,:),1),2); r1=r1(:);
m1=mean(r1); s1=sqrt(var(r1)); b1=m1-mean(mean(x(69:71,45:47)));
vv1=sqrt( mean((r1-m1).^4)/50 );

r2=mean(mean(y(37:39,62:64,:),1),2); r2=r2(:);
m2=mean(r2); s2=sqrt(var(r2)); b2=m2-mean(mean(x(37:39,62:64)));
vv2=sqrt( mean((r2-m2).^4)/50 );

r3=mean(mean(y(89:91,66:68,:),1),2); r3=r3(:);
m3=mean(r3); s3=sqrt(var(r3)); b3=m3-mean(mean(x(89:92,66:69)));
vv3=sqrt( mean((r3-m3).^4)/50 );

r4=mean(mean(y(29:31,64:66,:),1),2); r4=r4(:);
m4=mean(r4); s4=sqrt(var(r4)); b4=m4-mean(mean(x(29:31,64:66)));
vv4=sqrt( mean((r4-m4).^4)/50 );

r5=mean(mean(y(105:107,60:62,:),1),2); r5=r5(:);
m5=mean(r5); s5=sqrt(var(r5)); b5=m5-mean(mean(x(105:107,60:62)));
vv5=sqrt( mean((r5-m5).^4)/50 );

r6=mean(mean(y(66:68,70:72,:),1),2); r6=r6(:);
m6=mean(r6); s6=sqrt(var(r6)); b6=m6-mean(mean(x(66:68,70:72)));
vv6=sqrt( mean((r6-m6).^4)/50 );

r7=mean(mean(y(66:68,88:90,:),1),2); r7=r7(:);
m7=mean(r7); s7=sqrt(var(r7)); b7=m7-mean(mean(x(66:68,88:90)));
vv7=sqrt( mean((r7-m7).^4)/50 );

%[b1, s1, b2, s2, b3, s3, b4, s4, b5, s5, b6, s6, b7, s7]
[s1,s2,s6,s7,s5,s3,s4].^2/1e6
n1=mean(mean(x(69:71,45:47))); n2=mean(mean(x(37:39,62:64)));
n3=mean(mean(x(89:92,66:69))); n4=mean(mean(x(29:31,64:66)));
n5=mean(mean(x(105:107,60:62))); n6=mean(mean(x(66:68,70:72)));
n7=mean(mean(x(66:68,88:90)));
[s1/n1,s2/n2,s6/n6,s7/n7,s5/n5,s3/n3,s4/n4]
