[轉錄][大師] 程式設計苦惱

[f94 (→罪ㄞˋ 仁甫ａ婷←)]

※ [本文轉錄自 f94 信箱]

作者: [e-mail] ([e-mail])
標題: [大師] 程式設計苦惱
時間: Sat May  3 13:15:51 2008

作者: f94 (→罪ㄞˋ 仁甫ａ婷←) 看板: IamMaster
標題: [大師] 程式設計苦惱
時間: Sat May  3 13:15:07 2008

題目

4X+Y-Z=5
-X+3y+Z=-4
2X+2Y+5Z=1



C++++++++++++++++++++  NUM=3 IN THIS EXAMPLE +++++++++++++++++++++
C+++ COEFFICIENTS: C(1:NUM 1:NUM),CONSTANT TERMS:C(1:NUM,NUM+1)+++
C++++++++++++DIMENTION OF INPUT:C(NUM+1,NUM+1)++++++++++++++++++++
c++++++++++++++DIMENTION OF OUTPUT: S(NUM) +++++++++++++++++++++++
        DIMENSION C(4,4),S(3)
        C(1,1)=4
        C(1,2)=1
        C(1,3)=-1
        C(1,4)=5
        C(2,1)=-1
        C(2,2)=3
        C(2,3)=1
        C(2,4)=-4
        C(3,1)=2
        C(3,2)=2
        C(3,3)=5
        C(3,4)=1
        NUM=3
        CALL SOLMATR(NUM,C,S)
        WRITE (*,*) S(1),S(2),S(3)
        STOP
        END

        SUBROUTIN SOLMATR (NUM,C,S)
         DIMENSION C(4,4),S(3)
         Z=0
         DO 200 I=1, NUM
            IF(C(I,I).NE.Z) GOTO 220
               DO 230 J=I+1, NUM
                  IF (C(C,J.NE.Z) GOTO 250

230    CONTINUE
250    CALL PIVOT (I,J,NUM,C)
220    DIV=C(I,I)
       DO 270 J=1, NUM
          C(I,J)=C(I,J)DIV
270    CONTINU
       DO 290 II=1,NUM
          IF(II.EQ.I) GOTO 290
             RR=C(II,I)
             DO 310 K=1,NUM+1
                C(II,K)-RR*C(I,K)
310            CONTINUE
290   CONTINUE
200   CONTINUE
      DO 330 I=1, NUM
         S(I)=C(I,NUM+1)
330   CONTINUE
      RETURN
     END

     SUBROUTINE PIVOT(I,J,NUM,C)
       DIMENSION C(4,4)
       DO 400 K=1, NUM+1
          TRAN=C(I,K)
          C(I,K)=C(J,K)
          C(J,K)=TRAN

400   CONTINUE
      RETURN
    END



這是本週作業

我當然不會做

但是我會GOOGLE

我找到的範例是用Fortram寫的

我改寫之後當然還是Fortran(理論上....)

可是學校電腦沒裝Fortran

誰能救我一下

幫我改成C++ ><

--
愛情就像棒球
一支三殺打
會讓前面兩支安打的成果
煙消雲散
要得分只有靠連續安打
或是全壘打

--
※ 發信站: 批踢踢兔(ptt2.cc) 
◆ From: 140.137.9.123


--
愛情就像棒球
一支三殺打
會讓前面兩支安打的成果
煙消雲散
要得分只有靠連續安打
或是全壘打

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.137.9.123
推 skychenkid:很急嗎? 我這兩個禮拜之內就不了你=   ="
→ skychenkid:自己都快掛了... 要寫程式又要燒版子... 作業一堆阿 囧
→ f94:多謝sky
推 vince920741:我們這學期也是在燒板子~ 每週一Lab
推 kshs920742:FPGA???  再寫verilog???
推 skychenkid:樓上正解... 救命阿Q____Q
推 kshs920742:我們大一下有修過一點點  還要看過書才知道orz
→ kshs920742:如果是我   我會考慮用克拉瑪解耶
→ kshs920742:壤使用者輸入方程式  限制他ax+by+cz=d
→ kshs920742:剩下的用行列式去解  多點參數跑就是了XD