Программа, генерирующую координаты точек прямоугольника, равномерно заполненного точками и включающего заданную фигуру

#include <iostream> 
using namespace std;
 
double fi1(double x, double r)
{
    return sqrt(pow(r, 2) - pow((x + r), 2));
}
 
double fi2(double x, double r)
{
    return -sqrt(pow(r, 2) - pow((x - r), 2));
}
 
int main()
{
    setlocale(LC_ALL, "Russian"); // подключение русского языка в консоли 
 
    double R,zx,zy,x1,x2,y1,y2;
 
    cout << "Введите радиус полуокружности: "; cin >> R;
    cout << "Введите координаты зрачка X Y: "; cin >> zx; cin >> zy;
 
    double step = R / 20;       // шаг по оси 
 
    cout << "\n\n\nПЕРВАЯ ФУНКЦИЯ\n\n\n";
    for (x1 = -2*R; x1 < 0; x1 +=  step)
    {
        for (y1 = 0; y1 <= R; y1 += step)
        {
            if (fabs(fi1(x1,R) >= y1)) // проверка условия попадания в область 
            {
                cout << "[ + ] Точка с координатами x=" << x1 << ", y=" << y1 << " попадает в область" << endl;
            }
            else
            {
                cout << "[ - ] Точка с координатами x=" << x1 << ", y=" << y1 << " не попадает в область" << endl;
            }
            if (x1 == zx && y1 == zy) 
                cout << "===========================\n Попадение в зрачок [" << x1 << "][" << y1 <<"]\n===========================\n";
        }
    }
    cout << "\n\n\nВТОРАЯ ФУНКЦИЯ\n\n\n";
    for (x2 = 0; x2 <= 2*R; x2 += step)
    {
        for (y2 = 0; y2 >= -R; y2 -= step)
        {
            if (fabs(fi2(x2,R) <= y2)) // проверка условия попадания в область 
            {
                cout << "[ + ] Точка с координатами x=" << x2 << ", y=" << y2 << " попадает в область" << endl;
            }
            else
            {
                cout << "[ - ] Точка с координатами x=" << x2 << ", y=" << y2 << " не попадает в область" << endl;
            }
            if (x2 == zx && y2 == zy)
                cout << "===========================\n Попадение в зрачок [" << x2 << "][" << y2 << "]\n===========================\n";
        }
    }
 
    cin.get();
}

#include <iostream>
#include <math.h>
using namespace std;
int main()
{
    setlocale(LC_ALL, "Russian"); // подключение русского языка в 
    //консоли 
    double x = 0.0, y = 0.0; // координаты по оси Х и У 
    double step = 0.1; // шаг по оси
    const double r = 1.5;
 
    for (x = -1.5; x <= 1.5; x = x + step)
    {
        for (y = -1.5; y <= 1.5; y = y + step)
        {
            if (((x * y) < 0) && (y < sqrt(pow(r, 2) - pow(x + r, 2))) && (y > (-sqrt(pow(r, 2) - pow(x + r, 2))))) // проверка условия попадания в область 
            {
                cout << "Точка с координатами x=" << x << ", y=" << y << " попадает в область" << endl;
            }
            else
            {
                cout << "Точка с координатами x=" << x << ", y=" << y << " не попадает в область" << endl;
            }
        }
    }
    system("pause");
}

#include <iostream>
#include <math.h>
using namespace std;
int main()
{
    setlocale(LC_ALL, "Russian"); // подключение русского языка в 
    //консоли 
    double x = 0.0, y = 0.0; // координаты по оси Х и У 
    double step = 0.1; // шаг по оси
    double x1, x2, y1, y2, r;
    cout << "Введите координаты прямоугольника по оси х, х1: "; cin >> x1; cout << endl;
    cout << "Введите координаты прямоугольника по оси х, х2: "; cin >> x2; cout << endl;
    cout << "Введите координаты прямоугольника по оси y, y1: "; cin >> y1; cout << endl;
    cout << "Введите координаты прямоугольника по оси y, y2: "; cin >> y2; cout << endl;
    cout << "Введите радиус: "; cin >> r; cout << endl;
 
    for (x = x1; x <= x2; x = x + step)
    {
        for (y = y1; y <= y2; y = y + step)
        {
            if (((x * y) < 0) && (y < sqrt(pow(r, 2) - pow(x + r, 2))) && (y > (-sqrt(pow(r, 2) - pow(x + r, 2))))) // проверка условия попадания в область 
            {
                cout << "Точка с координатами x=" << x << ", y=" << y << " попадает в область" << endl;
            }
            else
            {
                cout << "Точка с координатами x=" << x << ", y=" << y << " не попадает в область" << endl;
            }
        }
    }
    system("pause");
}