| 查看: 1451 | 回复: 8 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】急需一个程序
|
|||
| 急需：一个求解一元多次方程的Fortran程序，谢谢！ |
回复此楼» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有25人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有70人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 需要用matlab拟合数据的外行人急需帮助
已经有10人回复
- MATLAB 编程求解普通实验数据时间序列的最大李雅普诺夫指数！！！
已经有12人回复
- 两本经典的入门教程C Primer plus（第五版）和C++ Primer plus（第五版）
已经有58人回复
- 求助~~!!请懂SQL SERVER的进.
已经有13人回复
- 【求助】急需最新全国行政区划数据和常见汉字及其拼音（也可以是声母的第一个字母）数
已经有4人回复
- 【求助】用matlab编程，需要将数据自检分类（高手进来指点下）
已经有5人回复
- 【求助】购买程序降温的烘箱求助！！
已经有11人回复
- 【讨论】哪位有VB6.0英文版啊，急需【已完成】
已经有3人回复
- 【分享】Origin 8.x 视频教程 （75个）
已经有177人回复
» 抢金币啦！回帖就可以得到:
-
祝福！基金中中中！！！！
+1/397
-
期待你的到来
+1/375
-
西安交大马伟/王亚洲课题组诚聘联培博士和硕士学生多名-有机生物电子学方向
+1/184
-
上海同城
+1/118
-
2026年生物学、大数据与信息工程国际会议 (BSBIE 2026)
+1/95
-
中国科学院空天院精密干涉测量团队 精密机械结构工程师招聘启事
+1/82
-
中国科学院空天院精密干涉测量团队招聘光学/算法/机械方向博士后/特别研究助理
+1/80
-
祝福！祈福！！基金中中中！
+5/55
-
加拿大/英属哥伦比亚大学曹彦凯课题组招收博后 [机器人/控制方向]
+1/36
-
加拿大/英属哥伦比亚大学曹彦凯课题组招收全奖博士/博后 [机器人/控制方向]
+1/34
-
华南理工大学招2027年生物材料、药物递送申请考核制博士生
+1/29
-
祈福！
+1/25
-
希望这篇论文能中
+2/12
-
华南理工，博士后
+1/12
-
南京医科大学-冯璋课题组-招收2027级生化与分子生物学专业博士和硕士研究生
+1/4
-
UTS 招收2027年入学博士生
+1/3
-
中国科学技术大学环境学院招聘招生
+1/2
-
海南大学化学院—功能分子器件团队2027博士招生+博士后招聘
+1/2
-
中国科学技术大学苏州高等研究院张振教授2026级硕士生招生（推免）
+1/1
-
KFUPM Turki Baroud组招全奖博士/博后｜下一代膜材料・水凝胶・清洁水
+1/1
8楼2010-12-17 20:41:58
2楼2010-11-30 10:05:25
★
余泽成(金币+1):谢谢参与应助，欢迎常来程序语言版！ 2010-12-24 16:15:28
余泽成(金币+1):谢谢参与应助，欢迎常来程序语言版！ 2010-12-24 16:15:28
从书上给你找的！ DIMENSION A(4,4),B(4),X(4),JS(4) DOUBLE PRECISION A,B,X DATA A/0.2368,0.1968,0.1582,1.1161,0.2471,0.2071,1.1675,0.1254 * ,0.2568,1.2168,0.1768,0.1397,1.2671,0.2271,0.1871,0.1490/ DATA B/1.8471,1.7471,1.6471,1.5471/ N=4 CALL AGAUS(A,B,N,X,L,JS) IF (L.NE.0) THEN WRITE(*,10) (I,X(I),I=1,4) END IF 10 FORMAT(1X,'X(',I2,' )=',D15.6) END SUBROUTINE AGAUS(A,B,N,X,L,JS) DIMENSION A(N,N),X(N),B(N),JS(N) DOUBLE PRECISION A,B,X,T L=1 DO 50 K=1,N-1 D=0.0 DO 210 I=K,N DO 210 J=K,N IF (ABS(A(I,J)).GT.D) THEN D=ABS(A(I,J)) JS(K)=J IS=I END IF 210 CONTINUE IF (D+1.0.EQ.1.0) THEN L=0 ELSE IF (JS(K).NE.K) THEN DO 220 I=1,N T=A(I,K) A(I,K)=A(I,JS(K)) A(I,JS(K))=T 220 CONTINUE END IF IF (IS.NE.K) THEN DO 230 J=K,N T=A(K,J) A(K,J)=A(IS,J) A(IS,J)=T 230 CONTINUE T=B(K) B(K)=B(IS) B(IS)=T END IF END IF IF (L.EQ.0) THEN WRITE(*,100) RETURN END IF DO 10 J=K+1,N A(K,J)=A(K,J)/A(K,K) 10 CONTINUE B(K)=B(K)/A(K,K) DO 30 I=K+1,N DO 20 J=K+1,N A(I,J)=A(I,J)-A(I,K)*A(K,J) 20 CONTINUE B(I)=B(I)-A(I,K)*B(K) 30 CONTINUE 50 CONTINUE IF (ABS(A(N,N))+1.0.EQ.1.0) THEN L=0 WRITE(*,100) RETURN END IF X(N)=B(N)/A(N,N) DO 70 I=N-1,1,-1 T=0.0 DO 60 J=I+1,N T=T+A(I,J)*X(J) 60 CONTINUE X(I)=B(I)-T 70 CONTINUE 100 FORMAT(1X,' FAIL ') JS(N)=N DO 150 K=N,1,-1 IF (JS(K).NE.K) THEN T=X(K) X(K)=X(JS(K)) X(JS(K))=T END IF 150 CONTINUE RETURN END |
3楼2010-11-30 16:41:18
4楼2010-11-30 21:52:48
