Принадлежность точки к заданой плоскости

#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<string>
#include<set>
#include<algorithm>
#include<vector>
#include<map>
#include<iomanip>
 
using namespace std;
 
#define DB(a) cout<<#a<<"="<<a<<" ";
#define DBN(a) cout<<#a<<"="<<a<<endl;
 
 
typedef long long ll;
typedef unsigned long long ull;
typedef vector<int> vi;
typedef pair<int, int> pii;
typedef vector<ll> vll;
typedef vector<pii> vpii;
typedef map<int, int> mii;
typedef set<int> seti;
 
//#define min(a,b) ((a)<(b)?(a):(b))
//#define max(a,b) ((a)>(b)?(a):(b))
//#define sqrt(a) (sqrt(1.0*(a)))
//#define sqr(a) ((a)*(a))
//#define all(a) (a.begin(),a.end())
 
 
//const int mod = 1000000007; // const int mod = 1000000009;
//const double EPS = 1E-7;
 
int main() {
    double x, y, fx, fr;
    cin >> x >> y;
    if (x > M_PI / 2) x -= M_PI / 2;
    fx = sin(2 * x);
    fr = asin(y) / 2;
    if ((x >= M_PI / 2 && x <= M_PI / 2) && (fabs(x - fabs(asin(y - 0.0001)) / 2) < 0.0001 || fabs(x - fabs(asin(y + 0.0001) / 2)) < 0.0001)) cout << "ON BORDER" << endl;
    else if (fx + 0.0001 < y || x<-0.0001 || y<-0.0001 || x > M_PI / 2 + 0.0001 || y > 1 + 0.0001) cout << "NO" << endl;
    else if ((x >= 0 && x <= M_PI / 2) && (fabs(y) < 0.0001 || fabs(y - fx) < 0.0001)) cout << "ON BORDER" << endl;
    else if (fabs(fx - y) < 0.0001 || fabs(fr - x) < 0.0001) cout << "ON BORDER" << endl;
    else cout << "YES" << endl;
    return 0;
}