時代に逆行してCOBOL(GnuCOBOL)を学んでみる(2)話題のMojoが例題としてマンデルブロ集合を使っているそうなので、COBOLでマンデルブロ集合を描いてみよう。ただしテキストベースで! 多重ループの練習。
さてCOBOLをいまさら学ぶシリーズ第二弾。第一弾は二次方程式を桁落ちなしで計算することをやってみた。
今回は、、、そういやPythonより35000倍速いというMojoが例題にマンデルブロ集合をつかっていたということで!
マンデルブロ集合を描いてみよう。まあグラフ表示機能はないけど、、、大昔ラインプリンタで書いた人がいるので、
13世紀の僧侶はマンデルブロ集合を描かないが、1978年にマンデルブロより先にラインプリンタでマンデルブロ集合を描いた人たちがいる。
それに倣ってテキストベースで描く。
今回、ハマったのは多重ループを回すとき。COBOLは基本、文の最後にピリオド”."をつけるが、
多重ループの時は一番最後のend-performにだけ"."を付けないとだめということ。
結果がこちら。おお、ちゃんとマンデルブロ集合っぽい。
次は何をしようか…
※さらにフォントサイズを小さくして画面サイズを大きくするとこれくらいは描ける。
« 爆弾(呉勝浩さん)を読んだ。圧倒され、恐ろしいが面白い…些細な傷害事件で逮捕された卑屈な中年男が取調室で爆弾が爆発する、と予言し次々と本当になる。男の言動に振り回される刑事たち。人間の悪意とはと考えさせられ、そして衝撃の結末へ。2冠も納得。 | トップページ | やよい軒でGW付近限定のミックスグリル定食(サラダ付き)をいただく。デミハンバーグとチキンステーキ、ウインナーという豪華なもの。 »
「パソコン・インターネット」カテゴリの記事
- 出遅れましたがFizzBuzzをExcelのMAP,LAMBDA,SEQUENCE,IFS関数を使って一行(というかただの1セル入力)で作る。(2023.09.24)
- iOS17の新機能、Object CaptureをアプリHyper Captureで試す。ゴッホのスヌーピーとフェルメールの牛乳を注ぐミッフィーを3D スキャンでぐるぐる回す。3回角度変えてキャプチャが必要だった。プロセスはiPhone 12 Pro Maxなんで時間掛かるけどかなり綺麗に。(2023.09.21)
- 離散リアプノフ方程式 AXAᴴ - X + Q = 0がクロネッカー積とvecで計算できることを思い出した!(Vec Trickというそう)。Pythonのscipy.linalg.solve_discrete_lyapunovとnumpy.kronの両方で計算してちゃんと合うことを確認。(2023.09.22)
- Google検索の生成AI,SGEが使えるようになっていた。じゃあということでRIS(Reconfigurable Intelligent Surface)について聞いてみた。ChatGPT, BingのAIチャット、Bardと比較。ZTEが先日発表していたので聞いたが、SGEはその話にも触れた。詳しいのはChatGPTだった。(2023.09.02)
- PythonのmatplotlibにXKCD(棒人間の理系コミック、ホワット・イフのランドール・マンローさん作)の手書き風グラフが描ける機能があった!せっかくなのでルンゲクッタ法のDormand Prince8次(DOP853)でローレンツ方程式を描く。(2023.09.01)
「学問・資格」カテゴリの記事
- 浜村渚の計算ノート 10さつめ ラ・ラ・ラ・ラマヌジャンを読んだ。九章算術、ベクトル、四元数、電卓、そしてラマヌジャン!タクシー数も1+2+3+…=-1/12もいろんな公式も出てきます。カプレカー数も。高精度計算サイトkeisan.casio.jpにUPしているものとも関連していてよかった。(2023.09.26)
- 出遅れましたがFizzBuzzをExcelのMAP,LAMBDA,SEQUENCE,IFS関数を使って一行(というかただの1セル入力)で作る。(2023.09.24)
- (速報続報)iPhone15 Proが早くも分解。USモデルなのでミリ波アンテナ3つが見えてる。USモデルとその他でMLBの形から違う!(USはeSIM、他はSIMカード)、なのでスペースがなくて他モデルはミリ波アンテナ部分に部品乗せてる!(2023.09.23)
- 離散リアプノフ方程式 AXAᴴ - X + Q = 0がクロネッカー積とvecで計算できることを思い出した!(Vec Trickというそう)。Pythonのscipy.linalg.solve_discrete_lyapunovとnumpy.kronの両方で計算してちゃんと合うことを確認。(2023.09.22)
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2023年9月19日 Microwave Magazineの特集はRFIDや氷を検出する話、Microwave Journalで車載アンテナ評価でRanLOSというのを初めて知る、Gapwavesの多層導波管、PythonのRFライブラリScikit-RFに高木分解を使うTUG multiline TRLが。(2023.09.19)
「日記・コラム・つぶやき」カテゴリの記事
- 浜村渚の計算ノート 10さつめ ラ・ラ・ラ・ラマヌジャンを読んだ。九章算術、ベクトル、四元数、電卓、そしてラマヌジャン!タクシー数も1+2+3+…=-1/12もいろんな公式も出てきます。カプレカー数も。高精度計算サイトkeisan.casio.jpにUPしているものとも関連していてよかった。(2023.09.26)
- 出遅れましたがFizzBuzzをExcelのMAP,LAMBDA,SEQUENCE,IFS関数を使って一行(というかただの1セル入力)で作る。(2023.09.24)
- (速報続報)iPhone15 Proが早くも分解。USモデルなのでミリ波アンテナ3つが見えてる。USモデルとその他でMLBの形から違う!(USはeSIM、他はSIMカード)、なのでスペースがなくて他モデルはミリ波アンテナ部分に部品乗せてる!(2023.09.23)
- iOS17の新機能、Object CaptureをアプリHyper Captureで試す。ゴッホのスヌーピーとフェルメールの牛乳を注ぐミッフィーを3D スキャンでぐるぐる回す。3回角度変えてキャプチャが必要だった。プロセスはiPhone 12 Pro Maxなんで時間掛かるけどかなり綺麗に。(2023.09.21)
- 離散リアプノフ方程式 AXAᴴ - X + Q = 0がクロネッカー積とvecで計算できることを思い出した!(Vec Trickというそう)。Pythonのscipy.linalg.solve_discrete_lyapunovとnumpy.kronの両方で計算してちゃんと合うことを確認。(2023.09.22)
« 爆弾(呉勝浩さん)を読んだ。圧倒され、恐ろしいが面白い…些細な傷害事件で逮捕された卑屈な中年男が取調室で爆弾が爆発する、と予言し次々と本当になる。男の言動に振り回される刑事たち。人間の悪意とはと考えさせられ、そして衝撃の結末へ。2冠も納得。 | トップページ | やよい軒でGW付近限定のミックスグリル定食(サラダ付き)をいただく。デミハンバーグとチキンステーキ、ウインナーという豪華なもの。 »
コメント