 
                 if ((!autonomous) )
                 {
-                    HG = direction*H*ros_Gamma[2];
+                    HG = direction*H*(0.7592633437920482E+00);
                     for (int i=0; i<NVAR; i++){
                         K(index,istage,i) += dFdT(index,i)*HG;
 		     }

         *
         */
        ros_Integrator_ros4(var, fix, Tstart, Tend, Texit,
                //  Rosenbrock method coefficients
                ros_M, ros_E, ros_A, ros_C, 
                ros_Alpha, ros_Gamma, ros_ELO, ros_NewF, 
                //  Integration parameters
                autonomous, vectorTol, Max_no_steps, 
                roundoff, Hmin, Hmax, Hstart, Hexit, 

		for (int i=0; i<NVAR; i++)		
			K(index,istage,i)  = Fcn0(index,i);

                if ((!autonomous) )
                {
                    HG = direction*H*(0.5728200000000000E+00);
                    for (int i=0; i<NVAR; i++){
                        K(index,istage,i) += dFdT(index,i)*HG;
		     }
