Циклические вычислительные процессы. Затабулировать функцию Z=f(x,y) с шагами x и y внутри прямоугольника и области существования данной функции

#include <iostream>
#include <fstream>
#include <stdlib.h>
#include <math.h>
#include <stdio.h>
#include <iomanip>
 
using namespace std;
int main()
{ 
    setlocale(LC_ALL, "Russian"); // функция локализации вызывается только один раз
    cout<<" Задание: Затабулировать функцию Z=f(x,y) с шагами dx и dy внутри прямоугольника и области существования данной функции. Составить универсальную программу   для произвольных значений параметров   a, b, c, d, dx и dy. "<<endl;
    cout<<" "<< endl;
    cout<<" Мое задание: Вид функции: sqrt((x-y)-(y-x)). "<< endl;
    cout<<" Рабочий набор исходных данных: a=-1; b=1; c=-1; d=1; dx=0.2; dy=0.2; "<< endl;
    cout<<" "<< endl;
    cout<<" "<< endl;
    double a, b, c, d,dx, dy, x, y, z;
    cout << "Введите a: ";
    cin >> a;
    cout << "Введите b: ";
    cin >> b;
    if (a <= b)
        {
            cout << "Введите c: ";
            cin >> c;
            cout << "Введите d: ";
            cin >> d;
    if (c <= d)
        {
            cout << "Введите dx: ";
            cin >> dx;
            cout << "Введите dy: ";
            cin >> dy;
            system("cls");
            cout << "Выберите метод вывода (1- на экран, 2 - в файл) ";
            int num;
            cin >> num;
            ofstream out;
            switch (num)
                {
                    case 1:
                        out.open("out.txt");
                        cout  << setprecision(2) << "X = [" << a <<", " << b << "] Y = [" << c << ", " << d << "] dx = " << dx << " dy = " << dy<<endl;
                        //шапка таблицы
                        cout  << "\\       |"<<endl;
                        cout  << " \\ X    |";
                        for (x = a; x <= b;  x+=dx) 
                            {
                                cout <<fixed <<setw(6) << setprecision(2) << ceil(100*x)/100.0;
                            } 
                        cout  <<endl<< "   \\    |"<< endl;
                        cout  << "  Y \\   |" <<endl;
                        cout  << "     \\  |" <<endl;
                        cout  << "      \\ |" <<endl;
                        for (int i =0; i<80; i++)
                        cout<<"-";
                        cout<<endl;
                        for(y = c; y <= d; y+=dy);
                            {
                                //while (abs(x)>0.0001 || abs(y)>0.0001 || ((x/y)-(y-x)>0))                                     
                                //continue;         --------------------------  
                                if(y <0)
                                cout <<fixed<<setw(3) << setprecision(2) << "  "<< ceil(100*y)/100.0<<" |";
                                else
                                cout <<fixed<<setw(4) << setprecision(2) << "  "<< ceil(100*y)/100.0<<" |";
                                if (((x/y)-(y-x))>0)
                                z = sqrt((x/y)-(y/x));
                                if(z <0)
                                cout<<fixed <<setw(7) << setprecision(2) << ceil(100*z)/100.0<<"|";
                                 else
                                cout <<fixed<<setw(7) << setprecision(2) << ceil(100*z)/100.0<<"|";
                                cout <<fixed<<setw(3) << setprecision(2) << "   *   |";
                                out << endl;
 
                            }
 
 
    
                                
 
 
                    break;      
                    case 2: //на эту кнопку в будущем будет происходить сохранение в файл 
                    break;
                    default: cout<< "Данные введены неверно.";
                    break;
                }
            out.close();
        }
            else
            cout << "Некорректный диапазон для Y" << endl;
        }
            else
            cout << "Некорректный диапазон для X" << endl;
    system("pause");
    return 0;
}