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

1 名前：デフォルトの名無しさん [2006/01/24(火) 09:48:23 ] このスレッドは、他のスレッドでは書き込めない超低レベル、  もしくは質問者自身何が何だが分からない質問を勇気を持って書き込むスレッドです。  FORTRAN使いが優しくコメントを返しますが、  お礼はFORTRANの布教と初心者の救済をお願いします。  581 名前：デフォルトの名無しさん [2006/07/01(土) 19:14:04 ] >>579 関数の中に関数は定義できる。 引数の形が任意の関数を渡すことはできない。 数学でいう汎関数のようなものはできない。 できないのは、ほとんど必要とされていないから。 普通は、以下のようにサブルーチンは固定のままで話が済む。 program test S = f(3.0, 4.0, 4.0, 4.0, 1.0) END function f(x, a, b, c, d) f = 3.0 * x + g(x, c, d) + a + b return end function g(x, c, d) g = x**2 + h(x, d) + c end function h(x, d) h = 4.0 * x + d end 582 名前：デフォルトの名無しさん mailto:sage [2006/07/03(月) 13:36:18 ] どもできましたｗ function　fの引数にhを入れるのを忘れてましたｗ でもひとつうんこなしつもんなんだけどぉ functionの中に配列渡すにはどうしたら良いの？ real Un(100) integer Ko Ko=36 X=UNK(Un,Ko) end function UNK(Un,Ko) real Un integer Ko UNK=Un(Ko) end 583 名前：デフォルトの名無しさん mailto:sage [2006/07/03(月) 13:54:35 ] >>582 ネタでなかったらそんな聞き方はしない筈 584 名前：デフォルトの名無しさん [2006/07/03(月) 14:19:38 ] >>582 少しは前レス嫁ｗｗ　いいけどなｗｗ 受け取る側で配列であることを示せばいい。 function UNK(Un,Ko) real Un(*) UNK=Un(Ko) end たとえば、こういうの。配列の最後の次元については、＊で示してもいい。 ただ、普通はもう少し上品に配列の定義長を渡す。 多次元の場合は最後の次元以外でどうしても必要になる。 real Un(100) Ko=36 X=UNK(100,Un,Ko) end function UNK(Mn, Un, Ko) real Un(Mn) UNK=Un(Ko) end FORTRAN77では配列の先頭要素のアドレスしか渡しておらず、配列のサイズに関する情報は 別個に明示的に渡さないといけない。いろいろ便利だが、バグの原因にもなるので、やれるなら F90で書いてきっちりやったほうがいい。 けど、好きにやればいいさｗ 585 名前：デフォルトの名無しさん mailto:sage [2006/07/04(火) 14:06:28 ] ども access violationとかでる（＞＜） 変数宣言の順番変えたらうまくいったりする場合もあったり・・ これってバグってやつ？ f90でやったほうがいいの？それともCってやつのほうがいいの？ fortran以外だとdelphi pascalしかしらない・・・ 586 名前：デフォルトの名無しさん mailto:sage [2006/07/04(火) 14:19:48 ] つまり配列を渡す場合にはいったん外部で配列の値を入力した変数を渡してやれば良いってことジャン？？ めんどくせーーｗ こんなところで10時間近く悩んでたなんてあほだなｗｗｗｗｗｗｗ fortranってしょぼくね？ｗ real Un(100) integer Ko Ko=36 Un=Un(ko) X=UNK(Un) end function UNK(Un) real Un UNK=Un 587 名前：デフォルトの名無しさん mailto:sage [2006/07/04(火) 14:22:14 ] あ　これじゃだめかU=Un(Ko)とかにしないとだめなんだな 588 名前：デフォルトの名無しさん mailto:sage [2006/07/04(火) 14:45:26 ] fortranまじしょぼいｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗ 引数の順序変えたらまともに動かなくなったり動くようになったりｗｗｗ なんでこんなに不安定なの？ 使い物にならないジャン？ 589 名前：デフォルトの名無しさん [2006/07/04(火) 15:05:57 ] >>588 しょぼいのはFORTRANではない、お前だｗｗｗ 単に使い方を誤解しているだけだ。 引数の順番を変えてうまく行ったり行かなかったりするというのは、 じつはどっちにしろうまく行ってないということだｗ PASCALができるなら、F90の方が向いている。Module言語的になっているから。 しかし、書き込みを見ている限りPASCALで出来るとも思われぬがｗうひょｗｗｗ 590 名前：デフォルトの名無しさん mailto:sage [2006/07/04(火) 18:06:51 ] ごめんやっぱしょぼいのおれだったｗｗｗｗ ってか72文字しかかけないんだって？ だったら７２文字以上書いた時点でエラー出せよｗ コンパイルスンナｗ しかも変数の宣言間違ってたらエラー出せよｗ 591 名前：デフォルトの名無しさん mailto:sage [2006/07/04(火) 18:21:33 ] まあ、F77 がしょぼいのは確かだけどな。 592 名前：デフォルトの名無しさん mailto:sage [2006/07/04(火) 18:45:36 ] >>590 FORTRAN７７には予約語とかはないし、空白とかをトークンとして区切るわけでもない。 コンパイラはまず空白を削除して頭から解釈してゆく。 したがって７２文字より先の部分が切り捨てられても、解釈が通ればエラーにはならない。 大体、コンパイラが厳しければ厳しいで文句を言うんだから、ゆるいのを喜べ！ FORTRAN６６なんかDO　LOOPから飛び出して、なんかやった後またそのLOOPの中に 戻っていけんるんだぞ！こんなことできる言語はなかなかないぜ！ 593 名前：デフォルトの名無しさん mailto:sage [2006/07/04(火) 19:02:38 ] それは利点じゃねーw 594 名前：デフォルトの名無しさん mailto:sage [2006/07/04(火) 20:10:31 ] >>590 下品でない変数名にしようよ。 貴婦人が見ているかもしれないんだから。 595 名前：デフォルトの名無しさん [2006/07/05(水) 03:24:52 ] >>594 KISS MY ASS！　F*CK　U！ FORTRANの変数名が６文字の制限下では、Four Letter Words がぴったりだぜ！ 596 名前：デフォルトの名無しさん [2006/07/05(水) 05:08:37 ] さぁテポドンの弾道をFORTRANで計算しろ！ 597 名前：デフォルトの名無しさん mailto:sage [2006/07/05(水) 07:31:54 ] >>596 do i=1.10 宇宙の果てまで、ちゅぼーん 598 名前：デフォルトの名無しさん mailto:sage [2006/07/05(水) 14:07:32 ] >>597 ほんとに１０発撃ったらしいぞｗ　４発は不発だったらしいがｗｗｗ 599 名前：デフォルトの名無しさん [2006/07/05(水) 19:01:15 ] 問題：L行M列とM行N列の行列の積を計算するサブルーチン副プログラムを作れ。 サブルーチン副プログラムを呼び出す主プログラムも作成し、結果を出力すること。 主プログラムでは二つの行列のデータ読み込みと結果の出力を行うこととする。 (この後プログラムが正しく動作することを確認するための例題が続きます 途中の文章をいれたプリント文はそれを入力させるためのものです) real function SEKI(A,B,C,J,K,L) real A(2,2),B(2,3),C(2,3) c do 30 I=1,L do 20 J=1,N C(I,J)=0 do 10 K=1,M C(I,J)=C(I,J)+A(I,K)*B(K,J) 10 continue 20 continue 30 continue c return end 600 名前：デフォルトの名無しさん [2006/07/05(水) 19:02:58 ] c do 70 K=1,N print*,C(1,K) 70 continue do 80 K=1,N print*,C(2,K) 80 continue end c print*,'1,-2,3,0,4,-1,2,1,0,-3を入力して下さい' do 50 I=1,L do 40 K=1,M read*,A(I,K) 40 continue 50 continue 601 名前：デフォルトの名無しさん [2006/07/05(水) 19:03:53 ] do 70 K=1,M do 60 J=1,N read*,B(K,J) 60 continue 70 continue c end エラーが出るんですがもうどこが間違ってるのか全然わかりません。お手上げです。 誰か分かる方、ご教授願います。 602 名前：デフォルトの名無しさん [2006/07/05(水) 19:26:34 ] >>599 function ではなく　SUBROUTINE だべｗ 配列引数のサイズを、固定長で与えては一般性のあるプログラムにならない。 （これを次元の宣言と割り切る、はみ出し無視のFORTRAN６６時代の猛者なら別だがｗｗｗ） あとは、まぁまぁ大丈夫でないの？　 ７０番が2回出てくる上に、妙なところにend文が割り込んで理解に苦しむが、 こぴぺ時の間違いとみなしておくｗ 603 名前：デフォルトの名無しさん mailto:sage [2006/07/05(水) 20:05:21 ] できたおｗ program hoge implicit none real, allocatable :: A(:,:), B(:,:), C(:,:) integer :: i,j, l,m,n read *, l, m, n allocate(A(l,m), B(m,n), C(l,n)) read *, ((A(i,j), j=1,m), i=1,l) read *, ((B(i,j), j=1,n), i=1,m) call seki(A, B, C, l, m, n) do i=1,l print *, (C(i,j), j=1,n) end do deallocate(A, B, C) stop contains subroutine seki(A,B,C, l,m,n) integer, intent(in) :: l,m,n real, intent(in) :: A(1:l,1:m), B(1:m,1:n) real, intent(out) :: C(1:l,1:n) C = MATMUL(A, B) end subroutine seki end program hoge 604 名前：デフォルトの名無しさん [2006/07/05(水) 20:19:03 ] >>603 おまｗ　 ９０でやったらわざわざMATMUL相当のサブルーチン作る意味無いｗ ついでに７７でも作ってやれｗ 605 名前：デフォルトの名無しさん mailto:sage [2006/07/05(水) 20:53:04 ] 77だとこんな感じだお SUBROUTINE SEKI(A,B,C,L,M,N) REAL A(L,M), B(M,N), C(L,N) DO 30 I=1,L DO 20 J=1,N C(I,J) = 0. DO 10 K=1,M C(I,J) = C(I,J) + A(I,K)*B(K,J) 10 CONTINUE 20 CONTINUE 30 CONTINUE RETURN END C PROGRAM HOGE REAL A(3,3), B(3,2), C(3,2) L=3 M=3 N=2 READ(*,*) ((A(I,J), J=1,M), I=1,L) READ(*,*) ((B(I,J), J=1,N), I=1,M) CALL SEKI(A,B,C, L,M,N) DO 10 I=1,L WRITE(*,*) (C(I,J), J=1,N) 10 CONTINUE STOP END 606 名前：学生（アホ１６歳） mailto:hage [2006/07/05(水) 23:57:31 ] 以下の課題についてですが… 【自然対数 e は，1 + (1)-1 + (2*1)-1 + (3*2*1)-1 + ・・・ + (10*9*8*・・・・*2*1)-1 というように近似することができます．この値をディスプレイに表示するプログラムを作成しなさい．】 ちなみに７７で書かないと点くれないそうで…（苦笑 ７７で階乗計算って…どうすればいいのですか？ 607 名前：デフォルトの名無しさん mailto:sage [2006/07/06(木) 00:30:23 ] >>606 1! = 0!*1 = 1*1 = 1 2! = 1!*2 = 1*2 = 2 3! = 2!*3 = 2*3 = 6 4! = 3!*4 = 6*4 = 24 　　： n! = (n-1)! * n k = 1 do 10 n=1, 10 k = n * k write(*,*) 'n = ',n, ', n! =', k 10 continue end 608 名前：デフォルトの名無しさん [2006/07/06(木) 01:18:11 ] >>606 階乗は>>607でいいが、１項毎に定義どうりに計算するのではなく 一つ前の項に違いだけ掛ければいい。 y=1, y = y*x/1, y = y*x/2, y = y * x/3 .............. が、さらに順を変えて尻から進めるほうがもっといい。 1 + x (1 + x(1 + x/3) / 2 ) / 1 609 名前：学生（アホ１６歳） mailto:hage [2006/07/06(木) 03:22:58 ] >>607-608 ありがとうございます！ 出来ました。 610 名前：デフォルトの名無しさん [2006/07/07(金) 11:11:35 ] fortranでC言語の#defineみたいに、プリプロセッサで定数指定できるのってありますか？ 611 名前：デフォルトの名無しさん mailto:sage [2006/07/07(金) 13:02:41 ] 拡張子を .F にしたらプリプロセッサを解析してくれるコンパイラも多い。 612 名前：610 [2006/07/07(金) 13:35:41 ] >>611 返信どうもっす。じゃあ、fortranの標準では#defineみたいなのないんですね？ 613 名前：デフォルトの名無しさん [2006/07/07(金) 13:56:35 ] >>611 標準では存在しない。 #define は、便利な反面、批判も多いのでたぶん意図的に新しい規格でも不採用とおもわれ。 ただ問題点を無くした条件コンパイルなどのプリプロセッサへの対応は論じつづけられている。 614 名前：610 [2006/07/07(金) 18:17:48 ] >>613 なるほど、納得いたしました。どうもありがとうございます！ 615 名前：デフォルトの名無しさん mailto:sage [2006/07/07(金) 20:37:06 ] 定数は parameter だな。 616 名前：デフォルトの名無しさん [2006/07/11(火) 19:31:21 ] 学校の課題をやっていただきたいのですが、、、 次の数列の和を求めよ。（DO文） 1/1+x=1-x+x**2-x**3+x**4-x**5・・・ 打ち切り条件　1.0*10**-6 提出明日です。よろしくおねがいします。 617 名前：デフォルトの名無しさん mailto:sage [2006/07/11(火) 22:20:20 ] 割り込みすいません… 分数の足し算プログラムですが、 どこが違うかわからないので教えてくださいませんか？ PROGRAM FRACTION INTEGER I,J,K,L,M,N READ(*,*) I,J,K,L CALL TUUBUN(I,J,K,L) WRITE(*,*)M,N CALL REDUCTION(M,N) M=M/A N=N/A WRITE(*,*)M,'/',N END SUBROUTINE TUUBUN(A,B,C,D) INTEGER A,B,C,D M=A*D+B*C N=C*D RETURN END 618 名前：デフォルトの名無しさん [2006/07/11(火) 22:24:36 ] >>616 FORTRAN７７で適当に作った。適当に直してくれ。 数学的には　-1<x<1 で収束する式だが、１に近づくにつれて収束が遅くなる。 つまり、たくさんの項の和を取らねばならない。しかし、延々とやると数値誤差が たまるので意味がなくなる。ゆえに適当なところであきらめて打ち切る。 PROGRAM VIP PARAMETER (eps = 1.0e-6) WRITE(6, *) 'INPUT x (-1 < x < 1)' READ(5, *) x s = 1.0 t = 1.0 DO 10 i = 1, 10000 t = - t * x s = s + t IF (ABS(t) .LT. eps) GOTO 999 10 CONTINUE STOP 'NOT CONVERGED' 999 WRITE(6, *) ' sum =', s, ' iteration', i WRITE(6, *) ' 1/(1+x) =', 1.0 / (1.0 + x) STOP END 619 名前：デフォルトの名無しさん [2006/07/11(火) 22:32:24 ] >>617 サブルーチンの変数の範囲について根本的に誤解している。 メインプログラムとサブルーチンでは変数はお互いに見えない。 ゆえに、通分した分子と分母はM、Nはサブルーチンの引数で返す必要がある。 CALL TUUBUN(I,J,K,L, m, n) SUBROUTINE TUUBUN(iA,iB,iC,iD, m, n) M=iA*iD+iB*iC N=iC*iD RETURN END あと、暗黙の型は死んでも守るべし。周りがなんと言おうと守るべし。 かっぱに尻子玉を抜かれそうになっても寧ろ暗黙の型守るべし。 620 名前：617 mailto:sage [2006/07/12(水) 10:25:30 ] >>619 途中からなくなっていたので、 もう一度直して PROGRAM FRACTION INTEGER I,J,K,L,M,N READ(*,*) I,J,K,L CALL TUUBUN(I,J,K,L) CALL REDUCTION(M,N) M=M/A N=N/A WRITE(*,*)M,'/',N END SUBROUTINE TUUBUN(iA,iB,iC,iD) INTEGER iA,iB,iC,iD M=iA*iD+iB*iC N=iC*iD RETURN END 621 名前：617 mailto:sage [2006/07/12(水) 10:26:19 ] >>620 SUBROUTINE REDUCTION(P,Q) INTEGER P,Q,X,Y X=P Y=Q 10 IF(X.EQ.Y) THEN GO TO 20 ELSE IF(X.GT.Y) THEN X=X-Y ELSE Y=Y-X END IF END IF GO TO 10 20 CONTINUE A=X RETURN END 622 名前：デフォルトの名無しさん [2006/07/12(水) 14:25:40 ] >>621 PROGRAM FRACTION READ(*,*) I, J, K, L CALL TUUBUN(I, J, K, L, M, N) CALL REDUCTION(M, N, IA) M = M / IA N = N / IA WRITE(*,*) M, '/', N END SUBROUTINE TUUBUN(N1, N2, ID1, ID2, NUMERA, IDENOM) NUMERA = N1 * ID2 + N2 * ID1 IDENOM = ID1 * ID2 RETURN END SUBROUTINE REDUCTION(NUMERA, IDENOM, IGCD) IX = NUMERA IY = IDENOM 10 IF(IX .EQ. IY) THEN GO TO 20 ELSE IF (IX .GT. IY) THEN IX = IX - IY ELSE IY = IY - IX END IF GO TO 10 20 CONTINUE IGCD = IX RETURN END >>622 623 名前：さおり [2006/07/12(水) 16:17:47 ] すいません！無理を承知でカキコします。 私は今理系で某国立大学に通ってるんですけど、水商売のバイトにハマってしまって全然授業についていってない状態なんです。 それでプログラミングの授業があるんですけど課題が出て「放物の運動（空気抵抗を考慮する）」のプログラムをfortranで作れって課題が出たんです。。 初速度と投げる角度を入力して一定時間ごと（１秒ごととか）の物体の位置をｘとｙで出すらしいんです。 私、授業にほとんど出てなくて本当にまったくわかってない状態なんです。 少し真面目な大学なんで私みたいな子は浮いてしまって、学校には助けてくれる友達もいないんです。 誰かプログラム作ってください。お願いします。この授業を落としてしまうと本当に留年なんです。おねがいします！！ 624 名前：デフォルトの名無しさん [2006/07/12(水) 16:43:42 ] >>623 まずはおっぱいうｐ汁 625 名前：デフォルトの名無しさん mailto:sage [2006/07/12(水) 21:23:23 ] >>623 担当教官をお店で接待すれば優もらえるおｗ 626 名前：デフォルトの名無しさん mailto:sage [2006/07/13(木) 00:14:44 ] 金があるなら留年しろ。 627 名前：デフォルトの名無しさん [2006/07/13(木) 01:04:04 ] 俺は少年なら助けてやるのだが・・・・女は興味ないｗｗｗ 628 名前：デフォルトの名無しさん mailto:sage [2006/07/13(木) 07:18:28 ] >>623 課題だけ書く方が答えが得られるすれのようです。 629 名前：さおり [2006/07/13(木) 10:47:04 ] 本当にお願いします！ いろいろ事情があって留年できないんです。 630 名前：デフォルトの名無しさん mailto:sage [2006/07/13(木) 12:22:06 ] >>629 課題は省略しないで出題されたまま提示するべき ヒントになるかどうか知らん ttp://www12.plala.or.jp/ksp/mechanics/resistdown/ 631 名前：デフォルトの名無しさん mailto:sage [2006/07/13(木) 12:26:25 ] >>623 空気定数って物体のサイズによって二通りの計算方法があるけど、どっちを使うの？ 632 名前：631 mailto:sage [2006/07/13(木) 12:29:35 ] >>629 >>630 のヒントの方の計算方法だと楽でいいよね。さあどっちだろうか？ 633 名前：631 mailto:sage [2006/07/13(木) 16:25:24 ] >>631 ×空気定数→○空気抵抗 634 名前：さおり [2006/07/13(木) 21:43:35 ] えーっと空気抵抗が速度の２乗に比例する場合らしいです。本当にありがとうございます。 635 名前：デフォルトの名無しさん mailto:sage [2006/07/13(木) 22:08:31 ] 解決したみたいだから、次の宿題。 636 名前：デフォルトの名無しさん [2006/07/13(木) 22:56:22 ] これから宿題をストーリーつきで頼む人は、 『ボク宿題が解けなくて困ってるんだ。 お兄ちゃん様に一緒に解いてほしいなぁ〜。云々』 という、もっとやる気を出させるストーリーを考えること！ 637 名前：デフォルトの名無しさん mailto:sage [2006/07/13(木) 23:00:37 ] >>634 ああ、めんどくさい方ね。 それ、解析的に解が求まらないから、微分方程式出して数値解析するしかないかな？ 一応、微分方程式を立ててみるけど、チェックしてみて vx, vy を速度の x,y 成分として、物体の質量を m 、重力定数を g、空気抵抗の 係数を k とすると m dvx/dt = -k vx √(vx^2+vy^2) m dvy/dt = -k vy √(vx^2+vy^2) - mg となるかな？あとは、授業で習った数値解析で vx, vy をまず求めるのだが…… 638 名前：デフォルトの名無しさん mailto:sage [2006/07/13(木) 23:19:43 ] >>636 兄貴ぃ〜！ 俺、宿題分かんないッス！！ 徹夜で手取り足取り教えて欲しいッス！！ 639 名前：デフォルトの名無しさん [2006/07/14(金) 00:15:18 ] >>638 サブ！　手取り腰取り教えてやるぜ！！ 640 名前：デフォルトの名無しさん [2006/07/14(金) 01:05:54 ] >>637 最後まで解けよｗ 641 名前：デフォルトの名無しさん mailto:sage [2006/07/14(金) 02:10:46 ] >>637 言いだしっぺのh(ry 642 名前：641 mailto:sage [2006/07/14(金) 02:11:58 ] 間違った 637でなくて>>640 643 名前：637 mailto:sage [2006/07/14(金) 10:13:37 ] >>640 そうじゃなくて、学校で習った数値解析法で解かないと 合格になる気がしないんですけど。XX法とか言ってもらえれば 続きも出来るのですが…… 644 名前：デフォルトの名無しさん mailto:sage [2006/07/14(金) 11:08:41 ] >>634の書き込みを見ると、「らしい」なんて書いているから、課題の正確な内容を 把握してなさそう。仮に折角>>637が解いたところで「やっぱり違った」とかで別解を 要求してきそうだ。 645 名前：デフォルトの名無しさん [2006/07/14(金) 13:55:36 ] C++スレ見たいに宿題依頼のテンプレが必要なんでね？ 見たところ、宿題はいまだF77中心のようだし。ｆ９０かｆ７７かも選択してもらわないと。 646 名前：さおり [2006/07/14(金) 16:17:43 ] そんなこと絶対しません。よろしくおねがいします。 647 名前：デフォルトの名無しさん mailto:sage [2006/07/14(金) 18:27:20 ] つーことで、さおりタンの仕事は 1. 微分方程式の確認 2. 授業でならった微分方程式の数値解析法の(せめて)名前の確認 だね。 648 名前：デフォルトの名無しさん mailto:sage [2006/07/14(金) 19:44:38 ] >>646 「そんなことしません」じゃなくて、レポートとして期待されているプログラムを 作成するためには情報が足りねーの。 課題を入手してきなさい。あと、おっぱいうｐな。 649 名前：デフォルトの名無しさん [2006/07/14(金) 20:48:58 ] >>637に従って解いたｗ 微分方程式は、一番原始的な方法で積分したｗ　 大学１，２生ならこれでも合格ラインを超えるだろう。３，４年生ならルンゲ・クッタを使うべきだろう。 原点（０，０）から初速ｖｘ、ｖｙで玉を放り投げる。地面より下に座標がきたら（ｙ＜０．０）終了する。 PROGRAM yamajun PARAMETER( g = 9.8, fric = 0.3 ) t = 0.0 dt = 1.0 / 8.0 amass = 1.0 x = 0.0 y = 0.0 vx = 10.0 vy = 10.0 ax = 0.0 ay = 0.0 1 CONTINUE WRITE(9, *) x, y, t !, vx, vy !, ax, ay IF ( y < 0.0 ) GOTO 999 t = t + dt x = x + dt * vx y = y + dt * vy ax = - fric * vx * SQRT(vx**2 + vy**2) / amass ay = -g - fric * vy * SQRT(vx**2 + vy**2) / amass vx = vx + dt * ax vy = vy + dt * ay GOTO 1 999 STOP END FORTRANは、うほっの世界であるべきなので、さおりは次回からは カヲルという少年に性転換すること。 650 名前：637 mailto:sage [2006/07/14(金) 21:17:58 ] >>649 thx. いやー。数値解析の教科書とか見て気が重かったのだが、ありがとう。 さおりタン満足してくれるといいね。 651 名前：649 mailto:sage [2006/07/14(金) 21:40:27 ] >>650 >>649はあまりに手抜きｗ　誰かもっと良いの作ってやれｗ ２元２階の微分方程式なので、本のルンゲ・クッタルーチンにそのまま放り込めないｗ ちょっと脳みそ使う必要がある。数独１回分くらいｗ　 今日暑いから考えらんねー　冷房ほしいよー 652 名前：デフォルトの名無しさん [2006/07/14(金) 21:43:43 ] 2,3,5.6,2.7,2,-,9,3 という「-」のような記号の混じった数列データが入ったファイルがあって、 それを下のように読み込みたいのですが、 real*8 a,b,c,d,e,f,g,h read(10,*) a,b,c,d,e,f,g,h これで「-」を読み込むときエラーが出て終了します。 「-」が出てきた場合には-1.0d-10などの任意の数字を代入しておく ような方法はないでしょうか？ よろしくお願いします。 653 名前：デフォルトの名無しさん [2006/07/14(金) 23:37:38 ] >>652　ほい。　F９０で書いたが、主要なところはF７７でもおｋ。　質問あればおｋ。 やり方は色々あるが、基本的にはまずデータを１行文字列として読み込んで、 カンマで区切って１こづつ数値として読み取るのが標準だろう。 ここでは読み取りエラーが起きてから例外処理的に”−”を処理しているが、 ”−”が多いならはじめから一度文字列で読み込んで処理するほうがすっきりする。 PROGRAM vip IMPLICIT NONE CHARACTER (LEN = 136) :: text REAL(KIND = 8) :: x(8) INTEGER :: i, ipos, io DO READ(9, *, IOSTAT = io) x IF (io == -1) EXIT IF (io > 0) THEN BACKSPACE(9) READ(9, '(A)') text DO i = 1, 8 ipos = INDEX(text, ',') - 1 IF (ipos == -1) ipos = LEN(text) ! 8th number IF (TEXT(ipos:ipos) == '-') THEN x(i) = -1.0d-10 ELSE READ(TEXT(1:ipos), *) x(i) END IF text = text(ipos + 2:) END DO END IF PRINT '(8g9.2)', x END DO STOP END PROGRAM vip 654 名前：６５３ [2006/07/14(金) 23:39:00 ] ■入力データ 2,3,5.6,2.7,2,1,9,3 2,3,5.6,2.7,2,-,9,3 2,3,5.6,2.7,2,2,9,3 ■出力 2.0 3.0 5.6 2.7 2.0 1.0 9.0 3.0 2.0 3.0 5.6 2.7 2.0 -0.10E-09 9.0 3.0 2.0 3.0 5.6 2.7 2.0 2.0 9.0 3.0 Press any key to continue 655 名前：デフォルトの名無しさん [2006/07/15(土) 03:04:03 ] 　　　　　 |　|/ノ二__‐──ｧ　　 ヽﾆニ二二二ヾ } ,'⌒ヽ 　　　　 /⌒!|　 =彳ｏ｡ﾄ￣ヽ　　　　　'´ !o_ｼ｀ヾ　| i/ ヽ !　　　　ウホッ！　いい男・・・ 　　　　 ! ハ!|　　ー─　'　　i　　!　 　 `' 　 '' " 　 ||ヽ　l | 　　　　| | /ヽ!　　　　　　　　|　　　　　　　　　　　 |ヽ i　! 　 656 名前：652 mailto:sage [2006/07/15(土) 09:43:09 ] >>653 ウホッ！ 657 名前：デフォルトの名無しさん mailto:sage [2006/07/15(土) 09:44:53 ] >>623,>>629,>>634,>>646 >>637をもとに>>649が解いたものは、m、k、初期位置、時間の刻みに 適当な値を代入して解いている。本来はあんたが提示すべき情報。 >>630,>>647,>>648を良く読め。 658 名前：デフォルトの名無しさん mailto:sage [2006/07/15(土) 10:57:31 ] 教えて！goo にも同じような質問が・・・ 659 名前：637 mailto:sage [2006/07/15(土) 12:20:32 ] >>658 じゃあ我々の役目も終ったのかな？ みんなお疲れ様。 660 名前：デフォルトの名無しさん [2006/07/15(土) 17:35:26 ] 女が来るとふいんきが乱れるから、ホモOnlyスレにしないか？ 661 名前：デフォルトの名無しさん mailto:sage [2006/07/15(土) 18:59:43 ] プログラムの実行中に一瞬ポーズする関数ってありますか？ 662 名前：デフォルトの名無しさん mailto:sage [2006/07/15(土) 19:17:53 ] >>661 PAUSE 663 名前：デフォルトの名無しさん mailto:sage [2006/07/15(土) 19:30:41 ] >>662 ただPAUSEはFORTRAN2003あたりで廃止勧告が出されている。 TSS利用で便利だから漏れは好きなんだけどなｗ 664 名前：デフォルトの名無しさん [2006/07/15(土) 20:28:29 ] write(*,*) abs(-16.8-(6.8))>10.0 を出力すると「T」になってしまいます。 助けてください。 665 名前：664 [2006/07/15(土) 20:32:23 ] 間違えました、 real*8 a,b a=-16.80000000000000 b=- 6.80000000000000 write(*,*) abs(a-b)>10 とやったら「T」と出てしまいます。 何が原因でしょうか？ 666 名前：664 mailto:sage [2006/07/15(土) 20:53:41 ] 更に状況を書かせて頂きますと、 読み取りファイルに -16.8,-6.8 とあって、 real*8 a,b open(10,file="test") read(10,*)a,b write(*,*) abs(a-b)>10 とやったら「T」と出力されるのです。 667 名前：デフォルトの名無しさん mailto:sage [2006/07/15(土) 20:59:25 ] FORTRAN77は大文字でプログラム書くのがデフォルトなのか 668 名前：デフォルトの名無しさん [2006/07/15(土) 21:04:15 ] >>664 それは良くあることｗ　気にするな。　浮動小数点之罪。 669 名前：デフォルトの名無しさん [2006/07/15(土) 21:06:04 ] >>667 本来のFORTRAN77の規格では大文字のみ。 そもそも、昔は小文字が表示できない端末とか良くあったｗ 670 名前：デフォルトの名無しさん [2006/07/15(土) 21:22:32 ] >>664 参考までに、INTEL　FORTRANの結果を見てくれ。 PROGRAM unko IMPLICIT NONE REAL (KIND = 8) :: a, b a = -16.8 b = -6.8 PRINT *, a, b, ABS(a - b) < 10 a = -16.8d0 b = -6.8d0 PRINT *, a, b, ABS(a - b) < 10.0d0 STOP END PROGRAM unko ■結果 -16.7999992370605 -6.80000019073486 T -16.8000000000000 -6.80000000000000 F Press any key to continue この場合は倍精度（８バイト）の変数に単精度（４バイト）の常数を与えているので、 長さの足りない４バイトにゴミが入っている。 ■ とにかく浮動小数点の誤差とか桁落ちとか実数の比較とかでググレば説明が出るはず。 この辺は昔は電子計算機利用の一番最初に注意されて来た事柄だ。 671 名前：664 mailto:sage [2006/07/15(土) 21:45:51 ] >>670 どうもありがとうございます。 一応計算中の変数を画面に表示させてチェックしたのですが、 aは-16.8000000000000 bは -6.8000000000000 と表示されていました。 この上でabs(a-b)>10をやるとTになりました。 672 名前：デフォルトの名無しさん mailto:sage [2006/07/15(土) 22:03:21 ] >>671 (T_T) 30桁くらい出力してごらん。 write ( *, '(2f35.30)' ) a, b 673 名前：670 [2006/07/16(日) 00:28:13 ] >>671 漏れの書き方が誤解を招いたらしい。すまんこ。 数直線のうちFORTRANであらわされるのは、０を中心として正負にほぼ同じだけの 範囲の中の、有限個の有理数だけである。数直線の上に、定規の目のように切った点が あって、定規の目に挟まれた部分は近くの目で代表することになる。 （浮動小数点のばあい、この定規の目は体育で使う巻尺のように、０の近傍は細かく 区切ってあるが、遠くへ行くほど荒く区切るようになっている。） コンパイラにも依存するのだが、基本的には単精度や倍精度の浮動小数点は、 ２の（負を含む）べき乗の和であらわされるような有理数になっている。 このため、１０進表記で書いた有限小数が、かならずしも数直線の上に切った目に 一致するとは限らない。 つまり １．５＝１＋０．５＝２＾０＋２＾−１ ０．７５＝０．５＋０．２５＝２＾−１＋２＾−２ これらの数は、きっかり割り切れる（すなわち、定規の目の上にのっている） ところが ０．８＝０．５＋０．２５＋０．０６２５＋・・・・＝２＾−１＋２＾−２＋２＾−４＋・・・・・ となって、数直線の上に切った目の上に載っていない。 つまり、近くの別の数に四捨五入されている。 >>672のように、FORMATで小数点いかの数を長く出すと、以下のようになる。 -16.7999992370605 -6.80000019073486 T -16.8000000000000 -6.80000000000000 F -16.800000000000000710542735760100 -6.799999999999999822364316059975 つまり、倍精度の−１６．８は、実際は-16.800000000000000710542735760100・・・・・・ という数に四捨五入されている。−６．８も同様。 （ただし、有効桁からすれば、＊フォーマットで出力される分しか制度が無い。 つまり＊フォーマットで出力される数の最後の桁程度の範囲で四捨五入がされている。） 674 名前：デフォルトの名無しさん [2006/07/16(日) 14:38:42 ] 　積分区間の下限Aと上限B、分割数N（偶数）を読み込んで、 Simpson公式を用いて、定積分I=(A,B)∫f(x)dxを計算するプログラムを作れ。 　ただし、数値積分の部分のプログラムは、任意の関数に対して使用できる サブルーチン・サブプログラムを用いよ。その場合関数f(x)の計算は関数サブプログラムで行う。 　さらに上記の分割数Nが奇数に対しても使えるように拡張せよ。 Simpson公式を用いて、Nが奇数のときは1区画だけ台形公式を用いるようにすればよい。 手も足も出ません…、お願いできますか？ 675 名前：デフォルトの名無しさん mailto:sage [2006/07/16(日) 16:37:20 ] >>674 PROGRAM HOGE READ(*,*) A, B, N IF (MOD(N,2) .NE. 0) GOTO 9999　… ※１ C DX = (B - A) / N S = 0. DO 10 I = 1, N/2 X1 = A + 2*(I-1)*DX X2 = A + (2*I - 1)*DX X3 = A + 2*I*DX S = S + (F(X1) + 4*F(X2) + F(X3)) 10 CONTINUE S = S * DX / 3. C IF (MOD(N,2) .NE. 0) THEN　… ※２ C X1 = B - DX C X2 = B C S = S + (F(X1) + F(X2)) * DX / 2 C END IF WRITE(*,*) S 9999 CONTINUE STOP END 前半部分はそのまま。（※２の所はいらない） 後半部分は※１のIF文を消して※２のコメントアウトされてるIFブロックを追加する。 676 名前：デフォルトの名無しさん [2006/07/16(日) 21:55:10 ] ファイルを読み込んだり作ったりするとき、 実行ファイルと同じディレクトリの中に 別のディレクトリを作って保存したりするにはどうしたら良いでしょうか。 677 名前：デフォルトの名無しさん [2006/07/16(日) 22:07:35 ] >>676 既出の質問だが・・・・お答えしよう。 FORTRANの標準命令ではできない。 コンパイラメーカー固有の拡張命令でシステムの関数を呼び出すものがあるだろうから それらを通じてやるしかない。マニュアルの後ろのほうの章を探して味噌。 678 名前：676 mailto:sage [2006/07/16(日) 23:32:02 ] >>677 ｻﾝｸｽｺ 679 名前：デフォルトの名無しさん mailto:sage [2006/07/17(月) 12:21:13 ] 番組の途中ですが・・・・・・・・・・・・ 盧武鉉（ノ・ムヒョン）大統領が今月11日に行われた与党ヨルリン・ウリ党の指導部 および国会の統一外交通商委員会に所属する議員らとの晩さん懇談会で行った 発言が波紋を呼んでいる。 一部新聞は懇談会出席者の証言を引用し、盧大統領は「ブッシュ米大統領が 北朝鮮問題を善と悪の対立概念で見ているため、説得が難しくなっている。 米国は友邦なので厳しく責めることは出来ないが、 日本とは対決しなければならない」と語ったという。 ＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝ ソース:朝鮮日報 japanese.chosun.com/site/data/html_dir/2006/07/17/20060717000017.html 680 名前：デフォルトの名無しさん [2006/07/18(火) 05:15:28 ] 多変数の4次のルンゲクッタのプログラムを書いているのですが、どうもうまくいきません。 今、厳密解を y(1) = exp(x) + exp(-x) y(2) = exp(x) - exp(-x) とおいて、その微分方程式 dydx(1) = exp(x) - exp(-x) dydx(2) = exp(x) + exp(-x) を数値的に解いて、厳密解と比較しています。 サブルーチンのderive中で、 dydx(1) = exp(x) - exp(-x) dydx(2) = exp(x) + exp(-x) とおいた場合は、きちんと数値解と厳密解が一致するのですが、 dydx(1) = y(2) dydx(2) = y(1) とおくと、両者は一致しません。これはなぜでしょうか？ 考えているのですが、一向にわかりません。 プログラムは以下になります。よろしくお願いします。 681 名前：デフォルトの名無しさん [2006/07/18(火) 05:17:34 ] program rk4 implicit none double precision x0, y0, x, h, y(2), dydx(2), yout(2) external derive h = 1.0d-2 x0 = 0.0d0 y(1) = 2.0d0 y(2) = 0.0d0 do x = x0, 2.0d0, h write(*, *) x, dabs((dexp(x) + dexp(-x)) - y(1)) write(*, *) x, dabs((dexp(x) - dexp(-x)) - y(2)) write(*, *) call derive(x, y, dydx) call rk4(y, dydx, x, h, yout, derive) y(1) = yout(1) y(2) = yout(2) end do stop end subroutine derive(x, y, dydx) implicit none double precision x, y(2), dydx(2) c dydx(1) = y(2) c dydx(2) = y(1) dydx(1) = dexp(x) - dexp(-x) dydx(2) = dexp(x) + dexp(-x) return end

---

---

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

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