diff --git a/include/teqp/finite_derivs.hpp b/include/teqp/finite_derivs.hpp index 5f6e37074ed9df2ce3e4d1ab20c90bef7abf379f..f95818dbe92ba5c809309bf2b4c20dee445e2b00 100644 --- a/include/teqp/finite_derivs.hpp +++ b/include/teqp/finite_derivs.hpp @@ -32,7 +32,7 @@ auto centered_diff(const Function &f, const Scalar x, const Scalar h) { using r = Scalar; static std::map<std::tuple<int, int>, DiffCoeffs> CentralDiffCoeffs = { {{1, 2}, {{-1,1}, {-r(1)/r(2), r(1)/r(2)}} }, - {{1, 4}, {{-2,-1,1,2}, {r(1)/ r(1), -r(2)/r(3), r(2)/r(3), -r(1)/r(12)}} }, + {{1, 4}, {{-2,-1,1,2}, {r(1)/ r(12), -r(2)/r(3), r(2)/r(3), -r(1)/r(12)}} }, {{1, 6}, {{-3,-2,-1,1,2,3}, {-r(1)/r(60), r(3)/r(20), -r(3)/r(4), r(3)/r(4), -r(3)/r(20), r(1)/r(60)}} }, {{2, 2}, {{-1,0,1}, {1, -2, 1} }},