--- imach/src/imach.c	2002/11/19 14:08:13	1.61
+++ imach/src/imach.c	2002/11/20 07:56:21	1.62
@@ -1,4 +1,4 @@
-/* $Id: imach.c,v 1.61 2002/11/19 14:08:13 brouard Exp $
+/* $Id: imach.c,v 1.62 2002/11/20 07:56:21 brouard Exp $
    Interpolated Markov Chain
 
   Short summary of the programme:
@@ -969,7 +969,7 @@ double func( double *x)
 	s1=s[mw[mi][i]][i];
 	s2=s[mw[mi+1][i]][i];
 	bbh=(double)bh[mi][i]/(double)stepm;
-	lli= (savm[s1][s2]>(double)1.e-8 ?(1.+bbh)*log(out[s1][s2])- bbh*log(savm[s1][s2]):log((1.-bbh)*out[s1][s2]));
+	lli= (savm[s1][s2]>(double)1.e-8 ?(1.+bbh)*log(out[s1][s2])- bbh*log(savm[s1][s2]):log((1.+bbh)*out[s1][s2]));
 	/*lli= (savm[s1][s2]>1.e-8 ?(1.+bbh)*log(out[s1][s2])- bbh*log(savm[s1][s2]):log((1.-bbh)*out[s1][s2]));*/
 	/*lli=(1.+bbh)*log(out[s1][s2])- bbh*log(savm[s1][s2]);*/
 	/*if(lli ==000.0)*/