Fix lyapunov

src/lyapunov.c

@@ -88,7 +88,7 @@ int lyapunov(
   int i;
   // init average
   for (i=0; i<MATSIZE; i++){
-    lyapunov_avg=0;
+    lyapunov_avg[i]=0;
   }
 
   // save times at which Lyapunov exponents are computed
@@ -133,6 +133,7 @@ int lyapunov(
       }
 
       // print
+      fprintf(stderr,"% .15e\n",time);
       for(i=0; i<MATSIZE; i++){
 	printf("% .15e % .15e % .15e\n",time, lyapunov[i], lyapunov_avg[i]);
       }