●●fortran90で二分法を●●

1 ：ルチャ ：2001/02/25(日) 03:19

FORTRAN90を覚えたての者です。
二分法を使って方程式の解を求めるプログラムを
大至急作らねばならなくなったのですが、
やっと見つけたサンプルプログラムがBASICで書かれていて困ってます。

どなたかこのサンプル↓↓をFORTRAN90に書き換えてくださいませんか。
77でもいいのですが初心者なのでなるべく分かりやすくしてもらえると助かります。
明日本屋に行って探してからでは間に合いそうにありません。
ほんと勝手ですけど助けてください。わらにもすがる思いです。(TΔT)

2 ：ルチャ ：2001/02/25(日) 03:19

DEF FNY(X)=X^3+X-1
PRINT "X^3+X+1=0 の二分法による数値計算"
A=0 ;B=1
LOCATE 10, 3:PRINT"初期値a=" ;;PRINT USING "##;A
LOCATE 10, 4:PRINT"初期値b=" ;;PRINT USING "##;B
GOSUB *NIBUN
LOCATE 15, 6:PRINT"近似解X=" ;;PRINT USING "##.###";C
END

*NIBUN: '二分法のサブルーチン
　　*AA:C=(A+B)/2
IF FNY(C)*FNY(A)<0 THEN B=C ELSE A=C
IF ABS(A-B)>.0001 THEN *AA
RETURN

3 ：名無しさん＠お腹いっぱい。 ：2001/02/25(日) 03:32

プログラム板とか理系板とか数学板とかその辺に行った方がいいのでわ。

4 ：F77 ：2001/02/25(日) 03:47

*234567
* 2-BUN-HOU
*
REAL F,X,EPS,A,B,C,Y
PARAMETER (EPS=1.0E-06)

F(X)=X**3+X-1
WRITE(6,*) '2-BUN-HOU NI YORU SUCHI KEISAN'
A=0
B=2
WRITE(6,'(A,F6.3)') 'SHOKICHI:A= ', A
WRITE(6,'(A,F6.3)') 'SHOKICHI:B= ', B
Y=F(A)

10 CONTINUE
C=(A+B)/2
IF(Y*F(C).LT.0) THEN
B=C
ELSE
A=C
END IF
IF(ABS(A-B).GE.EPS) GO TO 10

WRITE(6,'(A,F9.6)') 'KINJIKAI:X= ', C
END

5 ：F77 ：2001/02/25(日) 03:49

ずれてる

6 ：名無しさん＠お腹いっぱい。 ：2001/02/25(日) 04:09

http://mentai.2ch.net/test/read.cgi?bbs=prog&key=983038489&to=3

7 ：ルチャ ：2001/02/25(日) 04:10

>>5さん
(TΔT)ありがとうございます。
金払いたい気分です。
これがあれば自分のプログラムが組めそうです。
最後に参考書を見ても分からなかった部分を質問していいでしょうか。

�@　WRITE(6,'(A,F9.6)')←カッコ内の右側部分の意味
�A　最後の方のIF文に出てきたABSの意味
です。

8 ：F77 ：2001/02/25(日) 04:18

�@FORMAT文を中に入れる書き方をしただけ
�A絶対値
あと、>>4で変数Yはいらなかった。

# マルチポストはやめようね。

9 ：ルチャ ：2001/02/25(日) 04:24

>>F77さん
どうも御迷惑おかけしました　(TΔT)
お礼しようにもできないので、
いつか誰かが困っていたら今日の恩返しと思って助けてあげようと思います。

10 ：名無しさん＠お腹いっぱい。 ：2001/02/25(日) 05:24

あんた学生でしょ？。レポートは自分でやりましょう。

11 ：名無しさん＠お腹いっぱい。 ：2001/02/25(日) 08:16

マルチポストあげ
http://mentai.2ch.net/test/read.cgi?bbs=prog&key=983038489
http://piza.2ch.net/test/read.cgi?bbs=tech&key=983038656
http://cocoa.2ch.net/test/read.cgi?bbs=unix&key=983038745
http://mentai.2ch.net/test/read.cgi?bbs=infosys&key=983037683
http://mentai.2ch.net/test/read.cgi?bbs=sim&key=983037987
http://cheese.2ch.net/test/read.cgi?bbs=sci&key=983037846
http://cheese.2ch.net/test/read.cgi?bbs=math&key=983038350

12 ：名無しさん＠お腹いっぱい。 ：2001/02/25(日) 13:32

>>10
まあまあ、この時期は内定が出たのに単位が足りなくて焦ってる学生もいるだ
ろうからそう言わずに。

read.cgi ver5.26+ (01/10/21-)