霹靂旋風

應助: 3 (幼兒園)
貴賓: 0.95
金幣: 6428.7
散金: 100
紅花: 4
帖子: 1124
在線: 175.7小時
蟲號: 62568
注冊: 2005-03-17
性別: GG
專業(yè): 爆炸與沖擊動力學

[求助] 求高人幫忙看一下我這個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

  程序編譯通過，但是出現計算數據無法輸出的問題，請高手幫忙看一下。不勝感激。

回復此樓

» 猜你喜歡

26考研一志愿中國石油大學(華東)305分求調劑已經有4人回復
286求調劑已經有5人回復
材料與化工（0856）304求 B區(qū) 調劑已經有3人回復
268求調劑已經有8人回復
265求調劑已經有14人回復
350 本科985求調劑，求老登收留已經有3人回復
【考研調劑】化學專業(yè) 281分，一志愿四川大學，誠心求調劑已經有9人回復
初試 317 已經有3人回復
279分求調劑一志愿211 已經有15人回復
279求調劑已經有4人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

FORTRAN屬性聲明出錯已經有8人回復
有會看fortran程序的，能幫忙看看這個程序嗎？已經有10人回復
求助計算機高人，編寫個Fortran程序已經有13人回復
請教一個fortran小程序編譯出錯的問題，謝謝已經有9人回復
寫了一個fortran90的小程序，編譯通不過，請大俠幫忙已經有59人回復
【求助】讀文本數據出現段錯誤[fortran] 已經有14人回復
【求助】如何寫FORTRAN程序實現求平均最近鄰距離已經有34人回復
【求助】運行fortran程序出現錯誤 : error LNK2019 怎么回事啊已經有8人回復
【求助】用fortran怎么描述下面這樣的情況，我是新手，請大家?guī)兔?/a> 已經有8人回復

【求助】fortran程序未出錯，但運行不出結果，中途中斷！已經有6人回復
【求助】此FORTRAN90程序沒有錯誤，但是運行不出來。已經有4人回復
【求助】用fortran求解大型線性方程組時出現的錯誤【已解決】已經有11人回復
【求助】求助個Fortran計算，應該不難，半天搞不定，請高手幫忙已經有11人回復

1樓 2012-02-18 21:24:39

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖置頂 ( 共有2個 )

霹靂旋風

榮譽版主 (著名寫手)

應助: 3 (幼兒園)
貴賓: 0.95
金幣: 6428.7
散金: 100
紅花: 4
帖子: 1124
在線: 175.7小時
蟲號: 62568
注冊: 2005-03-17
性別: GG
專業(yè): 爆炸與沖擊動力學

引用回帖:

樓: Originally posted by ygl820 at 2012-02-21 10:31:22:
運行后，test.dat生成了沒有。

運行后，生成了test.dat

贊一下

回復此樓

8樓2012-02-22 13:05:55

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

霹靂旋風

榮譽版主 (著名寫手)

應助: 3 (幼兒園)
貴賓: 0.95
金幣: 6428.7
散金: 100
紅花: 4
帖子: 1124
在線: 175.7小時
蟲號: 62568
注冊: 2005-03-17
性別: GG
專業(yè): 爆炸與沖擊動力學

引用回帖:

樓: Originally posted by snoopyzhao at 2012-02-23 13:54:41:
這個程序用 gfortran 沒有辦法編譯，顯示 cotan 未定義……

你這個 test.dat 是空的很正常啊，說明 L .LE. 0 啊，呵呵……

謝謝樓上，
程序中的cotan（）的含義是求角度的余切值呀。請高人指點哦。

贊一下

回復此樓

10樓2012-02-23 19:37:37

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

夕陽西下

金蟲 (著名寫手)

應助: 138 (高中生)
貴賓: 0.041
金幣: 9247.1
散金: 73
紅花: 28
沙發(fā): 1
帖子: 1268
在線: 1164.4小時
蟲號: 220750
注冊: 2006-03-20
專業(yè): 凝聚態(tài)物性 II ：電子結構
管轄: 物理

【答案】應助回帖

感謝參與，應助指數 +1
jjdg(金幣+1): 感謝應助 2012-02-20 01:39:45

open (2000, file='test.dat', status='old' )
write(2000,*) X(I)
write 2000這樣寫是不是有問題呢？

贊一下

回復此樓

2樓2012-02-19 09:17:19

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

haizhfly

木蟲 (正式寫手)

應助: 16 (小學生)
金幣: 3134.8
帖子: 316
在線: 67.4小時
蟲號: 925614
注冊: 2009-12-11
專業(yè): 天氣學

【答案】應助回帖

感謝參與，應助指數 +1
jjdg(金幣+1): 感謝應助 2012-02-20 01:39:55

