くだすれFORTRAN（超初心者用）

1 名前：デフォルトの名無しさん [2006/01/24(火) 09:48:23 ] このスレッドは、他のスレッドでは書き込めない超低レベル、  もしくは質問者自身何が何だが分からない質問を勇気を持って書き込むスレッドです。  FORTRAN使いが優しくコメントを返しますが、  お礼はFORTRANの布教と初心者の救済をお願いします。  357 名前：デフォルトの名無しさん mailto:sage [2006/05/31(水) 19:49:52 ] >>356 program parking read(*,*) hour,fmin write(*,*) ifee(hour,fmin) stop end function ifee(hour,fmin) totalmin=hour*60.+fmin if(totalmin.le.60.) then ifee=300*int(totalmin/30) if(totalmin-30.*int(totalmin)/30.gt.0.) ifee=ifee+300 else ifee=600+200*int(totalmin-60.)/30 if(totalmin-30.*int(totalmin-60.)/30.gt.0.) ifee=ifee+200 endif return end 358 名前：デフォルトの名無しさん mailto:sage [2006/05/31(水) 21:09:21 ] >>357 ありがとうございます 359 名前：デフォルトの名無しさん mailto:sage [2006/05/31(水) 21:18:09 ] あ、忘れてたのですが最高課税額が4000円になるようにしてもらえませんか 本当に申し訳ない＞＜ 360 名前：デフォルトの名無しさん mailto:sage [2006/05/31(水) 22:26:55 ] >>359 357 じゃないけど。 せっかくプログラムがあるんだから、 いろいろいじってみようよ。 ヒント：if 文　をどこかに一行いれる。 361 名前：357 mailto:sage [2006/06/01(木) 11:32:39 ] 慌てて作ったら間違っていたみたい。 それを修正するのも良い勉強になるでしょう。 あと、変数に7文字以上使っているのも問題かも。 実は厳密なFORTRAN77を知らない。 362 名前：デフォルトの名無しさん [2006/06/01(木) 17:35:24 ] >>361 それ以前に小文字を使っているぞいｗ 363 名前：357 mailto:sage [2006/06/01(木) 20:21:52 ] そういやそうだったねorz。字下げもしない方が良かったかな。 364 名前：デフォルトの名無しさん [2006/06/01(木) 23:33:12 ] 1/√π×exp（-ｘ^2）は中央値0、偏差1の正規分布曲線を与える。 〔-1、1〕および〔-2、2〕の範囲に存在する確率はいくらか。台形法で求めよ。 学校の課題なのですがわかる方いたらお願いします。 先生いきなり難しいよ・・・ 365 名前：デフォルトの名無しさん mailto:sage [2006/06/01(木) 23:47:33 ] 小学生の算数じゃなかったか？＞台形面積 366 名前：デフォルトの名無しさん [2006/06/02(金) 01:38:52 ] >>363 いや字下げは文法にかなっているからおｋ FORTRAN７７までは大文字のみが文法で規定されていた。 小文字は見にくくてｷﾓｽｗ >>365 おまえは>>1を読み返せ！ このスレは糞既出FAQにも丁寧にマジレスし、宿題も文句も言わずに解いてやり、 FORTRAN依存症を作り出して初心者をFORTRAN Userに引き込む秘密結社の陰謀だ。 まぁ本スレの方ならなじった上で、こっちのスレに誘導してくれればおｋｗｗｗ 367 名前：1/2 mailto:sage [2006/06/02(金) 02:31:16 ] >>364　Fortran90で難しく書いてみた。が、少しイモくさいｗ MODULE mod_integ IMPLICIT NONE REAL :: pi CONTAINS REAL FUNCTION trapez(n, xmin, xmax) IMPLICIT NONE INTEGER, INTENT(IN) :: n REAL , INTENT(IN) :: xmin, xmax INTEGER :: i REAL :: h, x h = (xmax - xmin) / n trapez = 0.0 DO i = 1, n x = xmin + (i - 1) * h trapez = trapez + sub_trapez(x, h) END DO RETURN END FUNCTION trapez 368 名前：2/2 [2006/06/02(金) 02:32:02 ] REAL FUNCTION sub_trapez(x, h) IMPLICIT NONE REAL, INTENT(IN) :: x, h sub_trapez = 0.5 * h * ( func(x) + func(x + h) ) RETURN CONTAINS REAL FUNCTION func(x) IMPLICIT NONE REAL, INTENT(IN) :: x func = EXP( -x**2 ) / SQRT(pi) RETURN END FUNCTION func END FUNCTION sub_trapez END MODULE mod_integ ! PROGRAM oppai USE mod_integ IMPLICIT NONE INTEGER :: i REAL :: s, s1, s2 pi = 4.0 * ATAN(1.0) ! arctan(1.0) = pi / 4; tan(pi/4) = 1.0 ! [-1..1] s1 = trapez(100, -1.0, 0.0) ! [-2..2] s2 = trapez(200, -2.0, 2.0) ! [-inf..+inf] s = trapez(1000, -10.0, 10.0) PRINT *, s, s1, s2 STOP END PROGRAM oppai >>367-368 369 名前：367 [2006/06/02(金) 02:53:47 ] 実行結果 0.9999998 0.4213469 0.9953194 Press any key to continue 微妙に問題の偏差の定義などが一般的な統計の場合とずれているが、 ここでの定義を信じてゆくことにするｗ 正規分布の−∞〜＋∞の積分は１になるように規格化してあるはず。 （√π）で割っている部分。それを確かめたのが変数ｓ。 ただし積分範囲は−１０〜＋１０までにしている。 無限大じゃないじゃんか！というかもしれないが、これでも有効数字の範囲では もう無限大と変わりない。というか、関数は０になっている。 誤差は台形公式の誤差だろう。分割数ｎ（分割の間隔）を変えると調べられる。 -1..1の積分は約０．４２１、-2..2までの積分は約０．９９５になっている。 つまり、−∞から＋∞まで積分すると1.0になるはずものだが、 −２〜＋２までの範囲で約９９．５％が尽くされることが分かる。 プログラムでは、メインルーチンは積分範囲を指定して積分ルーチンを呼んでいる。 サブルーチンはモジュールの中に入れてある。 台形（ｔｒａｐｅｚｏｉｄ）積分は、等間隔にｘ軸を切ってそれを台形の高さとし、それぞれの 小さな台形の面積の和を返している。小さな台形の面先はさらに台形面積を返す ルーチンを呼び出している。小さな台形の面積を返すルーチンは正規分布関数の 値を返すルーチンを呼んでいる。 ・プログラム 　・積分ルーチン 　　　・台形面積ルーチン 　　　　　・正規関数ルーチン FORTRAN77で書いてほしければ、はじめから遠慮せずに言うようにｗ 370 名前：367 mailto:sage [2006/06/02(金) 02:55:02 ] なお、間違っていても謝罪も賠償もしないので４６４９！ 371 名前：364 [2006/06/02(金) 03:07:46 ] ありがとうございます！77とか90とかよくわかりません。学校のパソコンはどっちだろう・・・ これをテラパッドに入力してシグインで実行するわけですな 372 名前：367 [2006/06/02(金) 11:19:01 ] >>371 ! [-1..1] s1 = trapez(100, -1.0, 0.0) これまちがいｗ　ごめｗ　でも賠償は（ｒｙ 正しくは、s1 = trapez(100, -1.0, 1.0)　面積は０．８４２６７２９ 積分範囲がまちがって-1.0〜0.0になってたｗｗｗｗｗｗｗ なんか値がおかしいと思ったのだがｗ　標準偏差の定義が標準と異なるので 表で値を確かめるときに換算を間違えておｋと思ってしまったｗｗｗ ま、宴会から帰ってきた夜中だったので許せｗ ＊やや無関係な話題ｗ πの値はパラメータで与えてもいいが、ここでは４．０＊ＡＴＡＮ（１．０）を使った。 これは昔からよく使われている定石だが、どこかで１回計算しておかないと いけない。それがめんどい。だがこうしておくと精度を機械任せにできる。 今のプログラムではやや不自然で、結果醜くなっている。 Ｆｏｒｔｒａｎ２００３ならパラメータ文の中で関数が使えるようになるのだが・・・・ （パラメータは実行時ではなくコンパイル時に評価されるので、 今のところ四則演算とベキぐらいしか使えない） ＊正規分布の−∞から＋∞への積分 ∫exp(-x^2)dxは√πになるのだが、不定積分があらわせないので、 必死に積分しようとしないようにｗ 多重積分のテクニックで直交座標を円座標にして出すのが伝統芸。 ∫exp(-x^2)dx∫exp(-y^2)dy=∬exp(-(x^2+y^2))dxdy=∬exp(-r^2)ｒdrdθ =-1/2[exp(-r^2)]^∞_0 * 2π　=　1/2*2π = π つまり∫exp(-x^2)dx=√πだお。　 数値積分では無限大の積分は、普通変数変換して有限範囲にしてから実行する。 ここでの-10~10の積分は、あくまで確認のためのお遊びということで。 373 名前：364 [2006/06/02(金) 19:35:00 ] >>372 ありがとう。使わせてもらいます。 DAIKEIHOU gakusekibangou NAME sekibun.f conded at 29 May DIMENSION NJUN(J) DATA NJUN/10,50,100/ A=-1,0 B=1,0 WRITE(6,100)A,B FORMAT(5x,'SEKIBUN FOR 1/SQR(PI)*EXP(-X*X) A' ,F5,2,3X,'B',F5,2) DO 10 K=1,3 N=NJUN(K) H=(B-A)/N S=FX(A)+FX(B)/2,0 DO 20 J* J=1,N-1 X=A+BJ S*S+FX(X) ちなみにこれが学校の先生が見本としてやってたのを急いで写したものですが（途中まで） あってますか？ 374 名前：デフォルトの名無しさん [2006/06/02(金) 19:58:25 ] 微妙に写し間違えているｗ DIMENSION NJUN(3) ... S=( FX(A) + FX(B) ) / 2.0 DO 20 J = 1, N - 1 X = A + H * J S = S + FX(X) 20 CONTINUE S = S * H WRITE(6, *) 'BUNKATSU = ', NJUN(I), 'MENSEKI =', S 10 CONTINUE STOP END こんな感じになるのではないかな． 台形公式を>>364の様に素直に小さな台形の和を足して求めると， 関数の値を２回同じ値に対して求めなければならない． そこを改良すると，ここでの例題のように，両端だけを半分にして， 残りを一回計算して足し上げてゆけば良いことになる． 絵を描いて考えればすぐ分かると思う． NJUNは分割数だろう．何個の台形に分けるかを指定している． 数学的には分割が細かいほど正確な値になるが（積分は分割→０の極限で定義 されている）現実には数値誤差や計算時間の関係からほどほどの分割数がいい． いくつがいいかは時と場合によるので一概には言えない． 375 名前：デフォルトの名無しさん [2006/06/02(金) 22:03:08 ] 台形公式の改良版としてロンバーグ加速ってのがある。 台形公式の幅を半分ずつにしていって結果がどのように変化するかを調べて、 幅無限小の時の値を補外するというもの。 関数の形にもよるけど、台形公式とほぼ同じ処理量で凄い精度が出る。 376 名前：デフォルトの名無しさん [2006/06/04(日) 04:09:37 ] >>50 377 名前：デフォルトの名無しさん mailto:sage [2006/06/04(日) 04:22:25 ] ポインタが無いぶん最適化はしやすいと聞くが、 正直微々たる差だろうな。 378 名前：デフォルトの名無しさん mailto:sage [2006/06/04(日) 08:49:02 ] GCCの場合、どちらも一旦中間コードに直されてから最適化されるので、同じ構造のプログラムを書いた場合 基本的に最適化に差はほとんどない。 379 名前：デフォルトの名無しさん mailto:sage [2006/06/04(日) 09:27:59 ] >>378 ＧＣＣに限らず最近のコンパイラは、中間言語を出力する字句解析のフロントエンド部と、 機械語変換のバックエンド部に分離しているから、その理屈ではどんな言語でも 似たようなスピードになる。実際そういう傾向はある。 問題は今後のハードの発展が、ＣＰＵの並列化や、非連続な分散メモリー、 ＳＳＥのようなベクトル処理の方向に向かっていることで、言語自体が向いていないと 最適化がかかりにくくなってくる。 ＦＯＲＴＲＡＮの規格はハードウェアの進化のフィードバックを改定を受けているが、 Ｃ言語などは最近のネットセキュリティーの重視を受けてセキュアーさのフィードバックを 強く受けているようなので長期的にはどうなるか・・・ 最近のマルチコア・マルチスレッド・計算用専用プロセッサの利用など（ＧＰＵ流用や ＣＥＬＬとか）を見ていると、コンパイラによる自動並列化というより、ＭＰＩのような ＡＰＩ呼び出しで解決しようとしているようだし、どうなることかｗ 380 名前：デフォルトの名無しさん [2006/06/05(月) 01:24:19 ] すみませんフォートランの始め方からしてわかりません・・・ www.g.dendai.ac.jp/lecture/fortran/fort01.html このサイトにはFTN77 Personal Editionをダウンロードしろとあるのですが、 77でなく95しか置いてないみたいなんです。 Cpadが77向けなので95と77では合わずに不具合が生じてしまう気がするのですが・・・ どなたかご回答願います・・・ 381 名前：デフォルトの名無しさん [2006/06/05(月) 01:50:28 ] >>380 FORTRAN III pc8.2ch.net/test/read.cgi/tech/1104724162/ 861 名前：デフォルトの名無しさん[] 投稿日：2006/05/27(土) 21:37:28 FTN77 Personal Edition ってもうダウンロードできないの？ SalfordのHPやSilverfrostのHPからもリンクが消えてるんだけど… 862 名前：デフォルトの名無しさん[] 投稿日：2006/05/29(月) 01:52:14 >>861 FTN95がＦｒｅｅになったから消したんでないのかな？ 863 名前：デフォルトの名無しさん[] 投稿日：2006/06/03(土) 02:31:53 もうどこでもダウンロードできないんですかね? 864 名前：デフォルトの名無しさん[] 投稿日：2006/06/04(日) 00:35:52 >>863 ロシアか中国あたりにうｐられていそうな気もするがｗ 漏れは使ったこと無いけど、FTN77はそんなにいいものなのかい？ フリーのFTN95はバナーがでてウザかった気はするが使えたよ。 F77だったらOPEN WATFORやｇ７７などがあるので代替手段が無いわけではない。 WATFORもDOS EXTENDER時代に本来の開発が止まったものなので、昔懐かしい 感じを保っているはずだが・・・ｗ 865 名前：デフォルトの名無しさん[sage] 投稿日：2006/06/04(日) 01:38:40 「FTN77PE.EXE」でぐぐるとそれらしきものがいくつか見つかるが どれもなんか怪しい雰囲気がｗ >>864 単体で使う分には取り立ててよいというほどではないと思うけど CPad for Salford FTN77と一緒に使うとお手軽にプログラムの 作成・実行が出来るので重宝しています。 このセットで紹介しているHPも結構あるみたいだし。 382 名前：３８０ [2006/06/05(月) 01:58:19 ] およよ・・・ありがとうございますです。 がんばってみます。 すでにめげそう；； 383 名前：デフォルトの名無しさん [2006/06/05(月) 02:00:48 ] >>380 文法的にはFTN95はFTN77を包括している。 （完全に規約に準拠していると、７７の機能の一部が削られていることになるが） Cpadの設定を自分で書き直せば何とかなるだろう。 まぁ、そろそろ９０、９５に移行することをお勧めする。 といっても色々事情もあるだろうから、好きになさいｗｗｗ 安易な方法としては、違法くさいのがうｐされているので、具具ってそれを落としてくれば とりあえずは足りるべｗ　　 いつまであるかはわからんがｗ 384 名前：３８０ [2006/06/05(月) 02:04:19 ] みなさまありがとうございます。 がんばりますです。 385 名前：デフォルトの名無しさん [2006/06/05(月) 02:06:54 ] 95と77の違いってどんなところ？？ 386 名前：デフォルトの名無しさん [2006/06/05(月) 02:18:54 ] >>385 ９５はモジュール言語の特徴を大きく取り入れた。 ほとんど温泉旅館の新館と旧館のようになかば違う言語がくっついた感じになっている。 Fortran90の時点ではFORTRAN77を完全に包括していた。 Fortran95はFortran90のBugFix的な位置づけだが、一部Fortran90で廃止勧告が 出ていた機能を削っている。ただし、ほとんどのベンダーは過去互換性を重視して そういう機能を残している。 古い人なら、FORTRAN77がpascalと合体したと思えばいい。 387 名前：デフォルトの名無しさん [2006/06/05(月) 02:27:30 ] どもっす 388 名前：デフォルトの名無しさん [2006/06/05(月) 15:54:04 ] 行列がわかんねーorz 389 名前：デフォルトの名無しさん mailto:sage [2006/06/05(月) 16:47:31 ] >>388 それは FORTRAN の問題じゃなくて数学の問題ではなかろうか。 390 名前：デフォルトの名無しさん [2006/06/05(月) 22:07:14 ] DAIKEIHOU gakusekibangou NAME sekibun.f conded at 29 May DIMENSION NJUN(J) DATA NJUN/10,50,100/ A=-1,0 B=1,0 WRITE(6,100)A,B FORMAT(5x,'SEKIBUN FOR 1/SQR(PI)*EXP(-X*X) A' ,F5,2,3X,'B',F5,2) DO 10 K=1,3 N=NJUN(K) H=(B-A)/N S=FX(A)+FX(B)/2,0 DO 20 J* J=1,N-1 X=A+BJ S*S+FX(X) DIMENSION NJUN(3) ... S=( FX(A) + FX(B) ) / 2.0 DO 20 J = 1, N - 1 X = A + H * J S = S + FX(X) 20 CONTINUE S = S * H WRITE(6, *) 'BUNKATSU = ', NJUN(I), 'MENSEKI =', S 10 CONTINUE STOP END これ実行できないです。どこか間違ってますか？ 391 名前：デフォルトの名無しさん mailto:sage [2006/06/05(月) 22:26:40 ] >>390 コメントは行の先頭に"C"を付けろとか 小数点(".")とカンマ(",")を混同してるとか 文番号のないFORMATなんて何に使うんだとか 勝手に行を折り返すなとか 実行文の途中に配列宣言入れるなとか ツッコミどころ多数でどこもかしこも間違っている。 392 名前：デフォルトの名無しさん [2006/06/05(月) 23:39:37 ] >>390　もう少し自力でがんばれｗ PROGRAM DAIKEIHOU C gakusekibangou NAME C sekibun.f conded at 29 May DIMENSION NJUN(3) DATA NJUN/10,50,100/ A = -1.0 B = 1.0 WRITE(6,100) A,B 100 FORMAT(5X,'SEKIBUN FOR 1/SQR(PI)*EXP(-X*X) A', F5.2, 3X, 'B',F5.2) DO 10 K = 1, 3 N = NJUN(K) H = (B - A) / N S = ( FX(A) + FX(B) ) / 2.0 DO 20 J = 1, N - 1 X = A + H * J S = S + FX(X) 20 CONTINUE S = S * H WRITE(6, *) 'BUNKATSU = ', NJUN(K), 'MENSEKI =', S 10 CONTINUE STOP END C FUNCTION FX(X) PARAMETER (PI = 3.14159265) FX = EXP(- X * X) / SQRT(PI) RETURN END >>392 393 名前：392 [2006/06/05(月) 23:45:25 ] 小数点をカンマにするのは、ドイツ・フランス。中点にするのがイギリス。 あと日本の縦書きの新聞とかも中点にする。 しか〜しｗ　FORTRANやアメリカ・日本では小数点はピリオドだｗｗｗ FORMAT文はWRITE文の中に書き込んだほうがいいのだが、 学校教育では、まぁ文番号つきのFORMAT文でやっとけｗ 関数副プログラムを使ったが、今の場合FX(X)は文関数でも書ける。 質問があれば答えるから、まぁなんちゅーの、もう少し理解するようにしろｗｗｗｗ 394 名前：デフォルトの名無しさん [2006/06/05(月) 23:48:56 ] >>389 行列の表示の仕方がわからない・・・ 395 名前：デフォルトの名無しさん mailto:sage [2006/06/06(火) 00:14:17 ] >>394 real A(2,3) data A/1.,2., 3.,4., 5., 6./ do i = 1, 2 write(*,*) (A(i,j), j=1,3) end do end とか program hoge implicit none real :: A(2,3) = reshape((/1.,2., 3.,4., 5.,6./), (/2,3/)) integer :: i do i = 1, 2 print *, A(i,:) end do end program hoge とか 関係ないがFortran 2003では (/.../) のかわりに [...] と書けるので 可読性が少し良くなる（のかな？） 396 名前：デフォルトの名無しさん [2006/06/06(火) 00:42:54 ] >>395 解決しました！ありがとうございます。 397 名前：デフォルトの名無しさん [2006/06/06(火) 00:42:56 ] >>394 もっと具体的に質問すれば、華麗なるFORMATを教えて進ぜようゾナｗ 398 名前：364 [2006/06/06(火) 17:02:03 ] >>392 ありがとうございます。家のパソコンにはシグインないので学校行ったらチャレンジしてみます。 399 名前：392 [2006/06/07(水) 12:36:56 ] >>398 シグインってcygwinのことかな？ｗ 念のため言っておくけどｗ　変数のAとBが積分範囲の下限と上限なので、 −１，１というのを−２，２にすれば当初の課題が達成されるｗ 400 名前：デフォルトの名無しさん mailto:sage [2006/06/07(水) 13:10:18 ] シグインなんてDLすればいいジャンジャン 401 名前：デフォルトの名無しさん [2006/06/07(水) 15:16:53 ] Evca(Ω,x)を求めるプログラムを作りたいのですが教えてください。 Evca(Ω,x)=Ei(Ω,x)+E0(Ω,x)+E1(Ω,x)+E2(Ω,x) Ei(Ω,x)=-1.79186*Z**(5/3)/r E0(Ω,x)=Z(2.21/r**2-0.916/r+Ecorr) Ecorr=-(0.115-0.031logr) E1(Ω,x)=2πZe**2((1-X)Z**Cs*(Rm**Cs)**2(1+2/3*u**Cs)+xZ**K(Rm**k)**2(1+2/3u**K)) E2(Ω,x)=1/2��(4πe**2/ΩG**2)**-1(V(G))**2*1/(1-f(G))*(1/ε(G)-1) �狽ﾍG/=0 Gは逆格子ベクトルG(n1,n2,n3) V(q)=(1-x)V**Cs+xV**k(q) Vi(q)=-4πZe**2/Ωq**2((1+u)cos(qRm**i)-usin(qRm**i)/qRm) (i=Cs,K) u=qRmcos(qRm)/sin(qRm)-qRmcos(qRm) わからない値は後から自分で入れるので教えてください 状態方程式の計算をするらしいです 402 名前：デフォルトの名無しさん [2006/06/07(水) 15:24:39 ] E2(Ω,x)=1/2��(4πε**2/ΩG**2)**-1(V(G))**2(1/(1-f(G)))(1/ε(G)-1) �狽ﾍG/=0です Ｇは逆格子ベクトルです。G(n1,n2,n3) これを求めるプログラムはどうしたらいいか教えてください わからない値は自分でいれます。お願いします。 403 名前：デフォルトの名無しさん [2006/06/07(水) 17:42:42 ] >>401 質問が曖昧すぎるぞい！ｗ もう少しうまく質問しないと助けたくとも助けられないｗ 見た感じ結晶の電子相関込みの凝集エネルギーかなんかを計算しているようだが， どの辺のなにがわからんのか，求めたいものは何なのか，もうちょっと他人が分かる 質問をしろｗｗｗｗ これだけで答えられる奴はいねーよｗ 丸投げしたいなら持ってる情報を全部出せ！ｗ　 誰かが回答するかどうかは別としてｗ 404 名前：デフォルトの名無しさん [2006/06/07(水) 17:50:46 ] >>401-402 セシウムとカリウムのアルカリ金属合金の凝集エネルギーを求めて なにかしたいようであることは分かった． まずは欲張らずに単体のアルカリ金属の凝集エネルギーは求めろよｗ そしたら，ｘ＝１．０または０．０の時の極限にも対応するので あとからチェックにも使える． 一体どの辺が，どのくらいわかんないのよｗ　 405 名前：401 [2006/06/08(木) 11:05:05 ] すいません。あほなんで基本的には全然わかりません。 状態方程式のＰを求めるためにやってると思うんですけど。 圧力Ｐは結晶エネルギーＥvca(Ω,x)の体積Ωに関係する微係数で与えられるらしいのです。 a=5.225A A=10**-8 Ω=(a**3)/2 rs=(9π/4)**(1/3)/Kf Ef=h**2*Kf**2/2m c=2(1+0.0155π/Kf)**-1 f(q)=q**2/2(Q**2+cKf**2) ε(q)=1+(1-f(q))*(4πZe**2/Ωq**2)*(3/2Ef)*(1/2+((4Kf**2-q**2)/8Kfq)*log(2Kf+q/2Kf-q)) CsのＨ近似の場合Rmは4.610 uは-0.9605 KのＨ近似の場合Rmは3.019 uは-0.5682 せめて�濡/=0 のＧベクトルが０以外の時を全部足す方法が知りたいです。 お願いします。 406 名前：401 [2006/06/08(木) 13:57:06 ] k=5.225 a=k/(0.529177) o=a**3/2 z=1 rs=(3*o/z*4*3.14)**(1/3) Kf=(9*3.14/4)**(1/3)/rs Ef=Kf**2 c=2(1+(0.0155*3.14)/Kf)**-1 Ei=(1.79186*z**(5/3))/3*o*rs E0=z*(2.21/rs-0.916/rs+0.031/3*o) Rmk=3.019 uk=-0.5682 E1=-2*3.14*z*2/o*(z*(Rmk)**2*(1+2/3*uk)) k=5.225 a=k/(0.529177) o=a**3/2 z=1 rs=(3*o/z*4*3.14)**(1/3) Kf=(9*3.14/4)**(1/3)/rs Ef=Kf**2 c=2(1+(0.0155*3.14)/Kf)**-1 Ei=(1.79186*z**(5/3))/3*o*rs E0=z*(2.21/rs-0.916/rs+0.031/3*o) Rmk=3.019 uk=-0.5682 E1=-2*3.14*z*2/o*(z*(Rmk)**2*(1+2/3*uk)) 407 名前：401 [2006/06/08(木) 13:58:08 ] あとはE2(Ω,x)=1/2��(4πe**2/oG**2)**-1(V(G))**2*1/(1-f(G))*(1/ε(G)-1) �狽ﾍG/=0 を数値微分して全部足すと出るんですけどここがわからないです。 G=2π/a(n1*n1+n2*n2+n3*n3) 収束すると思うのである程度までの結果を出していきたいです。 n1=n2=n3/=0の時以外を計算させる方法がわかりません。数値微分もわかりません。 V(G)=-4*Z*2/oG**2((1+uk)cos(G*Rmk)-uk*sin(G*Rm)/G*Rm) ε(G)=1+(1-f(G))*(4*3.14*z*2/o*G**2)*(3/2*Ef)*(1/2+(4*Kf**2-G**2)/(8*Kf*G)*log((2*Kf+G)/2*Kf-G) f(G)=G**2/4*(G**2+c*Kf**2)+G**2/4*c*Kf**2 です。 408 名前：デフォルトの名無しさん [2006/06/08(木) 15:28:18 ] >>405 ＯＫ エネルギーを体積で微分すれば圧力になるのは熱力学の帰結だ． しかしこれはFORTRANというより物理の質問だなｗ ＞せめて�濡/=0 のＧベクトルが０以外の時を全部足す方法が知りたいです。 ＦＯＲＴＲＡＮスレ的には，これは３次元逆格子ベクトルに関する和なので， もっとも素朴には３重のＤＯ　ＬＯＯＰを回せばよいとい答えになる． DO n1 = -n, n 　DO n2 = -n, n 　　DO n2 = -n, n IF(n1==0 .and. n2==0 .and. n3==0) CYCLE　 云々 END DO END DO END DO 実際は結晶の対象性があるので，このループの１／８〜１／２４位回せばいいはず． 実際の計算には逆格子ベクトルを求める必要がある．これは結晶構造に依存する． もし今アルカリ金属が対象ならbcc構造ということになるだろう． その逆格子ベクトルはfccと同じものになる． これは固体物理の教科書に載っているから，一般式をつくるか， 手で求めて表を配列に入れることになる． DO　LOOPの中にn1,n2,n3の組み合わせが逆格子ベクトルにふさわしいかどうかを 判別する式を，G=0を飛ばすのと同じように，いれてやればいい． 409 名前：デフォルトの名無しさん [2006/06/08(木) 17:39:27 ] >>407 Gは逆格子ベクトルそのものではなくて、その長さだから（ルート忘れてるみたいだが） たんに逆格子点に相当するn1,n2,n3が分かっていればいい。 後は式の通りに書いて、逆格子点について足せばいい。 FORTRANの疑問ではないので、ここで聞く前に周囲の人に聞いた方が早いと思う。 数値微分は単に、格子間距離aを、特定のaとａ＋ｄａに取ったときに計算して 差分をｄａでわればいい。ｄａは小さい数だ。 410 名前：デフォルトの名無しさん [2006/06/09(金) 00:56:50 ] fortran90で、数式が二行に渡るとエラーがでるのですが、 どうしたらいいですか？ 411 名前：デフォルトの名無しさん [2006/06/09(金) 01:58:42 ] >>410 一行目のしりに＆をつけろ。 C言語やPASCALと違って、FORTRANは終端文字（；）が不要なかわり行センシティブ。 ゆえに継続行の記号を明示する必要がある。 漏れはFortran方式が好き。 412 名前：デフォルトの名無しさん mailto:sage [2006/06/09(金) 02:43:29 ] Ruby や Haskell のような方式もあるぞい 413 名前：401 [2006/06/09(金) 13:33:05 ] >>408 >>409 わざわざありがとうございます。物理の質問でしたねｗ ＣＹＣＬＥとはどういう意味ですか？ 結果にＮａＮとでるのはどういうことですか？ お願いします。 414 名前：デフォルトの名無しさん [2006/06/09(金) 14:08:18 ] >>413 CYCLEというのはFortran90で導入された命令で，DO文の中で使う． それ以降の文を実行すると無くDO行へ戻る． つまり，次のループ計算を実行する． これと類するものにEXIT命令があって，これはEND　DOの次の行に飛んで DO　LOOPを脱する． NaNとは数値に非ず．Not A　Numberの略号だ．配列がはみ出してとんでもない値を 書いているか，オーバーフロー等を起こしている． SUBCHKなどのデバッグオプションをつけてみそ． 415 名前：デフォルトの名無しさん mailto:sage [2006/06/09(金) 17:47:49 ] オーバーフローだと普通は ±inf じゃない？ NaN は実数での負値の平方根とか、未定義の演算結果かな。 416 名前：414 [2006/06/09(金) 21:27:34 ] ＞オーバーフローだと普通は ±inf じゃない？ おっしゃるとおり．漏れの間違いだ． 負の平方根とか関数の定義域外呼び出しでもでるね． 漏れが出会うNaNは，ゼロ割か，配列はみ出しでとんでもないメモリー領域の アクセスしたときが多いかな．個人的にはゼロ割で実行停止しなくなったのでとても助かる． NaNを厳密に追求したければ，IEEE745の浮動小数点の定義を調べればいいが気力那須ｗｗ 417 名前：初心者 [2006/06/09(金) 22:20:13 ] 最近、大学でFortranの授業が始まりました。 今日で授業が三回目だったのですが次のような課題がだされました。 おそらく基礎中の基礎だと思うのですが、授業はプリントに書かれた プログラムを説明もなくそのままうつして実行させるような授業ですので この課題を読んだ時日本語かどうかすら疑いました。 　いろいろと調べているのですがなかなかよくわかりません。よろしければ どなたかこの課題について解答していただけませんか？ 　あと、本当に自分勝手で申し訳ないのですが、少しでいいので解説もして いただけたらとても嬉しいです。よろしくお願いします。 （１）subroutineを呼ぶときの、変数の引き渡し方について説明せよ。 特に、配列変数を引き渡す方法についてプログラム名をあげて説明せよ。 （２）Fortranで文字列を取り扱うにはどうすればいいか調べて、プログ ラム例をあげて説明せよ。 （３）ＲＥＡＤ文やＷＲＩＴＥ文のＩ，Ｅ，Ｆ，Ａ形式について調べて 解説せよ。またプログラム例と印字例を示せ。 418 名前：初心者 [2006/06/09(金) 22:24:14 ] すいません。（１）のなかのプログラム名は間違いでプログラム例でした。 419 名前：デフォルトの名無しさん mailto:sage [2006/06/09(金) 23:56:51 ] 調べてって書いてるんだから、調べれば？ プリントか教科書あるんじゃないの？ 420 名前：初心者 [2006/06/10(土) 00:12:18 ] プリントは説明がなくて・・・教科書もないんです・・・やっぱりちょっと 調べればわかることみたいですね。 421 名前：デフォルトの名無しさん mailto:sage [2006/06/10(土) 00:20:38 ] FORTRAN で検索したら色々出てくるよ。 422 名前：初心者 [2006/06/10(土) 00:23:04 ] ありがとうございます！いろいろ調べて自分なりの解答作ってまた書き込み にくるんでよかったらみてください。 423 名前：デフォルトの名無しさん [2006/06/10(土) 00:29:05 ] 本物のプログラマって現存しますか? 424 名前：デフォルトの名無しさん mailto:sage [2006/06/10(土) 00:34:05 ] あれは最初からジョークだからw 色々示唆には富んでるけどね。 425 名前：デフォルトの名無しさん [2006/06/10(土) 00:44:15 ] >>419 教えてやれよｗ （１）なんて、参照渡しと値渡しの説明を求められていたら初心者には大変だぜ！ しかし、配列の先頭要素のアドレスしか渡していないことを知らないと７７では 実戦的なプログラムは書けない。 （２）はFORTRAN７７以降ならCHARACTERがあるが、FORTRAN６６だと整数型に ENCODE／DECODE文で文字列突っ込んだりしてるし。 （３）FORMAT書式は初心者がつまづく、難所中の難所だろう。 と、>>417をびびらせることを書いてみたりしてｗ　うふｗ 426 名前：デフォルトの名無しさん [2006/06/10(土) 03:03:51 ] data init/'12345'X,'23456'X,'34567'X,'45678'X/ ↑ '整数'X　が何を意味するのか分かりません。 よろしくお願いいたします。 427 名前：デフォルトの名無しさん [2006/06/10(土) 03:50:05 ] '整数'Xは16進数を表すものでした。 お手数をかけました。 428 名前：デフォルトの名無しさん mailto:sage [2006/06/10(土) 10:51:20 ] >>417 大阪大学工学部応用理工学科二回生乙。 429 名前：デフォルトの名無しさん mailto:sage [2006/06/10(土) 17:02:15 ] これが噂に聞いた「ソシャる」というやつか 430 名前：デフォルトの名無しさん [2006/06/10(土) 17:30:31 ] >>422 臆せずに書き込むように。 みんなの勉強になる。 431 名前：デフォルトの名無しさん [2006/06/11(日) 13:31:23 ] 次の宿題もってこいや！！ 432 名前：364 [2006/06/11(日) 16:21:12 ] >>392を使わせてもらったんですけど、実行できませんでした。 cygwinを起動してコンパイルするとこまでは行けたんですけど(g77 report.f)、 ./a.exewoをしてもcommand not found になります。どこかエラーがあるのでしょうか？ 433 名前：364 [2006/06/11(日) 16:22:43 ] ./a.exewo→./a.exeの間違いです。 434 名前：デフォルトの名無しさん [2006/06/11(日) 16:40:56 ] >>432 command not found ということは、a,exeが見つからないということなので、コンパイルに失敗している。 それ以前に何かエラーメッセージが出ていないか見て味噌。 まずはｌｓしてa.exeができているか見て味噌。 435 名前：ケロ [2006/06/11(日) 23:33:31 ] どなたか、教えてください;:il:il|;l|;il:i(-ω-｀；)ll|l|il|;:il|! マジでなやんでます…。 1/√(π)exp(-x^2) 中央値０　偏差１　の正規分布表をあたえる。 [−１．１]および[−２．２]の範囲に存在する確立はいくらか？ で台形法で解けと…。 プログラムがわかんないです…。 たすけてください… 436 名前：デフォルトの名無しさん mailto:sage [2006/06/11(日) 23:40:51 ] >>435 >>364 437 名前：ケロ [2006/06/11(日) 23:43:29 ] 一緒だ…。 438 名前：ケロ [2006/06/12(月) 00:02:52 ] なんとなく、わかりました。ありがとーございます<(_ _)> 439 名前：364 [2006/06/12(月) 00:09:02 ] >>434 ありがと。もっかいやってみます。Invalid first character at (^)とか出るんですよね・・・。 >>438 同じ授業を受けてる気がする・・・。特定はやめてね・・・ 440 名前：ケロ [2006/06/12(月) 00:12:48 ] ＞＞439 それ、思った…。 441 名前：364 [2006/06/12(月) 00:28:30 ] >>440 ここでもその名前でお願い ex11.2ch.net/test/read.cgi/campus/1149336703/ 442 名前：デフォルトの名無しさん mailto:sage [2006/06/12(月) 00:35:20 ] >>439 もしかしてプログラムを「１桁目から」書いてないか？ 2chの表示の都合上、ここではプログラムが行の先頭から書かれているように見えるが 実際のFORTRAN77のソースは行の各桁に何を書くかも文法上で決まっている。 そこら辺は初歩の初歩でどんな入門書にも載っているようなことなので 少しは教科書なりWeb上の入門記事なりを読んで勉強してくれ。 443 名前：364 [2006/06/12(月) 00:46:43 ] >>442 Σ(ﾟДﾟ)！すいません、そのままコピペしてました。 444 名前：デフォルトの名無しさん [2006/06/12(月) 02:21:28 ] >>443 君は２ｃｈ専用ブラウザを使ってるかい？ 使っているなら>>392の引用符を見れば一応ある程度インデントが効いて 本来の姿に近くなる。 >>392の中で自分を引用しているのは、そうやって見るためだ。 もうちょっと最小限のところは勉強しろｗｗｗｗｗｗ 他のスレならぶちキレられてるぞｗｗｗｗｗｗ このスレはVIPレベルなのでOKだがな（＾ω＾；） >>435のケロも遠慮せずに書き込めよｗ 445 名前：401 [2006/06/12(月) 13:09:42 ] >>414 わかりました。ありがとうございます。 一応作ってみたのですがエラーがでてコンパイルできません。 初めて作るのでめちゃくちゃかもしれませんがおかしいところを教えてください。 お願いします。 k=5.225 a=k/(0.529177) o=a**3/2 z=1 rs=(3*o/z*4*3.14)**(1/3) Kf=(9*3.14/4)**(1/3)/rs Ef=Kf**2 c=2*(1+(0.0155*3.14)/Kf)**-1 Ei=(1.79186*z**(5/3))/3*o*rs E0=z*(2.21/rs-0.916/rs+0.031/3*o) 　Rmk=3.019 uk=-0.5682 E1=-2*3.14*z*2/o*(z*(Rmk)**2*(1+2/3*uk)) E2 = 0 E4 = 0 da = 0.001 446 名前：401 [2006/06/12(月) 13:12:03 ] WRITE(*,*)'(n4)' READ(*,*) n4 DO 10 n1= -n4 , n4 DO 10 n2= -n4 , n4 DO 10 n3= -n4 , n4 　IF(n1==0.and.n2==0.and.n3==0) CYCLE G=2*3.14/a*(n1*n1+n2*n2+n3*n3) f=G**2/4*(G**2+c*Kf**2)+G**2/4*c*Kf**2 e=1+(1-f)*(4*3.14*z*2/o*G**2)*(3/2*Ef)*(1/2+(4*Kf**2-G**2)/ &(8*Kf*G))*log(2*Kf+G/2*Kf-G) V=(-4*Z*2/oG**2)*((1+uk)*cos(G*Rmk)-uk*sin(G*Rmk)/G*Rmk) E3=(4*3.14*2/o*G**2)**-1*V**2*(1/(1-f))*(1/e-1) E2 = E2 + E3 V1=(-4*Z*2/(o+da)G**2)*((1+uk)*cos(G*Rmk)-uk*sin(G*Rmk)/G*Rmk) e1=1+(1-f)*(4*3.14*z*2/(o+da)*G**2)*(3/2*Ef)*(1/2+(4*Kf**2-G**2)/ &(8*Kf*G))*log(2*Kf+G/2*Kf-G) E5 =(4*3.14*2/(o+da)*G**2)**-1*V1**2*(1/(1-f))*(1/e1-1) E4 = E4 + E5　 10 CONTINUE E6 = E2/2 E7 = E4/2 E8 = (E7-E6)/da P=Ei+E0+E1+E8 447 名前：401 [2006/06/12(月) 13:14:02 ] WRITE(6,601) P 601 FORMAT(1H , F15.8) STOP END エラー fortcom: Error: Illegal character in statement label field [�| fortcom: Error: Illegal character in statement label field [・ fortcom: Error: keita.f, line 55: Syntax error, found END-OF-FILE when expecting one of: <LABEL> <END-OF-STATEMENT> ; BLOCK BLOCKDATA PROGRAM TYPE COMPLEX BYTE CHARACTER ... . ^ compilation aborted for keita.f (code 1) あほすぎて迷惑かけます。 448 名前：素人 [2006/06/12(月) 14:38:08 ] -1以上1以下の値ｘを入力してe^xの近似式を出力するプログラム。e^xの近似式はe^x=1+x+x^2/2!+x^3/3!+..+x^n/n!を利用。 (最後の項に絶対値が0.00001より大きい間、項を増やして足し続ける)また、適正なxの値が入力される限り、何度もこの計算ができるようにすること。って課題が出されたんですけどまったく分かりません。誰か助けてください！ 449 名前：デフォルトの名無しさん [2006/06/12(月) 15:27:53 ] >>445 我慢して聞いてくれｗ 基礎が全くなっていないので，ちょっと見ただけでもまともに動くとは思われない． 周りに聞く人がいるなら聞いた方がいい．掲示板で説明出来る範囲を超えているｗｗｗｗ １．まず変数の型がむちゃくちゃ．第一行目から，整数型のＫに実数を入れている． 式の途中に（１／３）みたいなのがいっぱいあるが，整数の場合はこれは０になる． ようするに型に関する認識が全くなってないｗ　整数と実数は違うものだと思った方がいい． 実数に整数を代入する方は，一応型の変換が自動でなされるが，明示的に実数にしておいた ほうがいい． ２．有効桁に関する認識もなってない．円周率を３．１４で入れているが，これでは有効桁が 小数点以下２桁になってしまう．そもそも常数は，ＰＩとかちゃんとした名前の付いた変数 ないしパラメータ変数にすべきもの． ３．式も所々演算子が抜けているし，演算子の優先度に関して間違っているとしか 思えないところがある（割り算の分母）． ４．微分を取りにゆくときのセンスも悪い．まず特定の格子間距離に対してエネルギーを 求めるサブルーチンをつくって，それを２回異なる格子間距離について呼ぶべき． とにかく，なんというか，プログラミングの基本的センスがなっていないし， 問題自体に対する理解もなってないようなので， 今のままでは他者からのかなりの助けがないと完成はおぼつかない． 周りに聞く人がいれば，頭を下げて一から教えてもらえ． とにかく道はティッペラリーより遙かかなただｗ 450 名前：デフォルトの名無しさん [2006/06/12(月) 15:46:19 ] >>448 ほれ、７７形式で昔のテキスト風に芋くさく書いておいたぞい。 PROGRAM VIPPER 1 WRITE(*, *) 'INPUT x (-1.0<= x <= 1.0)' READ (*, *) x IF (ABS(x) > 1.0) STOP e = 1.0 f = x i = 1 10 CONTINUE e = e + f IF (ABS(f) < 0.00001) GOTO 999 i = i + 1 f = f * x / REAL(i) GOTO 10 999 CONTINUE WRITE(*, '( 2(a, f8.5) )') 'EXP(', x, ')', e WRITE(*, *) EXP(x) GOTO 1 END 451 名前：450 mailto:sage [2006/06/12(月) 15:48:58 ] すまんｗ　ついうっかりFortran90の癖を出してしまった。 ７７としては以下の部分を置き換えておいてくれｗ IF (ABS(x) .GT. 1.0) STOP IF (ABS(f) .LT. 0.00001) GOTO 999 452 名前：デフォルトの名無しさん mailto:sage [2006/06/12(月) 16:38:20 ] 大学はいってから、別にプログラム関係の学部でもないのに fortranをやらされ、まったくわからない… なんか時分秒の引き算のプログラムを作らなければいけない。 　開始時刻と終了時刻を指定すると、経過時間 を計算し端末に表示するプログラムをつくりなさい。 時刻はhh:mm:ssの形式で指定されるものとする。 終了時刻が開始時刻より早い場合には、その旨を画面に 表示して、再入力させること。 だそうだ。 PRPGRAM TIME REAL A,B,C,D,E,F,G,H,I READ(*,*)A,B,C,D,E,F X=A*60*60+B*60+C Y=D*60*60+E*60+F IF(X.GT.Y)THEN Z=Y-X G=Z/60/60 M=Z-G*60*60 H=M/60 I=M-H*60 WRITE(*,*)G,H,I ELSE WRITE(*,*)'ERROR' END IF STOP END 453 名前：452 mailto:sage [2006/06/12(月) 16:39:06 ] こんなのが限界で、hh:mm:ssと入力できないし、出力もできない。 ２４時間以降もERRORにしなければならないが、どうやるかわからない。 ;←を打つには書式を変えろって言われたけど習ってないし、 教えてもくれなかった＿|￣|○ 助けてください。こんな簡単なことを聞いてすいません… 454 名前：デフォルトの名無しさん mailto:sage [2006/06/12(月) 17:19:38 ] >>453 ほれ、出来るだけ藻舞のオリジナルを尊重しつつ77で作った。 変数は整数型にしないといかんべ。 MODは余りを求める関数だ。マニュアルみそ。 入力をhh:mm:ssの形式にしたければ、入力エラーは無いことを仮定すれば、 ＷＲＩＴＥ文のＦＯＲＭＡＴをちょっと直してＲＥＡＤ文の方にも書いてやればいい。 PROGRAM TIME INTEGER IH0, IH2, IM0, IM1, IS0, IS1, IT0, IT1 INTEGER IH, IM, IS, IT 1 READ(*, *) IH0, IM0, IS0, IH1, IM1, IS1 IT0 = IH0 * 60 * 60 + IM0 * 60 + IS0 IT1 = IH1 * 60 * 60 + IM1 * 60 + IS1 IF(IT1 .GT. IT0)THEN IT = IT1 - IT0 IH = IT / (60 * 60) IT = MOD(IT, 60 * 60) IM = IT / 60 IS = MOD(IT, 60) WRITE(*, '( I2.2, A1, I2.2, A1, I2.2 )' ) IH, ':', IM, ':', IS ELSE WRITE(*,*)'ERROR' GOTO 1 END IF STOP END >>454 455 名前：454 [2006/06/12(月) 17:41:03 ] ２４時間がどうのこうの謂う所は自分でやってちょｗ ＩＦ文で飛ばすだけだし、>>452を見る限り自力で出来る感じがするｗ 456 名前：410 [2006/06/12(月) 18:10:46 ] >>411 遅レスですいません！たすかります、どうもっす！ 457 名前：デフォルトの名無しさん [2006/06/12(月) 18:11:27 ] >>452 ＞大学はいってから、別にプログラム関係の学部でもないのに ＞fortranをやらされ、まったくわからない… 来るべき２１世紀の脱工業化社会には、電子計算機の知識が欠かせないのだ。 計算機言語のサンスクリット語たるＦＯＲＴＲＡＮを学べることを誇りに思って励め！

---

---

[ 続きを読む ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef