素数階段ならぬナベアツ階段(ナベアツ数 Nab(x)をxの関数としてプロット)を擬ナベアツ関数(quasi-Nabeatsu function)と比較(加藤文元さんの1ページ論文を見て。Pythonで計算)
加藤文元さんがナベアツ数の解析をされているのを初めて知った。
これですね。 pic.twitter.com/oSSUsQTaqs
— Fumiharu Kato 加藤文元(Bungen) (@FumiharuKato) May 15, 2021
おお、これは素数階段みたいだ。ということでPythonで描いてみよう。ついでに最後にquasi-Nabeatu function (nab(x))として記載されているものとも比較しよう。とりあえず10^5くらいで数分かかる。
Pythonのコードはこちら:
import matplotlib.pyplot as plt
import numpy as np
x=100000
count = 0
n = range(1, x+1)
Nab = []
qNab = []
for k in range(1, x+1):
if k % 3 == 0:
count += 1
else:
if '3' in str(k):
count += 1
Nab.append(count)
q = k - (2/3)*k**np.log10(9) - 1
qNab.append(q)
plt.figure(figsize=(12,8))
plt.bar(n, Nab, label = 'Nab(x)')
plt.plot(n, qNab, label='quasi-Nab(x)', color = 'Red')
plt.grid()
plt.title("Nabeatsu Numbers Nab(x)", fontsize=24)
plt.xlabel('x', fontsize = 18)
plt.ylabel('Nab(x)', fontsize = 18)
plt.tick_params(labelsize=18)
plt.legend(fontsize=18)
plt.xscale("log")
plt.yscale("log")
plt.show()
こちらも参照:
« 新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数&ロジスティック関数&Log-Logプロットでべき関数フィッティングした(5/16更新) 日本も韓国も全然増加スピードが落ちない。イギリスもなんかちょっとずつ増えて行ってる。 | トップページ | 松屋で極上グラナパダーノ チーズトマト厚切りロースかつ定食(コロッケ付き、豚汁変更のフルセット)でいただく。チーズが粉になって風味がよくなった。 »
「日記・コラム・つぶやき」カテゴリの記事
- 高周波・RFニュース 2025年5月22日 三星電機(SEMCO)が165℃対応の車載インダクタ発表、KYOCERA AVXがリップル電流についての技術文書発行、QualcommとXiaomiの契約15年目、OmdiaがNokiaをPrivate 5Gの2025の王者と決定(2025.05.22)
- Google ColabのJulia言語でFPUT問題(Fermi–Pasta–Ulam–Tsingou、非線形結合した振動子が最初に与えたモードに戻る再帰現象)をDifferentialEquations.jlの2階の常微分方程式ソルバーDynamicalODEProblemでシンプレクティック8次のKahanLi8で計算、振動子の動きも動画にしてみる。(2025.05.22)
- 高周波・RFニュース 2025年5月21日 TDKが0201のRFインダクタ発表、InfineonがUWBのFiraコンソーシアムの理事会に、ubloxがロボット用GNSSモジュール発表、FibocomがMediaTekのT930を使った5Gモジュール発表、Motolora Edge 60 Pro分解動画(2025.05.21)
« 新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数&ロジスティック関数&Log-Logプロットでべき関数フィッティングした(5/16更新) 日本も韓国も全然増加スピードが落ちない。イギリスもなんかちょっとずつ増えて行ってる。 | トップページ | 松屋で極上グラナパダーノ チーズトマト厚切りロースかつ定食(コロッケ付き、豚汁変更のフルセット)でいただく。チーズが粉になって風味がよくなった。 »
コメント