作者DrStein (交换关联)
看板Fortran
标题Re: [问题] 球的积分
时间Wed Apr 8 00:37:43 2009
建议直接用球座标做
int i, j , k
real*8 a = 1.0
real*8 b = 0.0
real*8 r, theta
real*8 dr, dtheta, dphi
real*8 V
real*8 pi
pi = 4.0*atan(1.0)
integer n = 100 //设切100片
dr = (a-b)/n
dtheta = pi/n
dphi = (2*pi)/n
V = 0
do k = 1, n
do j = 1, n
do i = 1, n
r = i*dr
theta = j*dtheta
V = (r*sin(theta)*dphi)*(r*dtheta)*(*dr) + V
enddo
enddo
enddo
※ 引述《b906161 (...)》之铭言：
: 是这样子的
: 要算半径为1的球的积分 先算一个象限後再乘以4和8
: 以长方形法来算圆的话就用do回圈(a=0 b=1)
: deltax=(b-a)/real(n)
: do i=1,n-1
: x=x+deltax
: y=f(x)
: sum=sum+y
: end do
: 而球的话我是想成
: do i=1,n-1
: do j=1,n-1
: x=x+deltax
: y=f(x)*f(x)
: sum=sum+y
: end do
: end do
: 不过答案错了
: 不知道错在哪
: 感谢指教~
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 59.117.66.246
1F：→ DrStein:啊。。注用到C++的方式。。 04/08 00:38