CP/M-80 (FN) ファイル名は大文字のこと。拡張子はFOR。 cpm.exe f80 (FN), (FN)=(FN) cpm.exe L80 (FN)/E, (FN)/N cpm.exe (FN) C Newton's Method PROGRAM SQROOT WRITE(*, *) 'input positive real' READ(*, *) A X = 0.1 * A DO 10 I = 1, 10 X = 0.5 * …

ブルーバックスの芹沢正三著『素数入門』6-5節に、古代バビロニアの三角関数表が載っているので、それを計算してみます。 数学の表現形式として、古代エジプト文明は測量術から幾何学を、古代メソポタミア文明は天体観測から整数論を好んだようです。古代ギ…

散乱に関するラザフォードの公式では、クーロン場の符号（斥力・引力）に依らずに同じ角度依存性が得られるわけですが、それについて計算結果を見てみることにします。 対称性に鑑みて、散乱体を含む二次元平面で切った計算をすることにします。散乱ポテンシ…

IDE の ドライブに ATA・ATAPI 制御コマンドを送って、シリアル番号やファームウェアバージョンなどのハードウェア固有の情報を得ることができます。 ■実行結果 ■Fortranソースコード 参考にしたページは、主にここ http://www.usefullcode.net/2007/02/hdd.…

原子核の存在を示唆することになった、ラザフォードによるα粒子散乱実験の、モデル計算を試みます。（ゆくゆくは散乱ポテンシャルを変えて散乱断面積の理論値を散乱体の影を見て実感したいという、おぼろげな希望が背景にあります。）固定されたクーロン型の…

コンウェイのライフゲームというものがあります。昔、マイコン上でよく遊ばれていたものです。あれは２次元でのものでしたが、OpenGL があるのでせっかくだから三次元でやってみることにしました。結論から言うと、生命発生・維持・死亡のバランスが難しく、…

以前の結果が間違っていたので訂正します。 http://d.hatena.ne.jp/fortran66/20090214 glutSolidCube を用いると、座標原点が動かされるのか、影を作る射影行列が正しく働かないようです。 そこで glVertex3f を用いて描画するようにしました。 射影面の法…

FORTRAN66 のプログラムを書いたことはあまりないんですが、制限がきつくて思いがけないところでしかられます。とりあえず、台形公式とシンプソンの公式でを求めるという典型的例題を解いてみました。IBM の16進浮動小数のためか単精度では意外に精度が出な…

JOB CARD //FORTHCLG JOB (001),'FORTRAN H FIBONACCI', // CLASS=A,MSGCLASS=A,MSGLEVEL=(1,1) //FIBONACC EXEC FORTHCLG,REGION.FORT=384K //FORT.SYSLIN DD UNIT=SYSDA //FORT.SYSABEND DD SYSOUT=A //FORT.SYSIN DD * C C FIBONACCI C READ(35, 100) N 1…

パソコンの興隆は FORTRAN77 の普及期に重なるのでパソコン用の信頼できる FORTRAN66 コンパイラは無いように思います。FORTRAN66 はコンピュータ言語初の標準化工業規格で、IBM の FORTRAN IV を元に制定されました。IBM の FORTRAN IV が使えればそれがも…

DEC 社の Visual Fortran (DVF) Ver.5 にひきつづき DVF Ver.6 が Win95/98 および WindowsNT 上の Fortran95 規格対応コンパイラとして出荷されました。MS Visual C++ Ver.6 と統合できるのが売りのひとつでしたが、DLL のバージョンがらみの微妙な問題が生…

OpenGL の範囲内で ARToolKit の像に影を付けることを目標とします。適当にググったところ、以下のページに行き当たりました。簡潔な説明で非常に参考になります。単純な幾何学での導出がされていて理解しやすいです。実際のプログラム例もあって助かりまし…

DEC社 の Windows 95/NT4 対応の Fortran90 コンパイラ。Microsoft は Fortran 言語開発から撤退するにあたり、引継ぎを DEC に任せました。バージョンが５から始まるのは、MS-FPS4 の後を継いだためだと思われます。MS-FPS4 からは優待割引の upgrade path …

