34 int m = m_initial_order;
35 double value_order_m = 0.;
36 double value_order_m_1 = m_quadrature(function, min, max, m);
40 value_order_m = value_order_m_1;
41 value_order_m_1 = m_quadrature(function, min, max, value_order_m, m);
42 diff = value_order_m_1 - value_order_m;
43 }
while (std::abs(diff / value_order_m) > m_relative_precion);
45 return value_order_m_1;