| 查看: 3353 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
为什么Fortran调试无错误运行不出结果 已有2人参与
|
|||
|
program main parameter n=600 double precision u(0:n),u1(0:n),x0,x1,x,t0,t1,h,dett,lambda,y(0:n) integer i,j x0=0.0d0 x1=6.0d0 t0=0.0d0 lambda=0.3d0 t1=0.5d0 h=0.01d0 dett=0.003d0 x=x0 !do i=0,n,1 !if(i<=250)then !y(i)=1.0 !else ! y(i)=0 !end if ! end do do i=0,n,1 if(x>=1 .and. x<=3) then u(i)=1.0d0 else u(i)=0.0d0 end if x=x+h end do open(1,file="yingfeng.dat"  t=t0+dett x=x0 do while(t<=t1) x=x0+h do j=1,n,1 u1(j)=u(j)-3*lambda*(u(j)-u(j-1)) if(abs(t-0.5)<=0.0001)then write(1,*)x,u1(j),y(j) end if x=x+h end do do j=1,n-1 u(j)=u1(j) u(0)=1 u(1)=0 end do t=t+dett end do stop end 求大神指导，为啥这个程序运行后yingfeng.dat文件里面啥数据都没有 |
回复此楼» 猜你喜欢
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
物理学I论文润色/翻译怎么收费?
已经有232人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有1人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Intel fortran与visual fortran
已经有11人回复
- FORTRAN
已经有0人回复
- Fortran
已经有7人回复
- 已经安装studio，为什么还是不能安装Fortran呢？
已经有0人回复
- 开始学习FORTRAN，结果杀毒软件总是认为生成的exe程序是木马，要求删掉，为什么呢？
已经有7人回复
- fortran
已经有4人回复
- Fortran77 代码调试中提示有Syntax error，但却不知道为什么出错。望高人指点！
已经有4人回复
- fortran循环求助
已经有3人回复
- Fortran 求解复数矩阵SVD时，用到cgesvd，结果与matlab不同，为什么？
已经有15人回复
- fortran里的intrisic function: precision()x结果表示什么意思
已经有5人回复
- fortran 程序修改后改回来，计算结果不同？？？为什么会这样
已经有2人回复
- data 给fortran 的二维数组如何命名？为什么老出错，解决送金币
已经有1人回复
- fortran输出结果为0
已经有3人回复
- fortran 程序有write和没有write两种情况计算的结果为什么会不一样?
已经有14人回复
- fortran输出求助
已经有6人回复
- 帮忙把公式转化为FORTRAN子程序
已经有6人回复
- fortran中计算数据写入文本为 NAN的问题
已经有4人回复
- 刚开始学Fortran 帮忙看一下这个程序 为什么算不出结果
已经有3人回复
- fortran断点
已经有2人回复
- fortran 运行问题
已经有2人回复
- 【求助】请问有什么软件可以将C语言代码转换为fortran吗？
已经有7人回复
|
dett相当于步长，t0+dett是初值，0.0001是阈值。 这些都要根据你要计算解决的问题来啊。 程序的关键是解决问题，所以要先理解透问题和你的需求。 按照t<=t1和abs(t-0.5)<=0.0001来看，你是希望得到最后一步或几步的结果。 那么你希望得到哪一步/几步？这几步的步长是希望大一点粗一点，还是小一点精一点？ 比如，如果其它条件都不变，只是希望得到最后5步的结果的话，甚至可以步长设为0.00002，这样最后会有6步的结果因符合(abs(t-0.5)<=0.0001)的条件而输出。 但是你要计算太多遍了，这样做是否值得？ 如果不需要如此精密的话，是否可以修改阈值为(abs(t-0.5)<=0.1)，这样就可以在少算很多的情况下得到结果？ 这些都是你自己需要去考虑的。 |
9楼2016-01-04 17:02:08
bennett_14
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 7204.8
- 红花: 4
- 沙发: 1
- 帖子: 932
- 在线: 421.5小时
- 虫号: 1073877
- 注册: 2010-08-13
- 性别: GG
- 专业: 计算机科学
2楼2016-01-02 16:09:59
bennett_14
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 7204.8
- 红花: 4
- 沙发: 1
- 帖子: 932
- 在线: 421.5小时
- 虫号: 1073877
- 注册: 2010-08-13
- 性别: GG
- 专业: 计算机科学
3楼2016-01-02 16:35:26
4楼2016-01-02 17:51:51
