
Я знаю, что существуют решения M или M-1 в зависимости от значения $\Delta$. Кто-нибудь знает какой-либо модуль или алгоритм, который я могу использовать в Python для этой задачи?
Я пробовал этот способ, но он очень чувствителен к допуску и не работает:
def f(k, N, d):
return np.tan(k*(N)) - d*np.sin(k)/(1+d*np.cos(k))
k = np.linspace(0, np.pi, 10000+1, endpoint=False)[1:]
def ksolutions(k,N,d):
solutions=[]
tol=4*1e-4
for i in k:
if abs(f(i,N,d))
Подробнее здесь: https://stackoverflow.com/questions/758 ... -in-python
Мобильная версия