24小时热门版块排行榜

>论坛更新日志 (3987)
>虫友互识 (1252)
>文献求助 (329)
>导师招生 (216)
>休闲灌水 (192)
>考博 (164)
>教师之家 (103)
>硕博家园 (96)
>博后之家 (93)
>论文投稿 (88)
>找工作 (79)
>基金申请 (76)
>论文道贺祈福 (70)
>考研 (65)
>公派出国 (57)
>绿色求助(高悬赏) (46)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

夜月丨噩梦

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 694.8
散金: 50
红花: 1
帖子: 20
在线: 17.4小时
虫号: 4040736
注册: 2015-08-28
性别: GG
专业: 半导体晶体与薄膜材料

[求助] 为什么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文件里面啥数据都没有

回复此楼

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有69人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有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人回复

1楼 2015-12-29 21:04:08

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

fmying

银虫 (小有名气)

应助: 29 (小学生)
金币: 618
散金: 50
红花: 6
帖子: 127
在线: 113.6小时
虫号: 1081069
注册: 2010-08-24
性别: GG
专业: 理论和计算化学

引用回帖:

8楼: Originally posted by 夜月丨噩梦 at 2016-01-04 09:18:34
大神，那应该把dett改成多少啊？我改了几个输出的结果都不对
...

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)，这样就可以在少算很多的情况下得到结果？
这些都是你自己需要去考虑的。