format short e
%load ../datalc2/myout4; y=myout4;
y=avsread('../huber/lc/out18.fld');
x=avsread('../huber/phtmtbl.fld');
m=mean(mean(x(85:89,50:54)));
z=mean(mean(y(85:89,50:54,:))); z=z(:); 
disp([m-mean(z) std(z)]);

