Распределение фон Мизеса с gnuplot

Question

Мне нужно подогнать некоторые данные к распределению фон Мизеса:

но прежде чем сделать это, я хотел бы сначала написать эту функцию в gnuplot. Я хочу воспроизвести кривые для разных значений k, вы можете увидеть их в википедии: https://en.wikipedia.org/wiki/Von_Mises_distribution.

Моя проблема опирается на член Io(k), который является модифицированной функцией Бесселя порядка 0. Io(k)=Jo(ik), где Jo - функция Бесселя для n=0, в основном из этого уравнения: ,

В gnuplot можно построить Jo(k) с помощью функции

> besj0(k)

но я не могу:

besj0(i*k) when the argument is imaginary.

У кого-нибудь есть идея, как построить распределение фон Мизеса в gnuplot?

Answer 1

Вы можете использовать версию разработки Gnuplot, которая в настоящее время поддерживает модифицированные функции Бесселя первого рода. В частности, I0 доступно как besi0,