2020-12-23 09:52:43 +00:00
|
|
|
"""
|
|
|
|
Implementing Secant method in Python
|
|
|
|
Author: dimgrichr
|
|
|
|
"""
|
|
|
|
from math import exp
|
|
|
|
|
|
|
|
|
|
|
|
def f(x: float) -> float:
|
|
|
|
"""
|
|
|
|
>>> f(5)
|
|
|
|
39.98652410600183
|
|
|
|
"""
|
|
|
|
return 8 * x - 2 * exp(-x)
|
|
|
|
|
|
|
|
|
|
|
|
def secant_method(lower_bound: float, upper_bound: float, repeats: int) -> float:
|
|
|
|
"""
|
|
|
|
>>> secant_method(1, 3, 2)
|
|
|
|
0.2139409276214589
|
|
|
|
"""
|
|
|
|
x0 = lower_bound
|
|
|
|
x1 = upper_bound
|
2022-10-13 16:03:06 +00:00
|
|
|
for _ in range(0, repeats):
|
2020-12-23 09:52:43 +00:00
|
|
|
x0, x1 = x1, x1 - (f(x1) * (x1 - x0)) / (f(x1) - f(x0))
|
|
|
|
return x1
|
|
|
|
|
|
|
|
|
|
|
|
if __name__ == "__main__":
|
2021-07-21 06:01:55 +00:00
|
|
|
print(f"Example: {secant_method(1, 3, 2)}")
|