南开大学计算物理讲义5.ppt

下载文档 降价啦

7
0
约9.89千字
约 117页
2019-09-14 发布于湖北
举报
版权申诉
保障服务

南开大学计算物理讲义5.ppt

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

Using three function evaluations: symmetry suggests for the interval x0- ? x to x0+?x the points should be evaluated at x0 and x0±?. The Taylor Series expansion of the integral 三点高斯-勒让德公式为积分区间为无限的积分变量代换法 x= -ln(t) dx= -dt/t 6.8 Example of numerical integration C Trapezium rule nx = 1 DO i=1,20 Value = TrapeziumRule(x0,x1,nx) WRITE(6,*)nx,Value,Value - Exact nx = 2*nx ENDDO REAL*8 FUNCTION TrapeziumRule(x0,x1,nx) C=====parameters INTEGER*4 nx REAL*8 x0,x1 C=====functions REAL*8 f C=====local variables INTEGER*4 i REAL*8 dx,xa,xb,fa,fb,Sum dx = (x1 - x0)/DFLOAT(nx) Sum = 0.0 DO i=0,nx-1 xa = x0 + DFLOAT(i)*dx xb = x0 + DFLOAT(i+1)*dx fa = f(xa) fb = f(xb) Sum = Sum + fa + fb ENDDO Sum = Sum * dx / 2.0 TrapeziumRule = Sum RETURN END FUNCTION f(x) C=====parameters REAL*8 x,f f = SIN(x) RETURN END REAL*8 FUNCTION MidpointRule(x0,x1,nx) C=====parameters INTEGER*4 nx REAL*8 x0,x1 C=====functions REAL*8 f C=====local variables INTEGER*4 i REAL*8 dx,xa,fa,Sum dx = (x1 - x0)/Dfloat(nx) Sum = 0.0 DO i=0,nx-1 xa = x0 + (DFLOAT(i)+0.5)*dx fa = f(xa) Sum = Sum + fa ENDDO Sum = Sum * dx MidpointRule = Sum RETURN END REAL*8 FUNCTION SimpsonsRule(x0,x1,nx) C=====parameters INTEGER*4 nx REAL*8 x0,x1 C=====functions REAL*8 f C=====local variables INTEGER*4 i REAL*8 dx,xa,xb,xc,fa,fb,fc,Sum dx = (x1 - x0)/DFLOAT(nx) Sum = 0.0 DO i=0,nx-1,2 xa = x0 + DFLOAT(i)*dx xb = x0 + DFLOAT(i+1)*dx xc = x0 + DFLOAT(i+2)*dx fa = f(xa) fb = f(xb) fc = f(xc) Sum = Sum + fa + 4.0*fb + fc ENDDO Sum = Sum * dx / 3.0 SimpsonsRule = Sum RETURN END REAL*8 FUNCTION GaussQuad(x0,x1,nx) C=====parameters INTEGER*4 nx REAL*8 x0,x1 C=====functions REAL*8 f C=====local variables INTEGER*4 i REAL*8 dx,xa,xb,fa,fb,Sum,dxl,dxr dx = (x1 - x0)/DFLOAT(nx) dxl = dx*(0.5D0 - SQRT(3.0D0)/6.0D0) dxr = dx*(0.5D0 + SQRT(3.0D0)/6.0D0) Sum = 0.0 DO i=0,nx-1 xa = x0 + DFLOAT(i)*dx + dxl xb = x0 + DFLOAT(i)*dx + dxr fa = f(xa) fb = f(xb) Sum = Sum + fa + fb ENDDO Sum = Sum * dx / 2.0 GaussQuad = Sum RETURN EN