 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術(shù)科研互動平臺 » 計算模擬區(qū) » 程序語言 » Fortran » 【求助】程序修改！
| 2 | 1/1 | 返回列表 |
| 查看: 883 | 回復: 10 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[交流]
【求助】程序修改！ 已有4人參與
|
|||
|
PROGRAM PENDULUM ! ! Main Program for a driven pendulum under damping solved with ! the fourth-order Runge-Kutta algorithm. ! Parameters: Q, B,and W (omega_0). ! PARAMETER (N=1000,L=100,M=1) DIMENSION Y(2,N) COMMON/CONST/Q,B,W ! PI = 4.0*ATAN(1.0) H = 3.0*PI/L Q = 0.5 B = 0.9 W = 2.0/3.0 Y(1,1) = 0.0 Y(2,1) = 2.0 ! ! Using the Runge-Kutta algorithm to integrate the ! equation ! DO 100 I = 1, N-1 T = H*I Y1 = Y(1,I) Y2 = Y(2,I) DK11 = H*G1(Y1,Y2,T) DK21 = H*G2(Y1,Y2,T) DK12 = H*G1((Y1+DK11/2.0),(Y2+DK21/2.0),(T+H/2.0)) DK22 = H*G2((Y1+DK11/2.0),(Y2+DK21/2.0),(T+H/2.0)) DK13 = H*G1((Y1+DK12/2.0),(Y2+DK22/2.0),(T+H/2.0)) DK23 = H*G2((Y1+DK12/2.0),(Y2+DK22/2.0),(T+H/2.0)) DK14 = H*G1((Y1+DK13),(Y2+DK23),(T+H)) DK24 = H*G2((Y1+DK13),(Y2+DK23),(T+H)) Y(1,I+1) = Y(1,I)+(DK11+2.0*(DK12+DK13)+DK14)/6.0 Y(2,I+1) = Y(2,I)+(DK21+2.0*(DK22+DK23)+DK24)/6.0 ! ! Bring theta back to the region [-pi,pi] ! IF(ABS(Y(1,I+1)).GT.PI) THEN Y(1,I+1) = Y(1,I+1)-2.0*PI*ABS(Y(1,I+1))/Y(1,I+1) END IF 100 CONTINUE ! WRITE(6,999)(Y(1,I),Y(2,I),I=1,N,M) STOP 999 FORMAT (2F16.8) END ! FUNCTION G1 (Y1,Y2,T) COMMON /CONST/Q,B,W G1=Y2 RETURN END ! FUNCTION G2 (Y1,Y2,T) COMMON /CONST/Q,B,W G2 = -Q*Y2-SIN(Y1)+B*COS(W*T) RETURN END ———————————————————————————————————— 以上程序如果把B=0.9修改為B在區(qū)間【1，1.5】，要想繼續(xù)運行的話該怎么修改此程序��？？？？？請賜教！ |
回復此樓» 猜你喜歡
333求調(diào)劑
已經(jīng)有10人回復
-
307求調(diào)劑
已經(jīng)有3人回復
-
考研化學308分求調(diào)劑
已經(jīng)有4人回復
-
298求調(diào)劑
已經(jīng)有7人回復
-
263求調(diào)劑
已經(jīng)有5人回復
-
276求調(diào)劑。有半年電池和半年高分子實習經(jīng)歷
已經(jīng)有6人回復
-
求調(diào)劑材料學碩080500，總分289分
已經(jīng)有21人回復
-
北科281學碩材料求調(diào)劑
已經(jīng)有8人回復
-
08工學調(diào)劑
已經(jīng)有9人回復
-
291 求調(diào)劑
已經(jīng)有4人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 寫了一個fortran90的小程序，編譯通不過，請大俠幫忙
已經(jīng)有59人回復
- 【求助】求高人幫我修改程序
已經(jīng)有5人回復
- 【求助】幫忙改寫展寬程序
已經(jīng)有11人回復
- 【原創(chuàng)】vasp的POSCAR轉(zhuǎn)換成cif格式程序修改版
已經(jīng)有7人回復
- 【求助】修改MATLAB支持向量機核函數(shù)程序遇到的疑問？求助
已經(jīng)有5人回復
hgqking
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 43.5
- 帖子: 40
- 在線: 11.7小時
- 蟲號: 991045
- 注冊: 2010-04-07
- 性別: GG
- 專業(yè): 計算機應(yīng)用技術(shù)
11樓2010-04-29 11:46:22
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 215215 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招08考數(shù)學 +4 | laoshidan 2026-03-20 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
|
[基金申請] 山東省面上項目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
|
[考研] 一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
|
[考研] 一志愿華南師大 070300（化學）304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿 西北大學 ，070300化學學碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學大學 ，080500材料科學與工程學碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿西安交通大學 學碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|