 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 876 | 回復(fù): 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ù)運(yùn)行的話(huà)該怎么修改此程序啊？？？？？請(qǐng)賜教！ |
回復(fù)此樓» 猜你喜歡
323求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑
已經(jīng)有4人回復(fù)
-
352求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
已經(jīng)有8人回復(fù)
-
277材料科學(xué)與工程080500求調(diào)劑
已經(jīng)有7人回復(fù)
-
317求調(diào)劑
已經(jīng)有18人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
-
280分求調(diào)劑 一志愿085802
已經(jīng)有7人回復(fù)
-
0854電子信息求調(diào)劑
已經(jīng)有3人回復(fù)
-
263求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 寫(xiě)了一個(gè)fortran90的小程序，編譯通不過(guò)，請(qǐng)大俠幫忙
已經(jīng)有59人回復(fù)
- 【求助】求高人幫我修改程序
已經(jīng)有5人回復(fù)
- 【求助】幫忙改寫(xiě)展寬程序
已經(jīng)有11人回復(fù)
- 【原創(chuàng)】vasp的POSCAR轉(zhuǎn)換成cif格式程序修改版
已經(jīng)有7人回復(fù)
- 【求助】修改MATLAB支持向量機(jī)核函數(shù)程序遇到的疑問(wèn)？求助
已經(jīng)有5人回復(fù)
snoopyzhao
至尊木蟲(chóng) (職業(yè)作家)
- 程序強(qiáng)帖: 16
- 應(yīng)助: 157 (高中生)
- 貴賓: 0.02
- 金幣: 18844.7
- 紅花: 29
- 帖子: 3803
- 在線(xiàn): 1422.4小時(shí)
- 蟲(chóng)號(hào): 183750
- 注冊(cè): 2006-02-13
- 專(zhuān)業(yè): 污染生態(tài)化學(xué)
2樓2010-04-26 14:36:37
holmescn
金蟲(chóng) (正式寫(xiě)手)
- 程序強(qiáng)帖: 37
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1918.8
- 散金: 275
- 紅花: 1
- 帖子: 699
- 在線(xiàn): 102.6小時(shí)
- 蟲(chóng)號(hào): 913482
- 注冊(cè): 2009-11-26
- 性別: GG
- 專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2010-04-26 17:06:38
4樓2010-04-27 13:39:01
5樓2010-04-27 13:41:28
snoopyzhao
至尊木蟲(chóng) (職業(yè)作家)
- 程序強(qiáng)帖: 16
- 應(yīng)助: 157 (高中生)
- 貴賓: 0.02
- 金幣: 18844.7
- 紅花: 29
- 帖子: 3803
- 在線(xiàn): 1422.4小時(shí)
- 蟲(chóng)號(hào): 183750
- 注冊(cè): 2006-02-13
- 專(zhuān)業(yè): 污染生態(tài)化學(xué)
6樓2010-04-27 15:29:57
7樓2010-04-28 12:22:18
snoopyzhao
至尊木蟲(chóng) (職業(yè)作家)
- 程序強(qiáng)帖: 16
- 應(yīng)助: 157 (高中生)
- 貴賓: 0.02
- 金幣: 18844.7
- 紅花: 29
- 帖子: 3803
- 在線(xiàn): 1422.4小時(shí)
- 蟲(chóng)號(hào): 183750
- 注冊(cè): 2006-02-13
- 專(zhuān)業(yè): 污染生態(tài)化學(xué)
8樓2010-04-28 15:28:26
9樓2010-04-28 22:07:10
snoopyzhao
至尊木蟲(chóng) (職業(yè)作家)
- 程序強(qiáng)帖: 16
- 應(yīng)助: 157 (高中生)
- 貴賓: 0.02
- 金幣: 18844.7
- 紅花: 29
- 帖子: 3803
- 在線(xiàn): 1422.4小時(shí)
- 蟲(chóng)號(hào): 183750
- 注冊(cè): 2006-02-13
- 專(zhuān)業(yè): 污染生態(tài)化學(xué)
10樓2010-04-29 11:24:01
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 215215 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 352求調(diào)劑 +3 | 大米飯！ 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 0854電子信息求調(diào)劑 +3 | α____ 2026-03-22 | 3/150 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +3 | niko- 2026-03-22 | 3/150 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 4/200 |
|
|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚(yú)ly 2026-03-21 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語(yǔ) 2026-03-20 | 3/150 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 材料專(zhuān)業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿武漢理工材料工程專(zhuān)碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù)
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3 | 小木蟲(chóng)200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|