--- imach/src/imach.c	2014/12/16 11:20:36	1.165
+++ imach/src/imach.c	2014/12/22 11:40:47	1.166
@@ -1,6 +1,9 @@
-/* $Id: imach.c,v 1.165 2014/12/16 11:20:36 brouard Exp $
+/* $Id: imach.c,v 1.166 2014/12/22 11:40:47 brouard Exp $
   $State: Exp $
   $Log: imach.c,v $
+  Revision 1.166  2014/12/22 11:40:47  brouard
+  *** empty log message ***
+
   Revision 1.165  2014/12/16 11:20:36  brouard
   Summary: After compiling on Visual C
 
@@ -577,11 +580,11 @@ typedef struct {
 #define ODIRSEPARATOR '\\'
 #endif
 
-/* $Id: imach.c,v 1.165 2014/12/16 11:20:36 brouard Exp $ */
+/* $Id: imach.c,v 1.166 2014/12/22 11:40:47 brouard Exp $ */
 /* $State: Exp $ */
 
 char version[]="Imach version 0.99, September 2014,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121)";
-char fullversion[]="$Revision: 1.165 $ $Date: 2014/12/16 11:20:36 $"; 
+char fullversion[]="$Revision: 1.166 $ $Date: 2014/12/22 11:40:47 $"; 
 char strstart[80];
 char optionfilext[10], optionfilefiname[FILENAMELENGTH];
 int erreur=0, nberr=0, nbwarn=0; /* Error number, number of errors number of warnings  */
@@ -695,7 +698,12 @@ static double maxarg1,maxarg2;
   
 #define SIGN(a,b) ((b)>0.0 ? fabs(a) : -fabs(a))
 #define rint(a) floor(a+0.5)
-
+/* http://www.thphys.uni-heidelberg.de/~robbers/cmbeasy/doc/html/myutils_8h-source.html */
+/* #define mytinydouble 1.0e-16 */
+/* #define DEQUAL(a,b) (fabs((a)-(b))<mytinydouble) */
+/* http://www.thphys.uni-heidelberg.de/~robbers/cmbeasy/doc/html/mynrutils_8h-source.html */
+/* static double dsqrarg; */
+/* #define DSQR(a) (DEQUAL((dsqrarg=(a)),0.0) ? 0.0 : dsqrarg*dsqrarg) */
 static double sqrarg;
 #define SQR(a) ((sqrarg=(a)) == 0.0 ? 0.0 :sqrarg*sqrarg)
 #define SWAP(a,b) {temp=(a);(a)=(b);(b)=temp;} 
@@ -4091,8 +4099,8 @@ fprintf(fichtm," \n<ul><li><b>Graphs</b>
 <img src=\"%s%d_2.png\">",stepm,subdirf2(optionfilefiname,"pe"),jj1,subdirf2(optionfilefiname,"pe"),jj1,subdirf2(optionfilefiname,"pe"),jj1); 
        /* Period (stable) prevalence in each health state */
        for(cpt=1; cpt<=nlstate;cpt++){
-	 fprintf(fichtm,"<br>- Convergence from cross-sectional prevalence in each state (1 to %d) to period (stable) prevalence in specific state %d <a href=\"%s%d_%d.png\">%s%d_%d.png</a><br> \
-<img src=\"%s%d_%d.png\">",nlstate, cpt, subdirf2(optionfilefiname,"p"),cpt,jj1,subdirf2(optionfilefiname,"p"),cpt,jj1,subdirf2(optionfilefiname,"p"),cpt,jj1);
+	 fprintf(fichtm,"<br>- Convergence to period (stable) prevalence in state %d. Or probability to be in state %d being in state (1 to %d) at different ages. <a href=\"%s%d_%d.png\">%s%d_%d.png</a><br> \
+<img src=\"%s%d_%d.png\">", cpt, cpt, nlstate, subdirf2(optionfilefiname,"p"),cpt,jj1,subdirf2(optionfilefiname,"p"),cpt,jj1,subdirf2(optionfilefiname,"p"),cpt,jj1);
        }
      for(cpt=1; cpt<=nlstate;cpt++) {
         fprintf(fichtm,"\n<br>- Life expectancy by health state (%d) at initial age and its decomposition into health expectancies in each alive state (1 to %d) : <a href=\"%s%d%d.png\">%s%d%d.png</a> <br> \