You may try
template <typename T>
double RombergIntegral(const double& a, const double& b, static T (*f)(T), const int& order)
{
}
class MyModel
{
public:
static double function(double t)
{
}
double integral_function(double t)
{
return RombergIntegral < double > (0.0, t, &MyModel::function, 5);
}
};