anchor.diff:    integer, intent(OUT) :: my_point_index, target_point_index
     real(kind=8) :: xa1, ya1, xa2, ya2
     type(coef_type), pointer :: c_ptr
  -  integer :: target_polygon_index
  -  integer :: i1, i2
  +  !integer :: target_polygon_index
  +  integer :: i1
  +  !integer :: i2
   
   
     c_ptr => cb
  
positive.code:       write(0,*) "cal_coefficient num_of_target error "
        stop
      end if
      if (target_polygon%mask) then
        call set_target_polygon_coef1_by_num(my_polygon, 1, 1.d0)
      else
  
negative.code:   integer, intent(IN) :: div_num
    real(kind=8) :: point_pos_x, point_pos_y
    integer :: i, j, p
    type(polygon_type), pointer :: rect_polygon, target_polygon
    real(kind=8), dimension(div_num+1) :: x12, x23, x43, x14, y12, y23, y43, y14
    real(kind=8) :: x_pos(4), y_pos(4), x1, y1, x2, y2, x3, y3
    real(kind=8) :: px1, px2, px3, px4, py1, py2, py3, py4
  
    ! set rectangular polygon
    if (pa%num_of_point == 4) then ! rectangular polygon
  