作者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