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
为什么不简单地使用:
相关问题 更多 >
编程相关推荐