      program Taylor
      implicit real*8 (a-h,o-z)
      parameter (N=10)
      dimension x(0:N),y(0:N),dydx(0:N),d2ydx2(0:N),d3ydx3(0:N)
      
      x(0)=0d0
      x(N)=1d0
      h=(x(N)-x(0))/N
      print*,'N   =',N
      print*,'x(0)=',x(0)
      print*,'x(N)=',x(N)
      print*,'h   =',h
      print*
      
      do i=0,N
         x(i)=i*h+x(0)
      enddo
      
      y(0)=1d0
      print '(3x,''i'',9x,''x'',9x,''y'',6x,''y_an'')'
      print '(i4,3f10.5)',0,0d0,y(0),1d0
      do i=0,N-1
         dydx(i)=x(i)+y(i)
c         d2ydx2(i)=?
c         d3ydx3(i)=?
         y(i+1)=y(i)+h*dydx(i)
c         y(i+1)=?
c         y(i+1)=?
c         y(i+1)=?
         print '(i4,3f10.5)',i+1,x(i+1),y(i+1)
     &        ,2d0*dexp(x(i+1))-x(i+1)-1d0
      enddo
      
      pause
      end
