| 5 | 1/1 | 返回列表 |
| 查看: 1920 | 回復(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 程序編譯通過,但是出現(xiàn)計(jì)算數(shù)據(jù)無法輸出的問題,請(qǐng)高手幫忙看一下。不勝感激。 |
至尊木蟲 (知名作家)
金蟲 (著名寫手)
木蟲 (正式寫手)
榮譽(yù)版主 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 332求調(diào)劑 +10 | 032500 2026-03-25 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 材料專碩 085600求調(diào)劑 +3 | BBQ233 2026-03-30 | 3/150 |
|
|
[考研] 297求調(diào)劑 +13 | 田洪有 2026-03-26 | 13/650 |
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +25 | RichLi_ 2026-03-25 | 25/1250 |
|
|
[考研] 318求調(diào)劑 +7 | 陳晨79 2026-03-30 | 7/350 |
|
|
[考研] 求調(diào)劑 +5 | 青春裁為三截 2026-03-29 | 5/250 |
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +6 | 深V宿舍吧 2026-03-29 | 6/300 |
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
10+4
|
chibby 2026-03-25 | 4/200 |
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩,求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分,英語已過四六級(jí),六級(jí)521,求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 324求調(diào)劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
[考研] 求調(diào)劑 +6 | 研研,接電話 2026-03-24 | 7/350 |
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|