時代に逆行してCOBOL(GnuCOBOL)を学んでみる(6) 31桁の10進演算で複素数を計算するライブラリを作ってみる。外部ファイルのサブルーチンを呼び出すCALL文とUSINGの練習。
今回は初心に戻って複素数の計算ライブラリ。
前にExcel VBAとかJavaScriptとかで作っていたのでCOBOLでもやってみよう。
複素数は変数宣言でレベル番号を使って上位が複素数、下位が実部、虚部としておくとサブルーチンを呼び出すとき便利。構造体っぽい。
外部ファイルを5つ準備する（四則演算＋表示）。linkage section.でサブルーチンの引数を定義しておき、procedure divisionのusingでそれを使うことを宣言する。
メインから呼び出すには call "ファイル名" using 引数でOK。
四則演算の結果はこちら。0がちかちかするがCOBOLっぽいのであえて消してない。
これまでのシリーズ：
時代に逆行してCOBOL(GnuCOBOL)を学んでみる(3) ロジスティック写像の分岐図を31桁まで10進計算ができることを生かして描いてみる。テキストベースで！今回は2次元配列の練習。
時代に逆行してCOBOL(GnuCOBOL)を学んでみる(4) 10進31桁まで計算できることを生かして4段4次のルンゲクッタ法でローレンツ方程式を計算して図示してみる。テキストベース（アスキーアート）で！ 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(5) 中二病のような名前のブルースカイカタストロフィ(Blue-sky catastrophe)を生じるGavrilov Shilnikov modelをドルマン・プリンス法（ode45)で計算する。テキストベース（アスキーアート）で！
« 松屋で平成レトロ トンテキ定食3枚盛（ご飯特盛）をいただく。ものすごくニンニクが強く、肉も結構分厚くて量も十分。 | トップページ | 大阪中之島美術館に行ったとき、あちこちにアンテナがついたトランシーバーのようなものが置いてあった。温度センサ？と思ったらやっぱりHanwellのRL4007 無線温度ロガー（440MHz帯）だったよう。でもその下には紙の記録計も。両方使っているのかな？ »
「パソコン・インターネット」カテゴリの記事
- 出遅れましたがFizzBuzzをExcelのMAP,LAMBDA,SEQUENCE,IFS関数を使って一行（というかただの１セル入力）で作る。(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関数を使って一行（というかただの１セル入力）で作る。(2023.09.24)
- （速報続報）iPhone15 Proが早くも分解。USモデルなのでミリ波アンテナ３つが見えてる。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関数を使って一行（というかただの１セル入力）で作る。(2023.09.24)
- （速報続報）iPhone15 Proが早くも分解。USモデルなのでミリ波アンテナ３つが見えてる。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)
« 松屋で平成レトロ トンテキ定食3枚盛（ご飯特盛）をいただく。ものすごくニンニクが強く、肉も結構分厚くて量も十分。 | トップページ | 大阪中之島美術館に行ったとき、あちこちにアンテナがついたトランシーバーのようなものが置いてあった。温度センサ？と思ったらやっぱりHanwellのRL4007 無線温度ロガー（440MHz帯）だったよう。でもその下には紙の記録計も。両方使っているのかな？ »
コメント