時代に逆行してCOBOL(GnuCOBOL)を学んでみる(番外編） ChatGPTにCOBOLでナベアツ数（3の倍数と3がつくときだけアホになる）を作ってもらう。3がつくかどうかの判定に再帰を使っているのがやるな！と思った。

Fizzbuzzは絶対どこかにコードがあるからそれじゃないが似てるもの…あ！ナベアツ数にしよう。

ChatGPT無償版の方を使っています。こんな感じで聞くと。。。

おお！それっぽいのを出してきた。しかも再帰的に判定してたりする！

ただこのままでは動かない。。。何回かやり取りをして、最終的にできたコードはこちら。

（ちょっとだけ手直ししてますが）

そして結果。ちゃんと30からアホが続く！

ちなみにナベアツ数はこんなのもやってました。

素数階段ならぬナベアツ階段（ナベアツ数 Nab(x)をxの関数としてプロット）を擬ナベアツ関数(quasi-Nabeatsu function)と比較（加藤文元さんの1ページ論文を見て。Pythonで計算） 

一般化ナベアツ数Nab_n(x) (xまでのうちnの倍数とnがつく数のときだけアホになる）のナベアツ数階段をPythonでプロットする。この前やった加藤文元さんのナベアツ数の図示の続き。どれもNab_n(x)/x→1？ 

 一般化ナベアツ率(一般化ナベアツ数Nab_n(x) (xまでのうちnの倍数とnがつく数のときだけアホになる）をxで割ったもの。xまでにいくらナベアツ数があるか？）をPython+Numbaで計算してプロット。

これまでのCOBOLシリーズ：

時代に逆行してCOBOL(GnuCOBOL)を学んでみる(1) Windows11のWSL2でインストール、Visual Studio Codeの拡張機能を使い、10進31桁まで計算できることを生かして二次方程式で桁落ちが落ちる例題をそのまま計算。言語そのものはそんな批判されるほどのこともない気がした。

時代に逆行してCOBOL(GnuCOBOL)を学んでみる(2)話題のMojoが例題としてマンデルブロ集合を使っているそうなので、COBOLでマンデルブロ集合を描いてみよう。ただしテキストベースで！ 多重ループの練習。

時代に逆行してCOBOL(GnuCOBOL)を学んでみる(3) ロジスティック写像の分岐図を31桁まで10進計算ができることを生かして描いてみる。テキストベースで！今回は2次元配列の練習。

時代に逆行してCOBOL(GnuCOBOL)を学んでみる(4) 10進31桁まで計算できることを生かして4段4次のルンゲクッタ法でローレンツ方程式を計算して図示してみる。テキストベース（アスキーアート）で！ 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(5)　中二病のような名前のブルースカイカタストロフィ(Blue-sky catastrophe)を生じるGavrilov Shilnikov modelをドルマン・プリンス法（ode45)で計算する。テキストベース（アスキーアート）で！

時代に逆行してCOBOL(GnuCOBOL)を学んでみる(6) 31桁の10進演算で複素数を計算するライブラリを作ってみる。外部ファイルのサブルーチンを呼び出すCALL文とUSINGの練習。

時代に逆行してCOBOL(GnuCOBOL)を学んでみる(7) 31桁の10進演算でLU分解で連立一次方程式を解くサブルーチンを作ってみる。Numerical Recipes（ニューメリカル・レシピ）のFORTRAN77(Cじゃなく）の移植版。

 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(8) 31桁の10進演算でFFT（高速フーリエ変換 Fast Fourier Transform ）サブルーチンを作ってみる。Numerical Recipes（ニューメリカル・レシピ）のFORTRAN77(Cじゃなく）の移植版。matlabと同じ規格化で例題はNumPyから。