implicit none
      real :: x,y,dy,eps=0.00001
      integer :: kaisuu

      x=1.0
      do kaisuu=1,1000
       y=x-cos(x)
       dy=1.0+sin(x)
        if(abs(dy)<eps)then
          write(*,*) 'NG'
          stop
        endif
       x=x-y/dy
      write(*,'(1x,A,F10.6)') 'x=',x
       if(abs(y)<eps) exit
       if(abs(y/dy)<eps) exit
      enddo

      write(*,'(1x,A,F10.6)') 'x=',x
      write(*,'(1x,A,F10.6)') 'y=',x-cos(x)
 
        
       end