首都大学東京物理学コースコミュのFortran

ログインして参加する

mixiユーザー 2006年05月27日 10:29

実験1のFortranのプログラムの質問が出たので、
プログラムの解説を載せましょうか？

> program picalc　　*プログラム名　picalc
> write(*,*) "How many random number?"　　　
> read(*,*) ntot　　*(多分、整数型で読み込まれてる)
> ncircle=0　　*変数(これも多分整数型)ncircleを0に初期化
> do i=1, ntot　　*iがntotになるまで繰り返し(1ずつ増やしていく)
> x=rand()　　*xに、発生した乱数を代入
> y=rand()　　*yに、発生した乱数を代入
> r=sprt(x**2+y**2)　　*r=√(x^2+y^2)（距離の計算）
> if(r .le. 1.0) ncircle=ncircle+1　　*rが1.0よりも小さいなら、ncircleを1増やす
> enddo　　*Doループ終了
> pi=float(ncircle)/float(ntot)*4　　*ncircleは整数型になってるので、
　　　　　　　　　　　　　　　　　　　　　　floatで実数型に直してから計算
　　　　　　　　　　　　　　　　　　　　　　じゃないと、円周率が整数で出てくる(+_+;)
> write(*,*) 'pi=', pi
> end

コメント(4)

最初
全て
最新の40件

1

[1] mixiユーザー 05月28日 19:51

ホントありがとうございます。

[2] mixiユーザー 05月29日 00:53

もう出しちゃったからちょっと遅いです（泣）
Ｃでもやってほしかったです～

[3] mixiユーザー 05月29日 23:46

つぶらんは遅かったかぁ～(^^;)
Ｃも分かるから、良かったらプログラムのっけて下さい(^^)b
他の実験でも、分かんないとこあったら、聞いてくださいなぁ～

ログインすると、みんなのコメントがもっと見れるよ

最初
全て
最新の40件

1

mixiユーザー: ログインしてコメントしよう！

首都大学東京物理学コース更新情報

最新のトピック
- 一覧へ
最新のイベント
- まだ何もありません
最新のアンケート
- まだ何もありません

首都大学東京物理学コースのメンバーはこんなコミュニティにも参加しています

星印の数は、共通して参加しているメンバーが多いほど増えます。

首都大学東京物理学コース

人気コミュニティランキング

困ったときには

関連ワード

mixiで趣味の話をしよう: 新規登録（無料）; ログイン

このページの上部へ