| 查看: 475 | 回复: 2 | |||
| 【悬赏金币】回答本帖问题，作者我也要搞科研将赠送您 100 个金币 | |||
[求助]
用fortran求解牛顿第二运动定理的问题已有1人参与
|
|||
|
当位移x的初值x0=-12,v0=3时，进行运动学求解，采用ivprk求解微分方程，此时加速度为2，但我想要当得到的位移xt>16时，停止计算，将停止时的时间t2变成新的t0，停止时的速度vt和xt变成新的速度和位移的初值v0,x0,然后再进行第二次求微分方程，此时加速度为10。我的程序已经列出，但是总是报错，求大神指导。 program main INTEGER*4 MTIME, MXPARM,NP, IDO, MC, istep PARAMETER (MXPARM=50,MTIME=5000000000, NP=50) REAL*8 Y(NP),YPRIME(NP),Y2(NP),YPRIME2(NP) REAL*8 PARAM(MXPARM) real*8 x(-300:300),xe(0:2000,-300:300) real*8 vex(-300:300),ve(0:2000) real*8 t(0:2000),a(0:3000) real*8 t_stop, t_start,del_t REAL*8 TOL,T1,T2,t11,t21 external DIVPRK,DSET,FCN1,fcn2 t_start=0d0 t_stop=5d0 del_t=0.1d0 TOL=1.0D-3 ccccc初值:y(1)=x0=-12,y(2)=v0=3 y(1)=-12 y(2)=3 open(unit=3725,file= 1'C:\USERS\YANGBO\DESKTOP\test_tang\test_xiabishi 1\test123.dat', 1status='unknown') CALL DSET(MXPARM, 0.0D0, PARAM, 1) PARAM(10)=1.0D0 PARAM(4)=MTIME+1 istep=0 IDO=1 T1=T_START mt=0 istep=istep+1 T2=T1+DEL_T CALL DIVPRK(IDO,NP,FCN1,T1,T2,TOL,PARAM,Y) if (y(1)>16) then y2(1)=y(1) !新的初值,求出来的y(1)变成新的x初值 y2(2)=y(2) CALL DSET(MXPARM, 0.0D0, PARAM, 1) PARAM(10)=1.0D0 PARAM(4)=MTIME+1 istep=0 IDO=1 T11=t2 !y(1)>16时，t的初值变为停止时的t2 mt=0 end if 110 CONTINUE istep=istep+1 T21=T11+DEL_T CALL DIVPRK(IDO,NP,FCN2,T11,T21,TOL,PARAM,Y2) write(3725,1345)t2,y(2),y(1),y2(2),y2(1) 1345 format(1x,4(x,d13.6)) IF(T2<=T_STOP) GOTO 110 print *,'end' stop end subroutine FCN1 (NP, T2, Y, YPRIME) INTEGER*4 NP,npv21,npz,npx REAL*8 Y(NP),YPRIME(NP) a=2 YPRIME(1) = Y(2) YPRIME(2) = a return end subroutine subroutine FCN2 (NP, T21, Y2, YPRIME2) INTEGER*4 NP,npv21,npz,npx REAL*8 Y2(NP),YPRIME2(NP) a=10 YPRIME2(1) = Y2(2) YPRIME2(2) = a return end subroutine |
回复此楼» 猜你喜欢
关于要不要三战
已经有4人回复
-
DMF怎么除水
已经有13人回复
-
24-25申博
已经有7人回复
-
酰胺合成
已经有6人回复
-
投稿某一SCI期刊，论文被拒绝，是否可以再次投稿到这个期刊？
已经有8人回复
-
七级岗，到手不到8000，带本科论文12教学分/生，额定教学分280
已经有28人回复
-
找工作千万不能来差学校
已经有13人回复
-
双非本科，学生又懒又笨。双一流本科好一些？入职门槛和考核要求比双非本科高多少？
已经有3人回复
-
想请问下IEEE Trans最后一轮，你们等了多久？
已经有8人回复
-
长烷基链氰基水解成羧酸 求助
已经有9人回复
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
nono2009: 金币-50, 屏蔽内容, 违规存档, 禁止回帖广告 2018-04-19 10:44:16
nono2009: 金币-50, 屏蔽内容, 违规存档, 禁止回帖广告 2018-04-19 10:44:16
|
本帖内容被屏蔽 |
2楼2018-04-09 23:14:43
1014511134
专家顾问 (正式写手)
-
专家经验: +82 - 应助: 102 (高中生)
- 金币: 5130.6
- 散金: 20
- 红花: 16
- 帖子: 421
- 在线: 182.1小时
- 虫号: 3220970
- 注册: 2014-05-20
- 性别: GG
- 专业: 应用地球物理学
- 管辖: 程序语言
3楼2018-04-10 16:39:23
40