Точка внутри или на границе многоугольника

Я использую алгоритм в http://www.ecse.rpi.edu/~wrf/Research/Short_Notes/pnpoly.html,
но когда входная точка находится на границе, этот алгоритм дает мне неверный результат. Может кто-нибудь помочь мне с точкой в ​​граничном случае?
Любая помощь приветствуется.

Это основная функция

#include <iostream>
#include <Polygon.h>
using namespace std;

int main()
{
vector<Point> v;
//v.push_back(make_pair(3.0,3.0));
v.push_back(make_pair(1.0,1.0));
v.push_back(make_pair(1.0,5.0));
v.push_back(make_pair(5.0,5.0));
v.push_back(make_pair(5.0,1.0));
Polygon *p = new Polygon(v);
cout << "A: " << p->IsInside(make_pair(1.0,3.0)) << endl;
cout << "B: " << p->IsInside(make_pair(3.0,1.0)) << endl;
cout << "C: " << p->IsInside(make_pair(5.0,3.0)) << endl;
cout << "D: " << p->IsInside(make_pair(3.0,5.0)) << endl;
delete p;
return 0;
}

Это функция проверки

bool Polygon::IsInside(Point p)
{
/*determine whether a point is inside a polygon or not
*  polygon's vertices need to be sorted counterclockwise
*  source :
*      http://www.ecse.rpi.edu/~wrf/Research/Short_Notes/pnpoly.html
*/
bool ans = false;
for(size_t c=0,d=this->vertices.size()-1; c<this->vertices.size(); d=c++)
{
if( ((this->vertices[c].y > p.y) != (this->vertices[d].y > p.y)) &&
(p.x < (this->vertices[d].x - this->vertices[c].x) * (p.y - this->vertices[c].y) /
(this->vertices[d].y - this->vertices[c].y) + this->vertices[c].x) )
ans = !ans;
}
return ans;
}