From 5c28e0abb061a32f61ad0626cbba64a1fbcd4576 Mon Sep 17 00:00:00 2001
From: hbui <hong-phuc.bui@htwsaar.de>
Date: Fri, 19 Jul 2024 15:23:57 +0200
Subject: [PATCH] fix unittest
---
num-int/src/numint/RiemannSum.py | 19 +++++--------------
1 files changed, 5 insertions(+), 14 deletions(-)
diff --git a/num-int/src/numint/RiemannSum.py b/num-int/src/numint/RiemannSum.py
index 12febc7..826ac3f 100644
--- a/num-int/src/numint/RiemannSum.py
+++ b/num-int/src/numint/RiemannSum.py
@@ -16,26 +16,22 @@
def numint(f, a, b, epsilon=1e-3):
dx = b - a
+ x = [a, b] # debug only
y = [f(a), f(b)]
s_left = dx * y[0]
s_right = dx * y[1]
- (n, parts) = (1, 1)
+ (n, parts) = (0, 1)
while abs(s_right - s_left) > epsilon:
- x = [a, b] # debug only
parts = 2 * parts
dx = dx / 2
n += 1
- print(n, " → ", parts)
+ # print(n, " → ", parts)
for i in range(1, parts, 2):
- x.insert(i, a + i * dx) # debug only
+ x.insert(i, a + i * dx)
y.insert(i, f(a + i * dx))
- # print(i, end=" ")
- # print()
- # print(" ", x)
- # print(" ", y)
s_left = sum(y[0:-1]) * dx
s_right = sum(y[1:]) * dx
- return s_left, s_right
+ return x, y, s_left, s_right, n
def numint_compact(f, a, b, epsilon=1e-3):
@@ -49,15 +45,10 @@
parts = 2 * parts
dx = dx / 2
n += 1
- print(n, " → ", parts)
for i in range(1, parts, 2):
- # print(i, end=" ")
y = f(a + i*dx)
y_l += y
y_r += y
- print()
- # print(" ", x)
- # print(" ", y)
s_left = y_l * dx
s_right = y_r * dx
return s_left, s_right
--
Gitblit v1.10.0-SNAPSHOT