:


if(test -r session.dat) then 
  n1=`grep Session_number session.dat|cut -f2 -d" "`
else 
  n1=1
fi 

n2=`echo $n1 1 - $3 + p| dc`

file=pcm_tab_"$n1"_"$n2"
echo FILE= $file 

lastpcm=`echo $3 1 - $2  \* $1 + p| dc`


echo  > $file
echo  x-axis: \"Pcm\"  from $1 to $lastpcm N_points= $3 >> $file
echo   Diff. cross section [pb] >> $file
 
 
pcm=$1
nstep=$3

while(test 0 -ne $nstep)
do

  ../change_mom $pcm $pcm
  ../vegas_run  5 10000 

  n=`grep Session_number session.dat|cut -f2 -d" "`


  res=`grep "< >"  prt_$n`
  res=` echo $res |cut -f3 -d" "`

  echo Pcm=$pcm[GeV]  sigma=$res[pb] 
  echo $res >>$file

  pcm=`echo  $pcm $2 +  p|dc`
  nstep=`expr $nstep - 1`
done   

