作者sjgau (sjgau)
看板Fortran
标题Re: [心得] fortran 做辛普森积分
时间Thu Sep 11 11:55:20 2008
请问你发表的 code, 有没有实际 验证过。
我用 肉眼观察的结果是:有很大的,很明显的错误。
在 do- loop 里面的
ll= 6.0 - ll
实在是 说不过去,无法解释。
请问你 发表的code, 有没有实际 验证过。
我所发表的 code, 一定验证过,才敢发表。
这是 科学精神。
※ 引述《DrStein (交换关联)》之铭言:
: ※ 引述《sjgau (sjgau)》之铭言:
: : ! do the simpson integrated
: : sum= 0.0
: : ! for i=1 to (ct-2) step 2 do ...
: : i= 1
: : do while (i .LE. (ct-2))
: : sum= sum + y(i) + 4*y(i+1) + y(i+2)
: : i= i + 2
: : end do
: 这地方有较好一点的写法:
: real*8 ll
: ll = 4.0
: sum = 2.0*y(1)
: i = 2
: do i = 2, ct-1
: sum = ll*y(i) + sum
: ll = 6.0 - ll
: enddo
: sum = 2.0*y(ct) + sum
: sum=sum/3.0
: ====
: 全加起来再除3.0 和2/3 4/3 2/3...差再那,我真的不知道。
: 反正老code都这样写,就学起来吧。
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 118.231.134.131
1F:推 noway501:? 4 2 4 2...跳 09/11 12:26
2F:→ sjgau:1, 4, 2, 4, 2, 4, . . . ., 2, 4, 1 09/11 12:32
3F:→ noway501:我想它应该只是把中间的地方改写而已 09/11 12:44
4F:推 DrStein:该敢成 sum/6.0 打错 sorry 09/11 13:31
5F:→ DrStein:喔 没错 是 /3 没错 09/11 13:32
6F:→ DrStein:是 1 4 2 4 2 4 ......2 4 1 跳 09/11 13:33
8F:→ sjgau:数值积分(Numerical Integration) 09/11 14:39