for i in RE*
do

echo "$i "

cat $i | gawk '{printf("%40.40f\n",$1) }' | sort -n |
gawk '{All[++N]=$1}
END {
     for(I=50; I<=(N-50) ;I++)
			print All[I]
}'  > sorted.$i

cat sorted.$i |
gawk '
{Sum += $1
 SumSq += $1*$1
 N++
}
END {
	srand();
	for(I=1;I<=N;I++)  {
	    X= normal(Sum/N,sd(SumSq,Sum,N)) 
	    
		printf("%40.40f\n",X)
	}
}
function sd(sumSq,sumX,n) {
	return sqrt((sumSq-((sumX*sumX) / n)) / (n-1))
}
function normal(m,s) {
   return m+box_muller()*s;
 }
 function box_muller(m,s,    n,x1,x2,w) {
   w=1;
   while (w >= 1) {
     x1= 2.0 * rand() - 1;
     x2= 2.0 * rand() - 1;
     w = x1*x1 + x2*x2};
   w = sqrt((-2.0 * log(w))/w);
   return x1 * w;
 }


' | sort -n > normal.$i

gnuplot<<EOF
set key top left 
set title "$i"
set terminal png large 
#set logscale y
#set noxtics
set ylabel "log(MRE)"
set xlabel "all MREs, sorted"
#set yrange [0.001:]
set output "png_$i.png"
#set pointsize 3
plot "sorted.$i" title "actual" with points 2,"normal.$i" title "normal" with linesp  12
EOF
done