素数階段ならぬナベアツ階段(ナベアツ数 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年2月14日 Microwave Journal今月号で透明アンテナ、6Gタワー、月面用アンテナなど、Keysightが6GでEUと提携、TDKが日本酒の味わいを視覚化&車載PoC用巻き線インダクタ、QualcommがSnapdragon 6 Gen 4 Mobile発表(2025.02.14)
- 高周波・RFニュース 2025年2月13日 iFixitがSamsung Galaxy S25 ultraを分解、ぐりぐり回せる3D CTモデルに、Spirentの5G Outlook レポート、QualcommのオンデバイスAIホワイトペーパー、InfineonがCoolSiC™ MOSFETs 650 Vの Q-DPAKとTOLLパッケージで発表(2025.02.13)
- 高周波・RFニュース 2025年2月12日 GSAが5G standaloneウェビナー開催、CTIAと5G americasがケーブルに対抗してSpectrum for Broadband Competitionで連携、NXPがエッジAIのKinaraを買収、PBKreviewsがSamsung Galaxy S25/plus/ultra全部分解していて5Gミリ波AiPは全部1個(2025.02.12)
« 新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数&ロジスティック関数&Log-Logプロットでべき関数フィッティングした(5/16更新) 日本も韓国も全然増加スピードが落ちない。イギリスもなんかちょっとずつ増えて行ってる。 | トップページ | 松屋で極上グラナパダーノ チーズトマト厚切りロースかつ定食(コロッケ付き、豚汁変更のフルセット)でいただく。チーズが粉になって風味がよくなった。 »
コメント