Нахождение точки внутри прямоугольника

#include <iostream>
#include <stdlib.h>
using namespace std;
 
bool inside(const int&px  , const int&py, const int&ax  , const int&ay, const int&bx  , const int&dy){
    if((px < bx)&&(px > ax)&&(py < dy)&&(py > ay)){
        return true;
    }
    return false;
 
}
int main()
{
    int px  , py;
    cout << "Koordinati tochki x i y: ";
    cin >> px>>py;
 
    int ax, ay, bx, by, cx, cy, dx, dy;
    /*
        D--------------------C
        |                    |
        |                    |
        |                    |
        |                    |
        A--------------------B
    */
 
    cout << "Koordinati prjamougoljnika ax i ay: ";
    cin >> ax>>ay;
    cout << "Koordinati prjamougoljnika bx i by: ";
    cin >> bx>>by;
    cout << "Koordinati prjamougoljnika cx i cy: ";
    cin >> cx>>cy;
    cout << "Koordinati prjamougoljnika dx i dy: ";
    cin >> dx>>dy;
    //(x1-x3)*(y2-y3)-(x2-x3)*(y1-y3);
    if(inside(px,py,ax,ay,bx,dy))
        cout << "Tochka prinadlezit!!!";
 
}

double Dlina(int x1,int x2,int y1,int y2){  // функция выч длины по координатам двух точек
          double Dlin=sqrt(pow((x2-x1),2)+pow((y2-y1),2); // корень суммы квадратов разницы координат))
          return Dlin;
 
}
 
double PLOCHAD( double a, double b, double c){ // функц вычисления площади по длинам сторон треуг.
        double p, s;
        p=(a+b+c)/2;
s=sqrt(p*(p-a)*(p-b)*(p-c))
   
  return s;
 
}

...............
#define _USE_MATH_DEFINES
#include <cmath>
...............
struct point{ double x,y,a; };
void main(int argc,char* argv[])
{
    setlocale(LC_ALL,"Rus");
    point ap[]={{4,1,0},{8,8,0},{2,5,0},{10,4,0}};
    double xmin,xmax,ymin,ymax;
    xmin=xmax=ap[0].x;
    ymin=ymax=ap[0].y;
    for(int i=1; i<4; i++)
    {
        if(xmin>ap[i].x)xmin=ap[i].x;
        if(xmax<ap[i].x)xmax=ap[i].x;
        if(ymin>ap[i].y)ymin=ap[i].y;
        if(ymax<ap[i].y)ymax=ap[i].y;
    }
    double x=(xmax+xmin)/2;
    double y=(ymax+ymin)/2;
    for(int i=0; i<4; i++)
    {
        ap[i].a=atan2(ap[i].y-y,ap[i].x-x);
        if(ap[i].a<0)ap[i].a+=2*M_PI;
    }
    for(int i=3; i; i--)
    for(int j=0; j<i; j++)
        if(ap[j].a>ap[i].a)
        {
            point t=ap[j];
            ap[j]=ap[i];
            ap[i]=t;
        }
    //for(int i=0; i<4; i++) cout<<ap[i].x<<" "<<ap[i].y<<endl;
    cout<<"x y через пробел: ";
    cin>>x>>y;
    bool in=true;
    if(ap[0].y==ap[1].y)
    {
        if(x<ap[1].x || x>ap[0].x || y<ap[2].y || y>ap[0].y) in=false;
    }
    else
    {
        //y=a*x+b
        int i;
        for(i=0; i<2; i++)
        {
            double a=(ap[i].y-ap[i+1].y)/(ap[i].x-ap[i+1].x);
            double b0=ap[i].y-a*ap[i].x;
            double b1=ap[i+2].y-a*ap[i+2].x;
            double b=y-a*x;
            if(b<b0 && b<b1 || b>b0 && b>b1) break;
        }
        if(i<2) in=false;
    }
    cout<<(in?"внутри\n":"снаружи\n");

#include <iostream>
#include <cmath>
 
using namespace std;
 
int main()
{
    int x,x1,x2,x3,x4,y,y1,y2,y3,y4;
    double stt,a,b,att,btt,ctt,ptt,sp,att1,btt1,ctt1,ptt1,stt1,att2,btt2,ctt2,ptt2,stt2,att3,btt3,ctt3,ptt3,stt3;
    cin>>x>>y>>x1>>y1>>x2>>y2>>x3>>y3>>x4>>y4;
    b=sqrt(pow(x2-x1,2)+pow(y2-y1,2));
    a=sqrt(pow(x3-x2,2)+pow(y3-y2,2));
    sp=a*b;
    //cout<<a<<" "<<b<<" "<<sp<<endl;
    att=sqrt(pow(x2-x,2)+pow(y2-y,2));
    btt=a;
    ctt=sqrt(pow(x3-x,2)+pow(y3-y,2));
    ptt=(att+btt+ctt)/2;
    stt=sqrt(ptt*(ptt-att)*(ptt-btt)*(ptt-ctt));
    //cout<<att<<" "<<btt<<" "<<ctt<<" "<<ptt<<" "<<stt<<endl;
    att1=ctt;
    btt1=b;
    ctt1=sqrt(pow(x4-x,2)+pow(y4-y,2));
    ptt1=(att1+btt1+ctt1)/2;
    stt1=sqrt(ptt1*(ptt1-att1)*(ptt1-btt1)*(ptt1-ctt1));
    //
    att2=ctt1;
    btt2=a;
    ctt2=sqrt(pow(x1-x,2)+pow(y1-y,2));
    ptt2=(att2+btt2+ctt2)/2;
    stt2=sqrt(ptt2*(ptt2-att2)*(ptt2-btt2)*(ptt2-ctt2));
    //
    att3=ctt2;
    btt3=b;
    ctt3=att;
    ptt3=(att3+btt3+ctt3)/2;
    stt3=sqrt(ptt3*(ptt3-att3)*(ptt3-btt3)*(ptt3-ctt3));
    cout<<stt+stt1+stt2+stt3<<" "<<sp;
}