diff options
| author | mzuenni <michi.zuendorf@gmail.com> | 2022-06-27 17:19:28 +0200 |
|---|---|---|
| committer | mzuenni <michi.zuendorf@gmail.com> | 2022-06-27 17:19:28 +0200 |
| commit | 5ab8a5088b729a9953b8dff1b2a985dc8fb2098b (patch) | |
| tree | ed40d6936c0e9eee40ba62751cbf99ecddbaddc2 /math/simpson.cpp | |
| parent | adabbad9c51cf7cd3874bfde8eac1fbcf84fec10 (diff) | |
updated tcr
Diffstat (limited to 'math/simpson.cpp')
| -rw-r--r-- | math/simpson.cpp | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/math/simpson.cpp b/math/simpson.cpp index dd887e2..a99b911 100644 --- a/math/simpson.cpp +++ b/math/simpson.cpp @@ -1,12 +1,12 @@ -double f(double x) { return x; } +double f(double x) {return x;} double simps(double a, double b) { - return (f(a) + 4.0 * f((a + b) / 2.0) + f(b)) * (b - a) / 6.0; + return (f(a) + 4.0 * f((a + b) / 2.0) + f(b)) * (b - a) / 6.0; } double integrate(double a, double b) { - double m = (a + b) / 2.0; - double l = simps(a, m), r = simps(m, b), tot = simps(a, b); - if (abs(l + r - tot) < EPSILON) return tot; - return integrate(a, m) + integrate(m, b); + double m = (a + b) / 2.0; + double l = simps(a, m), r = simps(m, b), tot = simps(a, b); + if (abs(l + r - tot) < EPS) return tot; + return integrate(a, m) + integrate(m, b); } |
