  */
 __device__ void integrate (const double t_start, const double t_end, const double pr, double* y) {
 	double Hmin = 0;
-	double Hmax = fabs(t_end - t_start);
 	double Hold = 0;
 #ifdef Gustafsson
 	double Hacc = 0;

#ifndef NEWTON_UNROLL
            if (NewtonDone) break;
#else //only break if it's at the end of the unroll
            if (NewtonDone && (sNewtonIter + 1) % sNewtonIter == 0) break;
#endif
        }
        if (!NewtonDone) {

#ifndef NEWTON_UNROLL
            if (NewtonDone) break;
#else //only break if it's at the end of the unroll
            if(NewtonDone && (NewtonIter + 1) % NEWTON_UNROLL == 0) break;
#endif
            if (NewtonIter >= NewtonMaxit)
            {
