Python-сообщество

Требуется найти корни уравнения вида :

(E^(200 I k0) - 3 E^(600 I k0))/(-3 + E^(400 I k0)) = 1

In[33]:= NSolve[(E^(200 I k0) - 3 E^(600 I k0))/(-3 + E^(400 I k0)) == 1 
&& -0.05067731239 <= Re[ k0 ] <= 0.05067731239, k0]
Out[33]={{k0 -> -0.0429185 + 0. I}, {k0 -> -0.0314159}, {k0 -> -0.0199133 +  0. I}, 
{k0 -> -0.0115026 + 0. I}, {k0 -> 0.}, {k0 -> 0.0115026 + 0. I}, {k0 -> 0.0199133 + 0. I}, 
{k0 -> 0.0314159}, {k0 -> 0.0429185 + 0. I}}

 from sympy import Symbol, sin, cos, I, exp, pi, solve, Eq, simplify, nsolve, solveset, pprint, ImageSet
k0 = Symbol('k0')
fun = (exp(200*I*k0) - 3*exp(600*I*k0))/(exp(400*I*k0) - 3) - 1
froot = solveset(fun, k0)
pprint(froot)

Отредактировано Simarglur (Ноя. 11, 2018 20:34:07)

ну ответ очевиден.
Код вольфрама не соответствует коду на питоне.
Немного почитайте докупи допишите правильно и все получится.
Но не угадывайте%)

Simarglur
Может вы знаете чем тут можно помочь?

  
>>> solveset((exp(200 * x * I) - 3 * exp(600 * x * I)) / (-3 + exp(400 * x * I)) - 1)
повисло

Отредактировано py.user.next (Ноя. 12, 2018 00:16:57)

py.user.next
На Wolfram вычисляет корни.
На simpy виснет

Отредактировано Simarglur (Ноя. 12, 2018 00:14:07)

Simarglur
вот и я о том же

Отредактировано py.user.next (Ноя. 12, 2018 00:23:10)

py.user.next
Можешь запулить

Отредактировано Simarglur (Ноя. 12, 2018 00:34:41)

А какие-то другие способы есть? просто по всякому пробовала ни в какую. Может есть еще что-то чего я не знаю ?

Simarglur
А какие-то другие способы есть?

Отредактировано py.user.next (Ноя. 12, 2018 00:40:47)

py.user.next
Есть клиент для Wolfram

пока выкрутилась тем что ищу через nsolve в окрестностях каждой точки, которая выбирается в заданном интервале с определенным шагом, а потом упорядочиваю результаты. Неудобно, но что поделать.

Отредактировано Simarglur (Ноя. 12, 2018 01:44:46)