作者DrStein (交换关联)
看板Fortran
标题Re: [问题] 如何用辛普森积分法解这题?
时间Mon May 18 21:23:30 2009
※ 引述《sawwas (saw)》之铭言：
: Q: 将 e^(-x^2) 用辛普森积分法从 0 积到 3
: ( 限定用Fortran语法编译 )
: 因为我从来没学过Fortran 所以只好来此求助
: 拜托各位大大了~~~
切1000格
integer :: i
real*8 :: y
real*8 :: x
real*8 :: dx
real ll
real sum
dx = 3.0/(1000)
ll = 4.0
sum = 1.0*exp(-(dx*1)**2 )*dx
do i = 2, i=999
x = i*dx
y = exp(-x**2 )
sum = ll*y*dx + sum
ll = 6.0 - ll
enddo
sum = 1.0*exp(-(dx*1000)**2 )*dx + sum
sum = sum/3.0
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 61.216.81.152
1F：推 dana79810:哈瞜!!刚好我也是做这一样功课 i因该从1开始才可以唷 05/29 17:34
2F：→ dana79810:因为我也跑了一两个数字才发现有不一样 ，这样就漏掉 05/29 17:35
3F：→ dana79810:F(x1)直接跳到F(x2) 05/29 17:35