| 5 | 1/1 | 返回列表 |
| 查看: 1904 | 回復(fù): 12 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
霹靂旋風(fēng)榮譽(yù)版主 (著名寫手)
|
[求助]
求高人幫忙看一下我這個(gè)Fortran程序的錯(cuò)誤出在哪里?
|
|
|
程序代碼如下: 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 程序編譯通過(guò),但是出現(xiàn)計(jì)算數(shù)據(jù)無(wú)法輸出的問(wèn)題,請(qǐng)高手幫忙看一下。不勝感激。 |
至尊木蟲 (職業(yè)作家)
金蟲 (著名寫手)
木蟲 (正式寫手)
榮譽(yù)版主 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +5 | 生物工程調(diào)劑 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +3 | 石瑞0426 2026-03-19 | 3/150 |
|
|
[考研] 317求調(diào)劑 +9 | 申子申申 2026-03-19 | 15/750 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +6 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無(wú)極 2026-03-20 | 5/250 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|