rm x_*.dat

for i in  10.11 1.04 13.29 15.67  249.0 9.0 7.33
do
echo $i | gawk '
BEGIN{ i=0.01
	   minpf=0.01
	   minpd=0.5
}
function prec(np,pd,pf) { return 1/(1 + np * pf/pd) }
{   np=$1
	for(pf=minpf; pf<=1; pf += i) 
		for(pd=pf; pd<=1; pd += i)  
				printf("%9.5f\t%9.5f\t%9.5f\t%9.5f\n", 
						prec(np,pd,pf), np,pd,pf) 
}' | sort  -n +0 | cat -n  > x_$i.dat
done
#gnuplot<<EOF
#	set logscale x
#set xrange [0.05:20]
#set xtics (0.01,0.1,1,2,3,4,5,6,7,8,10,20)
#plot 
#"x_1.dat" using 2:1 title "np", "x_1.dat" using 3:1 title "pd", "x_1.dat" using 4:1 title "pf"
#"x_5.dat" using 2:1 title "np", "x_5.dat" using 3:1 title "pd", "x_5.dat" using 4:1 title "pf"
#"x_10.dat" using 2:1 title "np", "x_10.dat" using 3:1 title "pd", "x_10.dat" using 4:1 title "pf"
#"x_20.dat" using 2:1 title "np", "x_20.dat" using 3:1 title "pd", "x_20.dat" using 4:1 title "pf"
#"x_100.dat" using 2:1 title "np", "x_100.dat" using 3:1 title "pd", "x_100.dat" using 4:1 title "pf"
#"x_200.dat" using 2:1 title "np", "x_200.dat" using 3:1 title "pd", "x_200.dat" using 4:1 title "pf"
#pause 100
#EOF
#
cat  <<EOF > dots.dat
2832 0.2025 PC3
2992  0.1960 KC3
2992 0.2064 CM1
2992 0.703  KC4 yes
3248 0.1527 PC1
3534 0.3155 PC4 yes
3988 0.269  MW1
4036 0.0202 PC2 yes
EOF
gnuplot<<EOF
set terminal postscript eps 20
set output "x.eps"
set xlabel ""
set size 1,1
#set yrange [0.005:1]
#set logscale y
set label "neg/pos= 1" at 200,0.6
set label "neg/pos= 7.33" at 200,0.17
set label "neg/pos= 15.67" at 200,0.04
set label "neg/pos= 249" at 2000,0.04
set xtics ("56" 2832,  "60" 2992,  "64" 3248,  "70" 3534, "80" 4000)
set ytics (0.02,0.15,0.20,0.26,0.32,0.7,1)
#set xrange [1:6500]
	set yrange [-0.01:1]
set ylabel "precsion"
set xlabel "calculations (sorted by precision)"
set key bottom right
#set arrow 4 from 0,0.32 to 4000,0.32 nohead linetype 0
#set arrow 5 from 4000,0.32 to 4000,0.004 nohead linetype 0
#set arrow 6 from 3700,0.2 to 3700,0.004 nohead linetype 0
#set arrow 7 from 3200,0.15 to 3200,0.004 nohead linetype 0
#set ytics (0.01,0.05,0.2,0.4,0.6,0.8)
#set xrange [0:25000]
set nokey
plot \
"x_10.11.dat"  with dots ,\
"x_1.04.dat" with dots  ,\
"x_13.29.dat" with dots ,\
"x_15.67.dat"  with dots,\
"x_249.0.dat"  with dots  ,\
"x_9.0.dat"  with dots  ,\
"x_7.33.dat"  with dots,\
"dots.dat" with points  pointtype 5 pointsize 2
EOF
epstopdf x.eps