作者alan77603 (阿叶)
看板Fortran
标题[问题] 为什跑出来值比课本上值小很多
时间Sun Oct 3 15:20:51 2010
我求的是(8/PI^2)sigma((1/n^2)*Sin(0.5n*PI))Sin(n*PI*x)*exp(-(n*PI)^2)*t)
sigma 底是n=1 到200
0<t<0.01 以0.005为间距
0< x<1 以0.1为间距
而我的程式码是
real*8 an,x,ua,t,u,ti,pi,n
do n=1,200
do x=1,9,0.1
xi=1/10.
do t=1,20,0.005
ti=t/20.
pi=3.1415926
ua=0.0
an=(1/n**2)*sin(0.5*n*pi*pi/180.)*sin(xi*n*pi*pi/180.)*EXP(-((n*pi)**2)*ti)
ua=ua+an
u=(8./pi**2)*ua
write(*,*)u
end do
end do
end do
end
课本上值大概为 0.X...
我跑出来大概0.000000000..............
那麽小
请问各位版上先进
是否小弟程式哪里写错了
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.121.221.176
1F:→ farewell2400:应该是sin函数里多除180. 10/03 15:33
2F:→ gazzy:ua=0.0 要放在do loop上面,u=(8./pi**2)*ua要放在do loop後 10/03 16:48
※ 编辑: alan77603 来自: 140.121.221.176 (10/03 17:20)