Вычислить значения интервально заданных функций

#include <iostream>
#include <math.h>
using namespace std;
 
int main()
{
    double y, x;
    for (double a = 0, b = -2; a <= 5, b <= 4; a += 0.5, b += 0.1) {
        if (b < 0) {
            x = pow(a, 5) + 6 * b;
            cout << "a = " << a << ", b = "<< b << ", x = "<< x;
        }
        else if (b >= 0 && b <= 1.5 && a != b) {
            x = 23 * pow(a, 2) + ((45 * a) / (a - b));
            cout << "a = " << a << ", b = " << b << ", x = " << x;
        }
        else if (b > 1.5 && b < 2.5) {
            x = pow(a, 8) + 2 * b * pow(a, 2) + pow(b, 8);
            cout << "a = " << a << ", b = " << b << ", x = " << x;
        }
        else if (b >= 2.5) {
            x = pow(a, 4) - 9 * b + 21 * a * b;
            cout << "a = " << a << ", b = " << b << ", x = " << x;
        }
 
        if (x > 2.3) {
            y = pow(pow((a - 24), 26) - 45 * x, 1 / 25);
            cout << ", y = " << y << endl;
        }
        else if (x <= 2.3) {
            y = pow((a - pow(x, 2)) / 2, -23 / 21);
            cout << ", y = " << y << endl;
        }
    }
    return 0;
}

#include <iostream>
#include <cmath>
using namespace std;
double x(double a,double b)
{
    if(b<0)return pow(a,5)+6.*b;
    if(b<=1.5)
    {
        if(a!=b)
            return 23.*a*a+45.*a/(a-b);
        return 0;
    }
    if(b<2.5)
        return pow(a,8)+2.*b*a*a+pow(b,8);
    return pow(a,4)-9.*b+21.*a*b;
}
double y(double a,double b)
{
    double xt=x(a,b);
    if(xt<=2.3)
    {
        double t=a-xt*xt/2.;
        if(t<0)return 0;
        if(t==0)return 1.;
        return pow(t,-23./21.);
    }else
    {
        double t=pow(a-24.,26)-45.*xt;
        if(t<0)return 0;
        if(t==0)return 1.;
        return pow(t,1./25.);
 
    }
}
int main()
{
    for(double a=0;a<=5.;a+=0.5)
    {
        for(double b=-2;b<=4.;b+=0.1)
        {
            cout<<a<<'\t'<<b<<'\t'<<x(a,b)<<'\t'<<y(a,b)<<endl;
        }
    }
    return 0;
}

if(t<-eps) return -pow(-t,1./25.);