From 82e6cf7d1046d6cee16f7e8b044ec33e7ec6c4b7 Mon Sep 17 00:00:00 2001
From: Nguyễn Gia Phong <vn.mcsinyx@gmail.com>
Date: Sun, 16 Feb 2020 14:26:55 +0700
Subject: [usth] Numerical Method is MATH2.4

---
 usth/MATH2.4/labwork/1/ratpoison.py | 8 ++++++++
 1 file changed, 8 insertions(+)
 create mode 100644 usth/MATH2.4/labwork/1/ratpoison.py

(limited to 'usth/MATH2.4/labwork/1/ratpoison.py')

diff --git a/usth/MATH2.4/labwork/1/ratpoison.py b/usth/MATH2.4/labwork/1/ratpoison.py
new file mode 100644
index 0000000..16aa65d
--- /dev/null
+++ b/usth/MATH2.4/labwork/1/ratpoison.py
@@ -0,0 +1,8 @@
+def ratpoison(f, df, x, es=10**-8, imax=20):
+    ea, i = 1, 0
+    while ea > es and i < imax:
+        i += 1
+        xold, x = x, x - f(x)/df(x)
+        if f(x) == 0: return x, 0, 0, i
+        if x: ea = abs((x - xold) / x)
+    return x, f(x), ea, i
-- 
cgit 1.4.1