Fortranのif文で出力が5になってしまう問題
- 以下のFortranのコードで、xとthetaに5,5と入力しても、出力が5となってしまいます。
- 500になってほしかったのですが、どこが悪いのか教えてください。
- 問題箇所はif文の条件判定とその後の処理です。正しい条件を設定していないため、出力が5になってしまいます。
- ベストアンサー
fortran if文
以下のコードで、xとthetaに5,5と入力しても、出力が5となってしまいます。500になってほしかったんですが・・・ どこが悪いのか教えて下さい。 ----------------------------------------------- read(5,*)x,theta if(0.0d0.LE.x.LE.4.0d0) then eta=theta go to 20 else if(4.0d0.LT.x.LT.7.0d0) then eta=100.0d0*theta go to 20 else if(7.0d0.LE.x.LE.10.0d0) then eta=sqrt(3.0d0*sqrt(2.0d0))*theta go to 20 end if 20 write(6,30) eta 30 format(F10.5) stop end -----------------------------------
- BOY12345
- お礼率12% (20/162)
- その他（プログラミング・開発）
- 回答数1
- ありがとう数1
- みんなの回答 （1）
- 専門家の回答
質問者が選んだベストアンサー
read(5,*)x,theta ! you cannot write 0 < x < 4 ! you need to write 0 < x and x < 4 ! to get intended behavior ! the two codes have the different meanings. ! tested with g95 compiler if(0.0d0 .LE. x .and. x .LE.4.0d0) then eta=theta go to 20 else if(4.0d0.LT. x .and. x .LT.7.0d0) then eta=100.0d0*theta go to 20 else if(7.0d0.LE. x .and. x .LE.10.0d0) then eta=sqrt(3.0d0*sqrt(2.0d0))*theta go to 20 end if 20 write(6,30) eta 30 format(F10.5) stop end
関連するQ&A
- fortran if文
fortranのif文で、以下の様に場合わけしてプログラムしたのですが、 このようなエラー文が出ます。どう対処したらいいのか教えて下さい。 Error: An unterminated block exists. プログラム↓ if(-0.5d0.LE.theta.LE.0.0d0) then eta=2 if(0.0d0.LT.theta.LT.0.5d0) then eta=5 if(0.5.LE.theta.LE.1.0d0) then eta=3 場合わけの範囲指定の仕方が不味いのではないかと思っていますが、 どうしていいものかわからず・・ 皆様のお力を貸して下さい。
- ベストアンサー
- その他（プログラミング・開発）
- fortran77でのif文
FORTRAN77の問題の質問です。 授業でこんな問題が出て困ってます。 あるクラスにN人います。 そのクラスの平均身長がXです。 そのN人の中で、X以上の身長の人数Aを印刷するプログラムを作れと言う問題です。 ここで、IF～THENを使いなさい。 分かる人がいたら教えてください。 PROGRAM TEST7 WRITE(*,*)\'x人分の身長データを入力 READ(5,*)x WRITE(*,*)\'身長のボーダーラインを入力\' READ(5,*)G A=A-G B=B-G C=C-G D=D-G E=E-G F=F-G IF(A.GT.0) THEN STOP END 途中までは、こんな感じかと思ったのですが、行き詰ってしまいました。
- 締切済み
- その他（プログラミング・開発）
- FORTRAN77のプログラミングについて教えてください。
・学校の課題で任意の数値の平均と、入力した数の個数を出力するプログラムを作 りなさいという課題が出たのですがうまくいきません。次にあげるプログラムの 中で修正する部分を教えてください。 *課題2 平均値 integer N real X,T1,M * 初期値 N=0 T1=0 * 累積 10 read(5,*,end=20) X N=N+1 T1=T1+X go to 10 * まとめ 20 if(N.GE.2) then M=T1/N write(6,*) '総数は' ,N, '平均は',M ELSE IF(N.EQ.1) THEN WRITE(6,*) '総数は',N, '平均は',T1 ELSE WRITE(6,*) '数値がない' END IF END IF END ・この問題は最初からよくわかりません。教えてください。 「3つの数a,b,cを読み込む。a,b,cを三角形の３辺の長さとしたとき三角形にな るかを判定しなさい。三角形にならない場合はその面積をヘロンの公式を求め て表示する。 s=(a+b+c)/2 , ss=s(s-a)(s-b)(s-c) として S=sqrt(ss) とする。」 以上の二つです。分かる方お願いします。
- 締切済み
- その他（プログラミング・開発）
- 二分法(FORTRAN)
大学の情報処理演習で、FORTRANのプログラミングをやっています。 で、二分法の問題なんですが、X^3-X-1の解で、[0.2] にあるものを計算するプログラムを作る課題が出されたのですが、うまくいきません。 そこで、自分の作ったプログラムを添削してください。 program nibunnhou1 c real x1,x2,c,x,f,r f(x)=x**3-x-1 c x1=0.000 x2=2.000 r=f(x1) c write(6,*)'start' c do n=1,100 c=(x1+x2)/2 if(r*f(c).gt.0) then x1=c else x2=c end if c if(f(c).lt.0.0000) go to 1 continue end do c 1 write(6,*)'ans=' write(6,*) c c stop end です。これだと、答えが「１．」とかになってしまいます。 お願いします。
- ベストアンサー
- その他（プログラミング・開発）
- fortranのif文で・・・
フォートランの if文で IF(M-1) 15,15,12 ってのがあるんですが，これはどういう意味を表しているのでしょうか ？ IF(R.LE.0.D0) THEN ・・・・ ならＲが0より小さいときは・・・・しなさい．って命令ですよね． でもうえのは，カッコの中がいったいどういうときならってのを表してるのかも分からないし，15，15，12 ってのも意味がわかりません． 基本的なことかもしれませんが，よろしくお願いします． 一応その部分のプログラムも載せておきます． X(N+1)=X(1) Y(N+1)=Y(1) DO 10 I=1,N XM(I)=(X(I)+X(I+1))/2 10 YM(I)=(Y(I)+Y(I+1))/2 IF(M-1)15,15,12 12 XM(NC(1))=(X(NC(1))+X(1))/2 YM(NC(1))=(Y(NC(1))+Y(1))/2 DO 13 K=2,M XM(NC(K))=(X(NC(K))+X(NC(K-1)+1))/2 13 YM(NC(K))=(Y(NC(K))+Y(NC(K-1)+1))/2 15 CONTINUE
- ベストアンサー
- その他（プログラミング・開発）
- fortranのif文にでるエラーがわかりません。
c termination test sorce=amax1(resorm,resoru,resorv) 310 if((niter.eq.20.and.sorce.gt.1.0e7*sormax) & .or.(niter.eq.maxit)) then write(6,320) 320 format & (1h ,'solution diverges or needs more iterations'/) if(sorce.gt.sormax) go to 300 c----- end of iterations 文献に載っていたプログラムのメインプログラムの一部です。 ３１０のif文に対してでるエラーの意味が どうしてもわかりません。 fortranのコンパイルのエラーに関してはだいたいはわかっているつもりなのですが・・・ 理由がわかるかたいらしたら教えていただけませんか？ よろしくお願いします。
- ベストアンサー
- その他（プログラミング・開発）
- fortran 繰り返しループから抜きたい
10 continue ： ： ： do 100 i=1,10 if (abs(x (i)-y(i)) .lt. 1.0) then r= r + dr goto 100 else goto 1000 end if 100 continue ： ： ： 1000 contiue end すべての条件（ｉ＝1～10）が満足されたら、ループから抜きたいのですが、出来ませんでした。 教えて頂けませんでしょうか？よろしくお願しますね！*<>*！
- 締切済み
- 数学・算数
- fortran 固有値を求めるプログラム
２x２の実行列の固有値を求めるモジュール関数の中で、d > 0 のときに出てくる sign(squt(d), -b)の値が示す意味と eval(2) = cmplx(c/e,0.0d0) でなぜ c/eの値が入るのかがわかりません。教えて下さい。よろしくお願いします。 プログラム module subprogs implicit none contains function eval2x2mat(a) result(eval) real(8), intent(in) :: a(:,:) complex(8) eval(2) real(8) b, c, d, e if (size(a,1) /= size(a,2)) stop ' not square ' if (size(a,1) /= 2) stop ' not 2x2 matrix ' b = -0.5d0*(a(1,1)+a(2,2)) c = a(1,1)*a(2,2)-a(1,2)*a(2,1) d = b**2-c if ( d < 0.0d0 ) then eval(1) = cmplx(-b,sqrt(-d)) eval(2) = conjg(eval(1)) else if ( d > 0.0d0 ) then e = -b+sign(sqrt(d),-b) ←ここの部分 eval(1) = cmplx(e,0.0d0) eval(2) = cmplx(c/e,0.0d0) else ↑ここの部分 eval(1) = cmplx(-b,0.0d0) eval(2) = eval(1) endif end function eval2x2mat end module subprogs program main use subprogs implicit none real(8), allocatable :: a(:,:) integer :: n = 2 allocate(a(n,n)) a(1,1:2) = (/-1,1/) a(2,1:2) = (/-1,-1/) write(*,*) a(:,:), eval2x2mat(a) end program main
- ベストアンサー
- その他（プログラミング・開発）
- FORTRAN77
助けてください～ --------------------------------------- 華氏゜Ｆが与えられた時、摂氏℃と絶対温度Ｚを求めよ C=(5/9)(F-32),Z=C+273.16 ここで華氏゜Ｆは書式カードF1.5でカードに用意されていて、数値999.9のデータを読み込んだら終わりとする。 ただし、華氏゜Ｆの値は、-50.0゜Fから200.0゜Fまで10.0゜Fごとに変化させる。 --------------------------------------- で、私が頑張って(っていっても教科書を見ながら)書いたのが以下です。 --------------------------------------- * write(6,200) 200 format(1h ,8x,1hf,10x,1hc, 9x,1hz) 10 read(5,100)f 100 format(f5.1) if(f.eq.999.9)stop c=(f-32.0)*5.0/9.0 z=c+273.16 write(6,201)f,c,z 201 format(1h ,3(5x,f5.1)) go to 10 end ------------------------------------- 100 formatの下からdo文を使うのが狙いかと思いますがよくわかりません。 ちなみに参考教科書は入力データが ------------------------------------- 32.5 135.5 68.3 999.9 ------------------------------------- となっていました。多分これは使わないです。 よろしくお願いします(><)
- 締切済み
- その他（プログラミング・開発）
- Fortran のエラー
9002700385647,2,0,1,0,915110373,20050524,0,2,4911 9004300501590,1,0,1,0,910410489,20050630,-1,1,4527 このようなデータ1万人分必要なところだけ読み取って 書き出そうとするのですがうまくいきません。 This application has requested the Runtime to terminate it in an unusual way. とエラーメッセージが出ます。なぜでしょうか。 ――――――――――――――――――――――――――――――― OPEN (1,FILE='OKAMOTODA.TXT',STATUS='OLD') OPEN (2,FILE='NEWNEW.TXT',STATUS='NEW', 1 FORM='FORMATTED') II=0 31 READ(1,*,END=900)K1,K2,K3,K4,K5,K6,K7,K8,K9,K10 II=II+1 WRITE(2,991) II,K1,K2,K3,K4,K5 991 FORMAT(2X,I3,5I20) IF(II .LT. 10000) THEN GOTO 31 ELSE GOTO 900 ENDIF 900 CLOSE(1) CLOSE(2) END ―――――――――――――――――――――――――――――――
- 締切済み
- その他（プログラミング・開発）
お礼
できました♪ ありがとうございます。 if(0.0d0 .LE. x .and. x .LE.4.0d0) then eta=theta go to 20 else if(x .LT.7.0d0) then eta=100.0d0*theta go to 20 else if(x .LE.10.0d0) then eta=sqrt(3.0d0*sqrt(2.0d0))*theta go to 20 end if でもいけることがわかりました！！ ありがとうございます☆