coeff_dx_m1 = x+RATIONAL(-1.0,2.0);
      coeff_dx_0 = RATIONAL(-2.0,1.0)*x;
      coeff_dx_p1 = RATIONAL(1.0,2.0)+x;