summaryrefslogtreecommitdiff
path: root/math/simpson.cpp
diff options
context:
space:
mode:
authormzuenni <michi.zuendorf@gmail.com>2022-06-27 17:19:28 +0200
committermzuenni <michi.zuendorf@gmail.com>2022-06-27 17:19:28 +0200
commit5ab8a5088b729a9953b8dff1b2a985dc8fb2098b (patch)
treeed40d6936c0e9eee40ba62751cbf99ecddbaddc2 /math/simpson.cpp
parentadabbad9c51cf7cd3874bfde8eac1fbcf84fec10 (diff)
updated tcr
Diffstat (limited to 'math/simpson.cpp')
-rw-r--r--math/simpson.cpp12
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);
}