From e461df7573c2b7b7e26c965d8cf2d8e175d67378 Mon Sep 17 00:00:00 2001
From: Nguyễn Gia Phong <vn.mcsinyx@gmail.com>
Date: Mon, 16 Dec 2019 21:31:18 +0700
Subject: [usth/MATH2.2] Numerical Methods

The future starts now.
---
 usth/MATH2.2/labwork/1/ratpoison.py | 8 ++++++++
 1 file changed, 8 insertions(+)
 create mode 100644 usth/MATH2.2/labwork/1/ratpoison.py

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

diff --git a/usth/MATH2.2/labwork/1/ratpoison.py b/usth/MATH2.2/labwork/1/ratpoison.py
new file mode 100644
index 0000000..16aa65d
--- /dev/null
+++ b/usth/MATH2.2/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