Python Root finding code

import mathdef derivative(f, x): h=1e-8 return (f(x+h)-f(x))/hdef solver(f, x0, epsilon, max_iter): xn=x0 for n in range(0,max_iter): y=f(xn) if abs(y)<epsilon: return xn slope=derivative(f,xn) if(slope==0): return None xn=xn-y/slope return Nonedef loop(f, L_bound, R_bound, increment): solutions=[] while L_bound<=R_bound: solution=solver(f, L_bound, 1e-10, 1000) if solution is not None: solution=round(solution,4) if solution not in solutions: solutions.append(solution) L_bound+=increment print(sorted(solutions)) print(“we found “+str(len(solutions))+” solutions!”)equation=””def f(x): try: y=eval(equation) except ZeroDivisionError: y= 1e-10 return yequation=”x**2–9"loop(f,-100,100,0.5)

Source: medium.com

Similar codes

147 votes188 views6 answers
Python if file exists
Python
1 vote215 views1 answer
Django print settings
Python Django
0 vote182 views1 answer
Modeltranslation
Python
24 votes412 views6 answers
Django update model
Python Django
-1 vote186 views1 answer
Example exponential distribution python
Python
0 vote186 views1 answer
How to display text on boxplot in python
Python
1 vote240 views1 answer
Model.predict([x_test]) error
Python
0 vote198 views1 answer
How to run django server outside world
Python
0 vote185 views1 answer
Check if a numpy array contains only 1's python
Python
5 votes204 views1 answer
Pandas version from python script
Python