NumCpp  2.12.1
A Templatized Header Only C++ Implementation of the Python NumPy Library
simpson.hpp
Go to the documentation of this file.
1
32#pragma once
33
34#include <functional>
35
36#include "NumCpp/Core/Types.hpp"
37
38namespace nc::integrate
39{
40 //============================================================================
41 // Method Description:
51 inline double
52 simpson(const double low, const double high, const uint32 n, const std::function<double(double)>& f) noexcept
53 {
54 const double width = (high - low) / static_cast<double>(n);
55
56 double simpson_integral = 0.;
57 for (uint32 step = 0; step < n; ++step)
58 {
59 const double x1 = low + static_cast<double>(step) * width;
60 const double x2 = low + static_cast<double>(step + 1) * width;
61
62 simpson_integral += (x2 - x1) / 6. * (f(x1) + 4. * f(0.5 * (x1 + x2)) + f(x2));
63 }
64
65 return simpson_integral;
66 }
67} // namespace nc::integrate
Definition: gauss_legendre.hpp:43
double simpson(const double low, const double high, const uint32 n, const std::function< double(double)> &f) noexcept
Definition: simpson.hpp:52
dtype f(GeneratorType &generator, dtype inDofN, dtype inDofD)
Definition: f.hpp:56
std::uint32_t uint32
Definition: Types.hpp:40