Windows95 で動作する Fortran90 対応のコンパイラでした。バージョンが 4.0 に飛んだのは、Visual C++ と開発環境統合が統合できたので Visual C++ のバージョンに合わせたからだったと記憶しています。Win32 API 呼び出しなども出来るようになっており、そ…

Windows3.1 上で動作する FORTRAN77 対応のコンパイラでした。DOS-EXTENDER を用いて 32 bit モードで動作するので、大型計算機用のソースをそのまま持ってきても動かせて感激でした。とはいえ、コンパイラには微妙なバグがあって色々困りました。非公式に 1…

Fortran90 風に書いたもの。 PROGRAM test1 USE mpi IMPLICIT NONE INTEGER :: istatus(MPI_STATUS_SIZE) INTEGER :: ierr, idest, itag = 0 INTEGER :: irank, nrank, isource, mess CHARACTER(LEN = 40) :: message CALL MPI_INIT( ierr ) CALL MPI_COMM_R…

しばらくぶりに、ARToolKitをいじってみます。 ローレンツ・アトラクターは複雑な三次元図形なので、ARToolKit で任意方向から眺め回せれば愉快かなとふと思う。去年作ったプログラムの描画部分を数行ばかり変えることで実現できるはず。アトラクター部分は…

MS-FORTRAN v5.1 (PC-9801 用)。FORTRAN77 規格に（ほぼ？）完全対応していました。マニュアルには数値フォーマット等について詳しく書いてあって感心しました。とはいえ当時定価が9万8000円だったのだから、そのくらい当然か。DOS の 640K の制約のせいで、…

昔学んだテキストが出てきたので、f90 で書き直してみたりしなかったりします。http://h50146.www5.hp.com/solutions/hpc/stc/soft/mpi.html MPIプログラミング入門 http://h50146.www5.hp.com/solutions/hpc/stc/soft/pdfs/mpi_training.pdf （PDF直リン） …

Fortran で、複数個の整数変数の AND や OR を取りたいことがあるのですが、Fortran90 以降で用意された IAND や IOR 関数は引数を二つしか取れず、IAND や IOR を入れ子式に用いなければなりません。例 m = IOR(I, IOR(J, K)) これを m = IOR(I, J, K) の様…

ネタ元 http://www.dolfyn.net/dolfyn/f03gl_en.htmlいわゆる Green Book からのサンプルのつづき。 元プログラムからの修正点は、RESHAPE ルーチンの引数の value 属性を消したこと。 module fgl08m use opengl_gl integer(kind=GLint) :: winWidth = 300, …

ネタ元 http://www.dolfyn.net/dolfyn/f03gl_en.htmlいわゆる Green Book からのサンプルのつづき。 元プログラムからの修正点は、RESHAPE ルーチンの引数の value 属性を消したこと。C Binding の削除。 program fgl06 use opengl_gl use opengl_glu use op…

ネタ元 http://www.dolfyn.net/dolfyn/f03gl_en.htmlいわゆる Green Book からのサンプルのつづき。 元プログラムからの修正点は、RESHAPE ルーチンの引数の value 属性を消したこと。C binding を消したこと。 program fgl04 use opengl_gl use opengl_glu …

森の妖精ことビリー・ヘリントン氏の来日がアナウンスされました！ http://anime-ch.nicovideo.jp/static/wonderful_hobby9 楽しみですね。 ネタ元 http://www.dolfyn.net/dolfyn/f03gl_en.htmlいわゆる Green Book からのサンプルのつづき。 元プログラムか…

ネタ元 http://www.dolfyn.net/dolfyn/f03gl_en.htmlいわゆる Green Book からのサンプルのつづき。 program fgl02 use opengl_gl use opengl_glu use opengl_glut external display integer(kind=GLint) :: iwin real(kind=GLdouble) :: window(4) = (/ 0.0…