作者masto (momo)
看板Programming
标题[问题] FORTRAN问题
时间Sat Mar 1 21:05:08 2008
f=0.5
do i=0,3.49,0.01
fz0=f+0.5*0.01*(((1./sqrt(2.*pi))*exp(-(z(i)**2./2.)))+((1./sqrt(2.*pi))*exp(-(z(i+1)**2./2.))))
fz(i)=fz0+0.5*0.01*(((1./sqrt(2.*pi))*exp(-(z(i+1)**2./2.)))+((1./sqrt(2.*pi))*exp(-(z(I+2)**2./2.))))
fz(I+1)=fz(i)+0.5*0.01*(((1./sqrt(2.*pi))*exp(-(z(i+2)**2./2.)))+((1./sqrt(2.*pi))*exp(-(z(I+3)**2./2.))))
fz(I+2)=fz(i+1)+0.5*0.01*(((1./sqrt(2.*pi))*exp(-(z(i+3)**2./2.)))+((1./sqrt(2.*pi))*exp(-(z(I+4)**2./2.))))
!fz(i)=fz(I+1)
end do
要跑到200可已帮我修改一下写法吗?
最後回到3.49
请指教 多谢
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.124.60.51
1F：→ MOONRAKER:就一直打到i+200就好了，不用别人帮忙 61.224.251.244 03/01 21:52
2F：→ MOONRAKER:才两百行而已！有什麽好怕的！ 61.224.251.244 03/01 21:53
3F：→ masto:恩!看看有没有更好的逻辑阿!卸卸大大分享 140.124.60.51 03/01 22:01
4F：→ MOONRAKER:…你知道do loop的用法吧？ 59.115.215.72 03/01 22:32
5F：→ MOONRAKER:你知道就不用打到i+200，不然就老实点 59.115.215.72 03/01 22:35
6F：→ MOONRAKER:这里面根本无逻辑可言，只是找一找 59.115.215.72 03/01 22:35
7F：→ MOONRAKER:基本指令怎麽用，就这样而已 59.115.215.72 03/01 22:35
8F：→ youling0725:do i=200,3.49,-0.01 <-这样吧 122.116.39.75 03/02 09:48
9F：→ youling0725:2年没碰Fortran 122.116.39.75 03/02 09:49