FORTRANのプログラムです。

質問

FORTRANのプログラムです。
下に記載したプログラムを、行列値の計算にするプログラムに改良するにはどうすればよいでしょうか？
アドバイスお願いします。

c ***************************
if(root.LT.x0) GOTO 11
x0=root; x1=xs0; x2=xs1
GOTO 99
11 x0=root; x1=xs0; x2=xs2

GOTO 99
c

c ***************************
OPEN(6,FILE='sum.dat')

1002 FORMAT(E12.5,3x,4(E12.5,4x))

CLOSE(6)
STOP

100 END

c *************************
SUBROUTINE FRES(x,F)
DOUBLE PRECISION x,F

c f=3.*x+DSIN(x)-DEXP(x)
f=(x-1.)*(x-2.)*(x-3.)*(x-4.)*(x-5.)*(x-6.)

RETURN
END

よろしくお願いいたします。

回答No.2

f272

補足によってプログラムの方は分かるけど，

| x-1 3 2-x |
| 5 x-2 x^3-4 |
| 3 x-1 x-x^2 | = 0
| 3 5 x^2-4 |

これが何を意味しているのか分かりません。見たところ行列式を意味しているようにもみえるが4x3行列に対しては無意味です。

回答No.1

f272

「下に記載したプログラムを、行列値の計算にするプログラムに改良するにはどうすればよいでしょうか？」
と言っているが，記載されているのはプログラムの断片にすぎないし，そのうえ行列値に計算にするプログラムというのが何のことか意味不明です。もう少し質問を明確にしてもらわなければ，どうしようもありません。

補足

あいまいな質問で申し訳ございません。

プログラム全文です。

DOUBLE PRECISION f0,f1,f2,x0,x1,x2,h1,h2,gam
DOUBLE PRECISION aa,bb,cc,f,root,b1,b2,BT
DOUBLE PRECISION xs0,xs1,xs2,EXPS,umin,umax

EXPS=1.e-3

umin=0.1; umax=5.5; x0=3.3

x2=umin; x1=umax

99 xs0=x0; xs1=x1; xs2=x2

h1=x1-x0; h2=x0-x2; gam=h2/h1

CALL FRES(x0,f0); CALL FRES(x1,f1); CALL FRES(x2,f2)

aa=(gam*f1-f0*(1.+gam)+f2)/(gam*h1**2.*(1.+gam))
bb=(f1-f0-aa*h1**2.)/h1
cc=f0

b1=bb+DSQRT(bb**2.-4.*aa*cc)
b2=bb-DSQRT(bb**2.-4.*aa*cc)
if (bb.GT.0.) BT=b1
if (bb.LT.0.) BT=b2
if (bb.EQ.0.) BT=2.*x0-x1
root=x0-2.*cc/BT

write(*,1002) x0,x1,x2
write(*,1002) f0,f1,f2
write(*,1002) h1,h2,gam
write(*,1002) aa,bb,cc,root

if(DABS((root-x0)/root).LT.EXPS) GOTO 100

if(root.LT.x0) GOTO 11
x0=root; x1=xs0; x2=xs1
GOTO 99
11 x0=root; x1=xs0; x2=xs2

GOTO 99
c

OPEN(6,FILE='sum.dat')

1002 FORMAT(E12.5,3x,4(E12.5,4x))

CLOSE(6)
STOP

100 END

SUBROUTINE FRES(x,F)
DOUBLE PRECISION x,F

f=(x-1.)*(x-2.)*(x-3.)*(x-4.)*(x-5.)*(x-6.)

RETURN
END

-------------------------------------------------------------------
このプログラムを、

| x-1 3 2-x |
| 5 x-2 x^3-4 |
| 3 x-1 x-x^2 | = 0
| 3 5 x^2-4 |

上式のxの値を求めるプログラムに改良したいです。