open (2000, file='test.dat', status='old' )
write(2000,*) X(I)

status='old' 要求輸出文本 test.dat是已經存在的文件，刪去這一文件狀態(tài)控制語句試一下看，個人覺得最大問題在這里

贊一下

回復此樓

3樓2012-02-19 11:21:31

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

霹靂旋風

榮譽版主 (著名寫手)

應助: 3 (幼兒園)
貴賓: 0.95
金幣: 6428.7
散金: 100
紅花: 4
帖子: 1124
在線: 175.7小時
蟲號: 62568
注冊: 2005-03-17
性別: GG
專業(yè): 爆炸與沖擊動力學

引用回帖:

樓: Originally posted by haizhfly at 2012-02-19 11:21:31:
open (2000, file='test.dat', status='old' )
write(2000,*) X(I)

status='old' 要求輸出文本 test.dat是已經存在的文件，刪去這一文件狀態(tài)控制語句試一下看，個人覺得最大問題在這里

這個我已經調試過。文件能編譯完，也能運行，就是沒有輸出結果。

贊一下

回復此樓

4樓2012-02-19 19:48:38

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

ygl820

至尊木蟲 (知名作家)

應助: 133 (高中生)
金幣: 52898.5
散金: 1007
紅花: 46
沙發(fā): 94
帖子: 7094
在線: 1047小時
蟲號: 81941
注冊: 2005-07-22
性別: GG
專業(yè): 有機合成

【答案】應助回帖

感謝參與，應助指數 +1
霹靂旋風(金幣+10): 謝謝 2012-02-21 10:29:01

按你的意思，將文件test.dat作為輸出了，運行完成后后，打開它。fortran文件操作性不強，看看test.dat是否只讀。

贊一下

回復此樓

5樓2012-02-21 10:24:17

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

霹靂旋風

榮譽版主 (著名寫手)

應助: 3 (幼兒園)
貴賓: 0.95
金幣: 6428.7
散金: 100
紅花: 4
帖子: 1124
在線: 175.7小時
蟲號: 62568
注冊: 2005-03-17
性別: GG
專業(yè): 爆炸與沖擊動力學

引用回帖:

樓: Originally posted by ygl820 at 2012-02-21 10:24:17:
按你的意思，將文件test.dat作為輸出了，運行完成后后，打開它。fortran文件操作性不強，看看test.dat是否只讀。

也不是只讀。運行可以就是沒有輸出數據。急死我了。

贊一下

回復此樓

6樓2012-02-21 10:28:44

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

ygl820

至尊木蟲 (知名作家)

應助: 133 (高中生)
金幣: 52898.5
散金: 1007
紅花: 46
沙發(fā): 94
帖子: 7094
在線: 1047小時
蟲號: 81941
注冊: 2005-07-22
性別: GG
專業(yè): 有機合成

運行后，test.dat生成了沒有。

贊一下

回復此樓

7樓2012-02-21 10:31:22

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學

余澤成(金幣+1): 謝謝參與應助！ 2012-02-25 04:29:52

這個程序用 gfortran 沒有辦法編譯，顯示 cotan 未定義……

你這個 test.dat 是空的很正常啊，說明 L .LE. 0 啊，呵呵……

贊一下

回復此樓

9樓2012-02-23 13:54:41

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 265求調劑 +12	梁梁校校 2026-03-19	14/700	2026-03-21 13:38 by lature00
[考研] 材料學學碩080502 337求調劑-一志愿華中科技大學 +4	順順順mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 材料 271求調劑 +3	展信悅_ 2026-03-21	3/150	2026-03-21 10:17 by 王笑宇一定是研?/a>
[考研] 265求調劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 328求調劑，英語六級551，有科研經歷 +4	生物工程調劑 2026-03-17	8/400	2026-03-21 02:12 by JourneyLucky
[考研] 271材料工程求調劑 +8	.6lL 2026-03-18	8/400	2026-03-21 00:58 by JourneyLucky
[考研] 考研調劑求學校推薦 +3	伯樂29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 一志愿蘇州大學材料求調劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 中南大學化學學碩337求調劑 +3	niko- 2026-03-19	6/300	2026-03-20 21:58 by luoyongfeng
[考研] 材料學碩297已過四六級求調劑推薦 +11	adaie 2026-03-19	11/550	2026-03-20 21:30 by laoshidan
[考研] 0817 化學工程 299分求調劑有科研經歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 求調劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 081700化工學碩調劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 梁成偉老師課題組歡迎你的加入 +9	一鴨鴨喲 2026-03-14	11/550	2026-03-19 17:22 by ！本暗一次！
[考研] 312求調劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 293求調劑 +11	zjl的號 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 東南大學364求調劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 277材料科學與工程080500求調劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運氣yunqi
[考研] 復試調劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen