Regula是一个在设定的时间间隔内找到最小值的函数(查找函数的根)。
def Regula(f,Tol,Nmax,a,b):
N = 1
while N < Nmax:
c = (b*f(a)-a*f(b)) / (f(a)-f(b))
if f(c) == 0 or abs(a-b) < Tol:
break
if (np.sign(f(c)) == np.sign(f(a))):
a = c
else:
b =c
N += 1
return c
这是我正在使用的两个功能。 M2呼叫G2。
def G2(T):
return (48*np.pi**2)/(11*N*(np.log(lam*T+Ts))**2)
def M2(T,func):
return ((N/6)*(func(T))*(T**2))**0.5
M2有效。我保证。
我的问题是我不知道如何推广Regula来调用M2,后者调用G2。
请帮忙!
Min_Mass = Regula (M2(T,G2), Tol, 50, a, b)
a = 1
b = 4
T = np.linspace(a,b,100)
Tol=1e-13
答案 0 :(得分:0)
为什么不简单地使用:
...
def M2(T,value):
return ((N/6)*( value )*(T**2))**0.5
Min_Mass = Regula (M2(T,G2(T)), Tol, 50, a, b)
...