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

くだすれFORTRAN（超初心者用）その4at TECH

くだすれFORTRAN（超初心者用）その4 - 暇つぶし2ch■コピペモード
□スレを通常表示
□オプションモード
□このｽﾚｯﾄﾞのURL
■項目テキスト

350:デフォルトの名無しさん
 09/06/21 10:45:18 
ループのところ 
do i=1, n 
  b(i) = 0 
  do k=1,n 
    b(i) = b(i) +M(i, k) * a(k) 
  end do 
end do 
 
直接関係ないけど、Mの行と列を入れかえたほうがいいかも。

351:346
 09/06/21 18:17:49 
>>350さん 
ありがとうございました。 
うまくいきました！

352:デフォルトの名無しさん
 09/06/23 11:08:29 
PROGRAM GAUSS_SOLVER 
IMPLICIT NONE 
 
REAL,DIMENSION(100,100)::A    !行列A（最大100元まで解ける） 
REAL,DIMENSION(100)::B,X      !右辺Bと解Xのベクトル 
INTEGER::N                    !元数 
INTEGER::I,J,K                !制御変数 
 
WRITE(*,*)'このプログラムは連立一次方程式AX=Bのい解Xを求めます' 
WRITE(*,*)'AandB.TXTの入力フォーマットは次のようです' 
WRITE(*,*)'たとえば,３元の連立方程式の場合' 
WRITE(*,*)'3               (元数を与えてください)' 
WRITE(*,*)'A11 A12 A13 B1  (データは空白やタブで別けてください)' 
WRITE(*,*)'A21 A22 A23 B2  (AIJとBIには実数を与えてください)' 
WRITE(*,*)'A31 A32 A33 B3' 
WRITE(*,*)'解XはX.TXTに出力されます' 
 
OPEN(1,FILE='AandB.TXT') 
READ(1,*) N 
DO I=1,N 
  READ(1,*) (A(I,J),J=1,N),B(I) 
    END DO 
CLOSE(1) 
 
DO I=1,N-1 
  DO J=I+1,N 
    DO K=I+1,N 
      A(J,K)=A(J,K)-A(J,I)*A(I,K)/A(I,I) 
     END DO 
     B(J)=B(J)-B(I)*A(J,I)/A(I,I) 
 END DO 
END DO

次ページ