No video
Root Finding in Python
Рет қаралды 12,174
Күн бұрынLearn how to numerically find roots of complex equations in python.
Script can be found here: www.hageslab.com/Resources.ht...
Here we are using "Spyder" IDE with the numpy, scipy, and matplotlib libraries
Script:
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import fsolve
Arbitrary Example
def fxn1(x):
return x**2+3*x - np.exp(x**(1/2))
xlist = np.linspace(0,1,num=1000)
plt.figure(num=1,dpi=120)
plt.plot(xlist,fxn1(xlist),label="Fxn1")
plt.plot(xlist,xlist*0,"--",label="y=0")
plt.legend()
root = fsolve(fxn1,0.5)
def fxn2(Ts, Tinf, Ti, A, Rtot, h, eps, sigma=5.67e-8): # SI units
return (Tinf - Ti)/(A*Rtot) - (h+eps*sigma*(Ts+Tinf)*(Ts**2+Tinf**2))*(Tinf-Ts)
Tlist = np.linspace(283.15,308.15,num=1000)
plt.figure(num=2,dpi=120)
plt.plot(Tlist,fxn2(Tlist,283.15,308.15,1.8,0.25,2,.95),label="Fxn2") # SI units
plt.plot(Tlist,Tlist*0,"--",label="y=0")
plt.legend()
Ts = fsolve(fxn2,290,args=(283.15,308.15,1.8,.25,2,0.95)) # [K]
Ts_C = Tsf-273.15 # [C]
@marocaineaucanada7538 Жыл бұрын
MERCI BEAUCOUP TRES INTERESSANT
@javierruiz3876 2 жыл бұрын
What if I'm getting the error "Result from function call is not a proper array of floats."
@charles3943 2 жыл бұрын
what if the equation doesnt equal to zero? where should we define this specific value? for the heat transfer example...
@mohamedyusufmohamud8193 2 жыл бұрын
should we include sigma variable in our function even though it is a constant but you didn't add it .
@bishalkarki8178 2 жыл бұрын
You should zoom the content,not clearly visible.
@asifahmadi116 15 күн бұрын
first thanks a lot for the excellent video and good description then I have a question how can you find the roots of this equation: R = ∑_(k=1)^(k=100)▒∑_(j=-100)^(j=100)▒(e^12 〖 (12)〗^k )/k! k/(2π(w^2+1/4)) √(1-〖(w/(LkR^2 ))〗^2 ) ∆w ∆k , ∆k=1,∆w=LkR^2/100