anchor.diff:  
       x=1.01; y=0.5
       passing = passing.and.test_poly(x,y,poly,.False.)
  -
   
       poly(1,:) = [0.5,0.0,0.5,1.0]
       poly(2,:) = [0.0,0.5,1.0,0.5]
  
positive.code:     endif
  
  
      poly(1,:) = [0,0,1,1]
      poly(2,:) = [0,1,1,0]
  
  
negative.code:     implicit none
  
      real :: x,y
      real :: poly(2,4), large_poly(2,8), triangle(2,3)
      real :: x0,y0,x1,y1
      integer :: i,j
  
      logical :: passing
      passing=.True.
  