版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

霹靂旋風

榮譽版主 (著名寫手)

應助: 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桿的內(nèi)力誤差
C(2)=cos(X(1))/sin(X(1)+X(2))*sin(20.0/180*3.14)-1.0 ! BC桿的內(nèi)力誤差
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

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

回復此樓

» 猜你喜歡

【考研調(diào)劑】化學專業(yè) 281分，一志愿四川大學，誠心求調(diào)劑已經(jīng)有11人回復
一志愿深大，0703化學，總分302，求調(diào)劑已經(jīng)有4人回復
材料學碩333求調(diào)劑已經(jīng)有3人回復
材料求調(diào)劑已經(jīng)有3人回復
333求調(diào)劑已經(jīng)有5人回復
求助已經(jīng)有6人回復
278求調(diào)劑已經(jīng)有9人回復
297求調(diào)劑已經(jīng)有12人回復
0703化學297求調(diào)劑已經(jīng)有3人回復
材料學碩301分求調(diào)劑已經(jīng)有6人回復

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

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

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

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

已閱回復此樓關(guān)注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

已閱回復此樓關(guān)注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

已閱回復此樓關(guān)注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 ：電子結(jié)構(gòu)
管轄: 物理

【答案】應助回帖

感謝參與，應助指數(shù) +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

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

haizhfly

木蟲 (正式寫手)

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

【答案】應助回帖

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

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

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

贊一下

回復此樓

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

已閱回復此樓關(guān)注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是已經(jīng)存在的文件，刪去這一文件狀態(tài)控制語句試一下看，個人覺得最大問題在這里

這個我已經(jīng)調(diào)試過。文件能編譯完，也能運行，就是沒有輸出結(jié)果。

贊一下

回復此樓

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

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

ygl820

至尊木蟲 (知名作家)

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

【答案】應助回帖

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

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

贊一下

回復此樓

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

已閱回復此樓關(guān)注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是否只讀。

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

贊一下

回復此樓

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

已閱回復此樓關(guān)注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

已閱回復此樓關(guān)注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

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主霹靂旋風的主題更新

返回列表

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

[考研] 一志愿深大，0703化學，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學員8dgXkO
[考研] 313求調(diào)劑 +4	肆叁貳壹22 2026-03-19	4/200	2026-03-21 17:33 by ColorlessPI
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	14/700	2026-03-21 13:38 by lature00
[考研] 求調(diào)劑 +6	Mqqqqqq 2026-03-19	6/300	2026-03-21 08:04 by JourneyLucky
[考研] 316求調(diào)劑 +6	梁茜雯 2026-03-19	6/300	2026-03-21 06:32 by Ecowxq666！
[考研] 070300化學319求調(diào)劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 一志愿西北大學，070300化學學碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-18	3/150	2026-03-21 00:46 by JourneyLucky
[考研] 南京大學化學376求調(diào)劑 +3	hisfailed 2026-03-19	6/300	2026-03-20 23:43 by hisfailed
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 287求調(diào)劑 +7	晨昏線與星海 2026-03-19	8/400	2026-03-20 22:19 by JourneyLucky
[考研] 一志愿西北大學，070300化學學碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 材料與化工 322求調(diào)劑 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] 0817 化學工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 一志愿南京航空航天大學大學，080500材料科學與工程學碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 0703化學調(diào)劑，六級已過，有科研經(jīng)歷 +13	曦熙兮 2026-03-15	13/650	2026-03-20 19:35 by Dream007008
[考研] 環(huán)境工程調(diào)劑 +9	大可digkids 2026-03-16	9/450	2026-03-20 17:38 by 醉在風里
[考研] 材料與化工求調(diào)劑 +7	為學666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 0703化學336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 334求調(diào)劑 +3	志存高遠意在機?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 326求調(diào)劑 +4	諾貝爾化學獎覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學獎覬?/a>

24小時熱門版塊排行榜

[求助] 求高人幫忙看一下我這個Fortran程序的錯誤出在哪里？

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【答案】應助回帖

【答案】應助回帖

【答案】應助回帖

[求助] 求高人幫忙看一下我這個Fortran程序的錯誤出在哪里？

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