#!/bin/bash -x
make md1
time md1
lags="1 2 3 5 10 20"
prefix="test"
for i in $lags
do
  time ../../analyze2 -s 50 -m 500 -n 200000 -a 0 -l $i -c 5 -d 10 $prefix.out
  grep -v '^"t' $prefix.cond > $prefix.cond.lag_$i
done
xgraph $prefix.distr $prefix.condref
../../invert.awk 0.1 $prefix.cdistr > $prefix.inv
xgraph $prefix.inv $prefix.cond.lag_{1,2,3,5,10,20,40}
/bin/echo "set xrange [0:30]\n\
set xlabel \"x: delay (sec)\"\n\
set ylabel \"y: probability\"\n\
set term fig\n\
#set term post land\n\
set output \"md1_lag.fig\"\n\
plot \"$prefix.inv\" t \"unconditional\" w lines, \\" > foo.plot
for i in $lags
do
  if [ $i = 20 ]
  then
    /bin/echo "\"$prefix.cond.lag_$i\" t \"lag=$i\" w lines" >> foo.plot
  else
    /bin/echo "\"$prefix.cond.lag_$i\" t \"lag=$i\" w lines, \\" >> foo.plot
  fi
done
gnuplot foo.plot
gv md1_lag.eps