大二下数值计算chapter212积分.pptxVIP

计算物理;带电细丝产生的静电势 ------定积分的近似计算;引入符号f(x)，令:;;想法：用插值函数的积分代替原函数的积分 ;;梯形公式;误差：;;抛物线积分;静电势问题; DIMENSION F(101) READ(*,40)A,C,X0,Y0,N IF(N.GT.100.OR.N.LT.3)STOP IF(N.EQ.2*(N/2))STOP H=2.0*A/FLOAT(N-1) X=-A V=0. DO I=1,N F(I)=C/SQRT((X-X0)**2+Y0*Y0) X=X+H ENDDO DO J=2,N-2,2 V=V+(F(J-1)+4.0*F(J)+F(J+1))*H/3.0 ENDDO WRITE(*,50)V,X0,Y0 GOTO 10 FORMAT(4F10.5,I4) FORMAT(8X,3F10.5) END;;