用python求方程的零点

Question

我试着写一段代码，在这个等式中，找到n。你知道吗

其余的作为用户定义的变量。你知道吗

from scipy.optimize import fsolve
from scipy.stats import t

def f(alpha, beta, sigma, delta, eps):
    n = ((t.ppf(1-alpha,2*n-2) + t.ppf((1-beta)/2,2*n-2))**2*sigma**2)/(2*      (delta-abs(eps))**2)

我还希望能够设置不同的参数场景，然后让它输出一个参数表和结果（例如，输入alpha1，alpha2，beta1，beta2等，然后输出[alpha1，beta1，…，n]，[alpha1，beta2，…，n]）。如果有人能给我指出正确的方向，我不太确定最好的方法是什么。你知道吗

Answer 1

通过你的方程的外观，你试图找到满足统计检验方程的观测数（n）。如果是这样的话，那么n是自然数（0，1，2…等），很容易计算。你知道吗

你可以自己设置一个解算器，其中n是可数，而result是方程的“result”：

for n in range(0, 1000):
    result = your_function(n, other_parameters)

然后您只需通过设置以下条件来检查是否满足方程式：

    if n >= result:
        print "result:", n
        break  # This will exit the loop

对于测试不同的用户给定参数，可以设置另一个循环，循环alpha、beta等的不同值。你知道吗