Scribd
Upload
6 views3 pages

Bisection

Uploaded by

Michael John P. Denque
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
6 views3 pages

Bisection

Uploaded by

Michael John P. Denque
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 3

import math

""" MICHAEL JOHN DENQUE'S WORKS,

CRDT: ONLINE TOTURIAL & SOURCE

"""

def bisection_method(func, lower, upper, max_iter):

#TABLE

print(f"{'Iter':<5}{'xl':<20}{'xu':<20}{'xr':<20}{'ea%':<25}")

print('-'*70)

for i in range(max_iter):

midpoint = (lower + upper) / 2 # Calculate midpoint/xr

f_mid = func(midpoint) #f(xr)

# Calculate the error/ea

error = abs((midpoint - lower) / midpoint) * 100

# Print the current iteration's results

print(f"{i + 1:<5}{lower:<20.5f}{upper:<20.5f}{midpoint:<20.5f}{error:<15.5f}")

if error <= 0.05: #KANI RAY GI DUNGAG

break

if f_mid == 0:

return midpoint
# Narrow down the interval

if func(lower) * f_mid < 0:

upper = midpoint # Root is in the lower half

else:

lower = midpoint # Root is in the upper half

return (lower + upper) / 2 # Return the best estimate after max iterations

# Example usage

if __name__ == "__main__":

# Input function as a string

while True:

func_input = input("Enter the function of x or 'q' to quit: ")

if func_input.lower() == 'q':

print("\n OKAY BYEE!!!")

break

else:

# Define the function using eval

def func(x):

return eval(func_input)

# Input lower and upper limits

lower = float(input("Enter the lower limit: "))

upper = float(input("Enter the upper limit: "))


# Input maximum iterations

max_iterations = int(input("Enter the maximum number of iterations: "))

try:

root = bisection_method(func, lower, upper, max_iterations)

print(f"\nThe root is approximately: {root:.5f}")

except ValueError as e:

print(e)

You might also like