有没有人帮忙解决一下舍入误差的问题？exp(1)的泰了展开

MRW

答案应该是从0-n积分与n-0积分结果是不一样的，可是我的程序里这两种方法是一样的


program main
      implicit none
      real(kind=8)::epsilion,i,j,k,a,n,y1,y2 ,z
      epsilion=1.0D-6
      y1=1.0d0
      y2=0.0d0
      do i=1.0d0,1.0d5,1.0d0
         a=1.0d0
         do n=1.0d0,i,1.0d0
            a=a*(1.0d0/i)
         end do
         if(a>epsilion) then
            continue
         else
            exit
         end if
      end do
      write(*,*)"n=",n
      do j=1.0d0,n,1.0d0
         a=1.0d0
         do k=1.0d0,j,1.0d0
            a=a*(1.0d0/k)
         end do
         y1=a+y1
      end do
      write(*,*)"from 0 to n make sum y1=",y1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      do j=n,1.0d0,-1.0d0
         a=1.0d0
         do k=1.0d0,j,1.0d0
            a=a*(1.0d0/k)
         end do
         y2=a+y2
      end do
      y2=y2+1.0D0
      write(*,*)"from n to 0 make sum y2=",y2
      z=exp(1.0d0)
      write(*,*)"exp(1)=",z
      pause
      end

胡文刚

为啥不一样？

MRW

胡文刚 发表于 2019-9-17 12:16
为啥不一样？

因为存在舍入误差

胡文刚

还好吧，可能没体现出来。
误差问题挺复杂的，不是所有情况都会体现出来。

能得到一致的结果，是个好事。乐吧您就。