| 查看: 741 | 回复: 4 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】大家帮忙看下错在哪！一直提示一个错误【已完成】已有2人参与
|
|||
|
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 --------------------Configuration: 单摆 - Win32 Debug-------------------- Compiling Fortran... E:\Program Files\Microsoft Visual Studio\MyProjects\aaa\单摆.f90 f90: Severe: No such file or directory ... file is 'E:\Program Files\Microsoft Visual Studio\MyProjects\aaa\μ￥°ú.f90' Error executing df.exe. 单摆.exe - 1 error(s), 0 warning(s) [ Last edited by wangen994 on 2010-4-12 at 15:47 ] |
回复此楼» 猜你喜欢
AI论文写作工具：是科研加速器还是学术作弊器？
已经有5人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
tjyl
金虫 (正式写手)
- 程序强帖: 2
- 应助: 0 (幼儿园)
- 金币: 3218.1
- 红花: 2
- 帖子: 576
- 在线: 156.6小时
- 虫号: 765184
- 注册: 2009-05-07
- 专业: 生物无机化学
- 管辖: 程序语言
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
jjdg(金币+1):感谢参与 2010-04-12 01:07
wangen994(金币+2):感谢参与讨论 2010-04-12 08:24
wangen994(金币+2):活动期间额外奖励 2010-04-12 08:24
小木虫(金币+0.5):给个红包，谢谢回帖交流
jjdg(金币+1):感谢参与 2010-04-12 01:07
wangen994(金币+2):感谢参与讨论 2010-04-12 08:24
wangen994(金币+2):活动期间额外奖励 2010-04-12 08:24
|
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 |
2楼2010-04-11 23:18:42
tjyl
金虫 (正式写手)
- 程序强帖: 2
- 应助: 0 (幼儿园)
- 金币: 3218.1
- 红花: 2
- 帖子: 576
- 在线: 156.6小时
- 虫号: 765184
- 注册: 2009-05-07
- 专业: 生物无机化学
- 管辖: 程序语言
3楼2010-04-11 23:19:25
4楼2010-04-12 13:05:45
5楼2010-04-21 12:06:26