| 5 | 1/1 | 返回列表 |
| 查看: 1921 | 回復(fù): 12 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
霹靂旋風(fēng)榮譽版主 (著名寫手)
|
[求助]
求高人幫忙看一下我這個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ù)無法輸出的問題,請高手幫忙看一下。不勝感激。 |
榮譽版主 (著名寫手)
金蟲 (著名寫手)
木蟲 (正式寫手)
榮譽版主 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 江蘇科技大學(xué)招材料研究生 +4 | Su032713. 2026-04-01 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +14 | 一樣YWY 2026-04-01 | 14/700 |
|
|
[考研] 一志愿鄭大材料工程290求調(diào)劑 +18 | Youth_ 2026-03-30 | 18/900 |
|
|
[考研] 324求調(diào)劑 +10 | hanamiko 2026-03-26 | 12/600 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-04-01 | 5/250 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +7 | 月月鳥木 2026-04-01 | 7/350 |
|
|
[考研] 085410人工智能 初試316分 求調(diào)劑 +3 | 殘星拂曙 2026-03-31 | 3/150 |
|
|
[考研] 一志愿北京科技大學(xué)材料學(xué)碩328分求調(diào)劑 +4 | 1段時間 2026-03-31 | 5/250 |
|
|
[考研] 299求調(diào)劑 +8 | 嗯嗯嗯嗯2 2026-03-27 | 8/400 |
|
|
[考研] 340求調(diào)劑 +4 | 希望如此i 2026-03-31 | 4/200 |
|
|
[考研] 一志愿浙江大學(xué)工科動力工程370,數(shù)一121,專業(yè)課135,現(xiàn)在能去哪里 +3 | 080700調(diào)劑 2026-03-30 | 4/200 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-29 | 3/150 |
|
|
[考研] 085601一志愿西北工業(yè)大學(xué)初試346 +4 | 085601初試346 2026-03-30 | 4/200 |
|
|
[考研] 285求調(diào)劑 +6 | AZMK 2026-03-29 | 9/450 |
|
|
[考研] 0703化學(xué)321分求調(diào)劑 +10 | 三dd. 2026-03-30 | 11/550 |
|
|
[考研]
|
gr哈哈哈 2026-03-28 | 3/150 |
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|