作者f94 (→罪ㄞˋ 仁甫ａ婷←)
標題[轉錄][大師] 程式設計苦惱
時間Sat May 3 13:17:16 2008
※ [本文轉錄自 f94 信箱]
作者:
[email protected] (
[email protected])
標題: [大師] 程式設計苦惱
時間: 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
1F：推 skychenkid:很急嗎? 我這兩個禮拜之內就不了你= =" 05/03 14:06
2F：→ skychenkid:自己都快掛了... 要寫程式又要燒版子... 作業一堆阿 囧 05/03 14:06
3F：→ f94:多謝sky 05/03 14:21
4F：推 vince920741:我們這學期也是在燒板子~ 每週一Lab 05/03 15:37
5F：推 kshs920742:FPGA??? 再寫verilog??? 05/03 23:44
6F：推 skychenkid:樓上正解... 救命阿Q____Q 05/04 09:55
7F：推 kshs920742:我們大一下有修過一點點 還要看過書才知道orz 05/04 21:25
8F：→ kshs920742:如果是我 我會考慮用克拉瑪解耶 05/04 21:26
9F：→ kshs920742:壤使用者輸入方程式 限制他ax+by+cz=d 05/04 21:26
10F：→ kshs920742:剩下的用行列式去解 多點參數跑就是了XD 05/04 21:26