diff options
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); } |
