北航数值分大作业3.docVIP

数值分析第三次作业 设计方案 对Fredholm积分方程，用迭代法进行求解： ，其中 对于公式中的积分部分用数值积分方法。 复化梯形积分法，取2601个节点，取迭代次数上限为50次。实际计算迭代次数为18次，最后算得误差为r= 0.97E-10。 复化Simpson积分法，取迭代次数上限为50次，取2*41+1，即83个节点时能满足精度要求。实际计算迭代次数为17次，最后的误差为 r= 0.97E-10。 Guass积分法选择的Gauss—Legendre法，取迭代次数上限为50次，直接选择8个节点，满足精度要求。实际计算迭代次数为24次，最后算得误差为r= 0.87E-10。 全部源程序 module integral implicit none contains !//////////复化梯形 subroutine trapezoid(m) implicit none integer :: i,j,k,m real*8 :: x(m+1),u(m+1) real*8 :: sum,sum1,g,r,h real*8 :: e=1.0e-10 h=2./m do i=1,m+1 x(i)=-1.+(i-1)*h end do u=0.02 do k=1,50 do i=1,m+1 sum1=0. g=dexp(x(i)*4.)+(dexp(x(i)+4.