5．Pythonのプログラム 5.1 プログラムの流れ 図1 各記号の定義と意味はnoteの別の記事[1]を参照 5.2 プログラムのポイント 5.2.1  scipyモジュールを使う 特殊関数である球ベッセル関数$${j_n(z)}$$や球ノイマン関数$${n_n(z)}$$が計算できる関数が収められています。Pythonを使うメリットのひとつです ...

I would expect linspace(X, X, N) to produce an array with N copies of X. This works for finite values, but not for ±inf. This is an edge case that probably doesn't affect many people, but I don't ...

Python programmers usually strives to have one obvious way. But as soon as you discover all three methods, choosing one of them becomes weird (especially for non-interactive/long-term code).