Я стремлюсь создать class NumericalAntiderivative
, который можно использовать так:
auto f = [](double x){
return std::sin(x); // in general this may be a complicated expression that is not analytically integrable
};
NumericalAntiderivative ad(f,minXofInterest,maxXofInterest);
// the line above may be costly, but the line below should be efficient
std::cout << "Integral of f from 0 to 5 is " << (ad(5.0)-ad(0.0)) << std::endl;
Я описал некоторые идеи в этом математическом вопросе , но этот вопрос может быть более подходящим для SO. Мне просто нужны некоторые идеи / указатели, которые могут привести к хорошему способу достижения этой цели.