| 查看: 1936 | 回复: 12 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
霹雳旋风荣誉版主 (著名写手)
|
[求助]
求高人帮忙看一下我这个Fortran程序的错误出在哪里？
|
|
|
程序代码如下： program ganeiliwucha implicit none External FS Integer I,L,N Dimension X(2),Y(2),C(2) Double Precision X, Y, EPS,C open (2000, file='test.dat', status='old' ) Data X/0.00,0.01/ EPS=1.0E-06 CALL DSNSE(2,EPS, X, Y, FS, L) IF(L.GT.0)then C(1)=cos(X(2))/sin(X(1)+X(2))*tan(20.0/180*3.14)-1.0 ! AC杆的内力误差 C(2)=cos(X(1))/sin(X(1)+X(2))*sin(20.0/180*3.14)-1.0 ! BC杆的内力误差 do I=1,2 write(2000,*) X(I) end do end IF close(2000) end SUBROUTINE DSNSE(N,EPS,X,Y,FS,L) DIMENSION X(N),Y(N) DOUBLE PRECISION X,Y,F,D,S L=500 5 CALL FS(X,N,F,Y) IF (F.GE.EPS) THEN L=L-1 IF (L.EQ.0) RETURN D=0.0 DO 20 J=1,N 20 D=D+Y(J)*Y(J) IF (D+1.0.EQ.1.0) THEN L=-1 RETURN END IF S=F/D DO 30 I=1,N 30 X(I)=X(I)-S*Y(I) GOTO 5 END IF RETURN END SUBROUTINE FS(X, N, F, Y) DIMENSION X(N), Y(N) Double precision X,Y,F,F1,F2,DF1,DF2,A,B A=20.0/180*3.14 K=0.5 F1=X(1)+X(2)-cos(A) F2=sin(X(1))*cos(A)*(1+K*cotan(A))+SIN(X(2))*(1.0-K/sin(A))-cos(A) F=F1*F1+F2*F2 DF1=1.0 DF2=cos(X(1))*cos(A)*(1+K*cotan(A)) Y(1)=2.0*(F1*DF1+F2*DF2) DF1=1.0 DF2=COS(X(2))*(1.0-K/sin(A)) Y(2)=2.0*(F1*DF1+F2*DF2) RETURN END 程序编译通过，但是出现计算数据无法输出的问题，请高手帮忙看一下。不胜感激。 |
回复此楼» 猜你喜欢
国自然面上和省基金B类撒花
已经有4人回复
-
华师大读博
已经有3人回复
-
急需审稿人！！！
已经有3人回复
-
申博/考博
已经有8人回复
-
一志愿A区211，22408 321求调剂
已经有8人回复
-
295分求调剂
已经有6人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有13人回复
-
085600材料与化工调剂
已经有6人回复
-
期刊推荐
已经有5人回复
-
有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- FORTRAN属性声明出错
已经有8人回复
- 有会看fortran程序的，能帮忙看看这个程序吗？
已经有10人回复
- 求助计算机高人，编写个Fortran程序
已经有13人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】读文本数据出现段错误[fortran]
已经有14人回复
- 【求助】如何写FORTRAN程序实现求平均最近邻距离
已经有34人回复
- 【求助】运行fortran程序出现错误 : error LNK2019 怎么回事啊
已经有8人回复
- 【求助】用fortran怎么描述下面这样的情况，我是新手，请大家帮忙
已经有8人回复
- 【求助】fortran程序未出错，但运行不出结果，中途中断！
已经有6人回复
- 【求助】此FORTRAN90程序没有错误，但是运行不出来。
已经有4人回复
- 【求助】用fortran求解大型线性方程组时出现的错误【已解决】
已经有11人回复
- 【求助】求助个Fortran计算，应该不难，半天搞不定，请高手帮忙
已经有11人回复
ygl820
至尊木虫 (知名作家)
- 应助: 133 (高中生)
- 金币: 54583.1
- 散金: 1007
- 红花: 48
- 沙发: 94
- 帖子: 7384
- 在线: 1081.2小时
- 虫号: 81941
- 注册: 2005-07-22
- 性别: GG
- 专业: 有机合成
7楼2012-02-21 10:31:22
夕阳西下
金虫 (著名写手)
- 应助: 138 (高中生)
- 贵宾: 0.041
- 金币: 9247.1
- 散金: 73
- 红花: 28
- 沙发: 1
- 帖子: 1268
- 在线: 1164.4小时
- 虫号: 220750
- 注册: 2006-03-20
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 物理
2楼2012-02-19 09:17:19
3楼2012-02-19 11:21:31
霹雳旋风
荣誉版主 (著名写手)
- 应助: 3 (幼儿园)
- 贵宾: 0.95
- 金币: 6428.7
- 散金: 100
- 红花: 4
- 帖子: 1124
- 在线: 175.7小时
- 虫号: 62568
- 注册: 2005-03-17
- 性别: GG
- 专业: 爆炸与冲击动力学
4楼2012-02-19 19:48:38
5