应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】大家帮忙看下错在哪！一直提示一个错误【已完成】
51/1返回列表
查看: 803  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

215215

银虫 (小有名气)

[交流] 【求助】大家帮忙看下错在哪！一直提示一个错误【已完成】 已有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 ]
回复此楼

» 猜你喜欢
nevergiveup!
1楼 2010-04-11 22:52:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjyl

金虫 (正式写手)
★ ★ ★ ★ ★ ★
小木虫(金币+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
赞一下(4人) 回复此楼
2楼2010-04-11 23:18:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjyl

金虫 (正式写手)
★ ★ ★
jjdg(金币+2):辛苦了 2010-04-11 23:25
wangen994(金币+1):活动期间额外奖励 2010-04-12 08:24
就改了

    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)


还有就是不要用中文文件名。
赞一下(2人) 回复此楼
3楼2010-04-11 23:19:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

215215

银虫 (小有名气)
★
jjdg(金币+1):欢迎新虫 2010-04-12 14:15
引用回帖:
Originally posted by tjyl at 2010-04-11 23:19:25:
就改了

    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)


还有就是不要用 ...

谢谢哈！
回复此楼
nevergiveup!
4楼2010-04-12 13:05:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

215215

银虫 (小有名气)

如果把b=0.9改为b属于区间【1，1.5】继续运行的话该怎么修改此程序？？谢谢赐教！！

引用回帖:
Originally posted by tjyl at 2010-04-11 23:18:42:
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= ...

如果把b=0.9修改为b属于区间【1，1.5】 继续运行此程序的话，该怎么修改此程序？？多谢赐教。
赞一下 回复此楼
nevergiveup!
5楼2010-04-21 12:06:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 215215 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 275求调剂 +7 1624447980 2026-04-08 8/400 2026-04-11 00:45 by dc19971224
[考研] 还有化工二轮调剂的学校吗 5+13 化工人999 2026-04-09 47/2350 2026-04-11 00:20 by 骑牛渡寒江
[考研] 287求调剂 +15 Fnhc 2026-04-07 21/1050 2026-04-10 19:09 by chemisry
[考研] 307求调剂 +8 tzq94092 2026-04-10 8/400 2026-04-10 17:33 by 286640313
[考研] 314求调剂 +14 weltZeng 2026-04-09 14/700 2026-04-09 23:14 by wolf97
[考研] 调剂 +8 只叙离别辞 2026-04-09 10/500 2026-04-09 20:25 by zl8213662
[考研] 085501机械英二77总分294求调剂，接受跨专业学习 +6 守法公民亓纪 2026-04-08 6/300 2026-04-09 15:55 by wp06
[考研] 349学科化学045106求调剂，化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 085801 总分275 本科新能源 求调剂 +8 bradoner 2026-04-08 9/450 2026-04-09 13:43 by only周
[考研] 280求调剂 +5 兮兮夜夜 2026-04-09 8/400 2026-04-09 11:15 by 兮兮夜夜
[考研] 308求调剂 +17 墨墨漠 2026-04-06 17/850 2026-04-09 09:25 by 壹往無前
[考研] 求调剂 +9 月@163.com 2026-04-07 11/550 2026-04-08 14:48 by qlm5820
[考研] 283求调剂 +19 A child 2026-04-04 19/950 2026-04-08 14:26 by xingguangj
[考研] 312求调剂 +4 Say Never 2026-04-04 4/200 2026-04-08 08:41 by barlinike
[考研] 0854电子信息319求调剂（接受跨专业调剂） +5 星星不眨眼喽 2026-04-05 6/300 2026-04-07 22:16 by hemengdong
[考研] 一志愿南科大生物学297分，求调剂推荐 +8 Y-yyusx 2026-04-06 9/450 2026-04-07 19:38 by biomichael
[考研] 285求调剂 +5 mapmath 2026-04-06 6/300 2026-04-06 17:18 by 蓝云思雨
[考研] 一志愿C9的化学工程（085602） 340分，感觉校内调剂无望，求调剂 +12 万事宜臻 2026-04-04 12/600 2026-04-06 07:46 by 无际的草原
[考研] 材料调剂 +14 壹贰贰亿 2026-04-04 14/700 2026-04-05 23:31 by 来看流星雨10
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-04-04 5/250 2026-04-05 14:06 by imissbao
信息提示
关闭
请填处理意见
关闭