gnuplot - gnuplot によるミーゼス分布

Question

いくつかのデータをフォン ミーゼス分布に適合させる必要があります。

しかし、その前に、この関数を最初に gnuplot でプロットしたいと思います。さまざまな k 値の曲線を再現したいのですが、ウィキペディアで見ることができます: https://en.wikipedia.org/wiki/Von_Mises_distribution .

私の問題は、次数 0 の修正ベッセル関数である項 Io(k) に依存しています。

gnuplot では、関数で Jo(k) をプロットできます。

> besj0(k)

でもできないんです：

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

gnuplot でミーゼス分布をプロットする方法を知っている人はいますか?

Answer 1

Gnuplot の開発バージョンを使用できます。これは現在、第 1 種の修正ベッセル関数をサポートしています。具体的にI0は、 として利用できますbesi0。