Вычислить значения функции на заданном интервале изменения аргумента. С++

#include <iostream>
#include <cmath>
 
using namespace std;
 
 
int main()
{
 const double b = 0.01078;
 const double deltax = 0.1;
 const double a[] = {3.89, 1.3, 0.44};
 
 double y = 0.0;
 
 for(double x = -0.1; x <= 1.1; x+=deltax)
  {
   if(x < 0.2)
   {
    y = (a[0] * a[0] + b * b) * x *x + sqrt((x * x * x *x) + 1);
   }
   else
   if(x == 0.2)
    {
     y = 5 * x * x + a[1] * a[1] * sqrt((x * x) + 2);
    }
   else y = a[2] * x + exp(x - 0.5);
 
   cout << "y (" << x << ")=" << y << endl;
  }
return 0;
}

#include <iostream>
#include <iomanip>
#include <cmath> 
using namespace std;
 
void Tabulation (const double a, const double b)
{
    cout <<"if a="<<a<<"  b="<<b<<endl;
    cout <<"  x    y\n";
    double dx=0.1, y;
    for(double x = -0.1; x <= 1.1; x+=dx)
      {
      if(x < 0.2) y = (a*a + b*b)*x*x + sqrt(x*x*x*x + 1);     
      else if(x == 0.2) y = 5*x*x + a*a* sqrt(x*x + 2);      
      else y = a*x + exp(x - 0.5); 
      cout <<setw(4)<< x <<"  "<< y << endl;
      } 
      cout << endl;
}
 
int main()
{
    const double a1 = 3.89, a2 = 1.3, a3 = 0.44, b=10.78e-3;
    Tabulation (a1,b);
    Tabulation (a2,b);
    Tabulation (a3,b);
system("pause");
return 0;
}