北航数值分析第三次作业.docVIP

数值分析第三次作业

算法的设计方案

对Fredholm积分方程，用迭代法进行求解：

，其中

对于公式中的积分部分用数值积分方法。

复化梯形积分法，取2601个节点，取迭代次数上限为50次。实际计算迭代次数为18次，最后算得误差为。

复化Simpson积分法，取迭代次数上限为50次，取2*41+1，即83个节点时能满足精度要求。实际计算迭代次数为17次，最后的误差为10。

Guass积分法选择的Gauss—Legendre法，取迭代次数上限为50次，直接选择8个节点，满足精度要求。实际计算迭代次数为24次，最后算得误差为。

全部源程序

moduleintegral

implicitnone

contains

!//////////复化梯形

subroutinetrapezoid(m)

implicitnone

integer::i,j,k,m

real*8::x(m+1),u(m+1)

real*8::sum,sum1,g,r,h

h=2./m

doi=1,m+1

x(i)=-1.+(i-1)*h

enddo

dok=1,50

doi=1,m+1

sum1=0.

g=dexp(x(i)*4.)+(dexp(x(i)+4.)-dexp(-4.-x(i)))/(x(i)+4.)

doj=2,m

sum1=sum1+dexp(x(i)*x(j))*u(j)

enddo

sum=h/2.*(dexp(x(i)*-1.)*u(1)+dexp(x(i)*1.)*u(m+1)+2*sum1)

u(i)=g-sum

enddo

r=h/2.*((dexp(x(1)*4)-u(1))**2+(dexp(x(m+1)*4)-u(m+1))**2)

doi=2,m

r=r+h*(dexp(x(i)*4)-u(i))**2

enddo

if(dabs(r)=e)exit

enddo

write(*,*)k

open(1,file=trapezoid.txt)

doi=1,m+1

write(1,(3(f18.12)))x(i),u(i),dexp(x(i)*4.)

enddo

write(1,(4x,a2,e9.2))r=,r

close(1)

return

endsubroutinetrapezoid

!///////////复化simpson

subroutinesimpson(m)

implicitnone

integer::i,j,k,m

real*8::x(2*m+1),u(2*m+1)

real*8::sum,sum1,sum2,g,r,h

h=2./(2.*m)

doi=1,2*m+1

x(i)=-1.+(i-1)*h

enddo

dok=1,50

doi=1,2*m+1

sum1=0.

sum2=0.

g=dexp(x(i)*4.)+(dexp(x(i)+4.)-dexp(-4.-x(i)))/(x(i)+4.)

doj=1,m

sum1=sum1+dexp(x(i)*x(2*j))*u(2*j)

enddo

doj=1,m-1

sum2=sum2+dexp(x(i)*x(2*j+1))*u(2*j+1)

enddo

sum=h/3.*(dexp(x(i)*-1.)*u(1)+dexp(x(i)*1.)*u(2*m+1)+4*sum1+2*sum2)

u(i)=g-sum

enddo

r=h/3.*((dexp(x(1)*4)-u(1))**2+(dexp(x(2*m+1)*4)-u(2*m+1))**2)

doi=1,m

r=r+4.*h/3.*(dexp(x(2*i)*4)-u(2*i))**2

enddo

doi=1,m-1

r=r+2.*h/3.*(dexp(x(2*i+1)*4)-u(2*i+1))**2

enddo

if(dabs(r)=e)exit

enddo

write(*,*)k

open(2,file=simpson.txt)

doi=1,2*m+1

write(2,(3(f18.12)))x(i),u(i),dexp(x(i)*4.)

enddo

write(2,(4x,a2,e9.2))r=,r

close(2)

return

endsubroutinesimpson

!///////////Gauss_Legendre法

subroutineGauss

imp