第七讲物理学中定积分是数值计算方法.pptVIP

* * 第七讲 Fortran中的数值积分 物理问题的求解，常常涉及到积分，因此本节介绍在Fortran中进行数值积分的方法。 一、基本的数值积分方法 数值求该面积：一般都是将积分区间[a,b]分成有限多个小段（每段宽h=(b-a)/n），然后“近似”求出每小段的“面积”，再把它们加起来，以此作为积分的近似值。 定积分： 表示f(x)正下方与x轴所围的“面积”。 * * 主要有3种对小段面积的近似方法： （1）矩形法 这种方法把每小段用宽h、高f(a+ih)的“矩形”来近似，即： 据此，计算的函数子程序程序如： ！矩形积分子程序： function sjx(f,a,b,n) h=(b-a)/n sjx=0 do i=1,n sjx=sjx+h*f(a+i*h) end do end * * （2）梯形法 每小段用“梯形”---上底f(xi-1)、下底f(xi)、高h，“近似”其面积，即： 于是,得到计算的子程序： !梯形积分子程序： subroutine stx(f,a,b,n,s) h=(b-a)/n s=0.5*h(f(a)+f(b)) do i=1,n-1 s=s+h*f(a+i*h) end do end 注意：利用了x0=a、xn=b，且xi=a+ih。 * * （3）辛卜生（Simpson）法 该四边形的面积为： 这种方法把“每2个小块”用“1个抛物线为顶