Табулирование функции двух переменных

#include <iostream>
#include <cmath>
 
double F_x_y(double x, double y)
{
    if( x >= y ) return log( x / (1 + y));
     else
    if( y > x ) return exp(-(x + y)) * ((1 + x) / (1 +y));
}
 
using namespace std;
 
int main()
{
  double x_start = 0.2;       double y_start = 0.0;
  double x_end   = 0.6;       double y_end   = 0.4;
  double x_step  = 0.1;       double y_step  = 0.05;
  
  for(double x1 = x_start; x1 <= x_end; x1+=x_step)
   {
     for(double y1 = y_start; y1 <= y_end; y1+=y_step)   
      {
        cout << "F(" << x1 << "," << y1 << ")= " << F_x_y(x1 , y1) << endl;  
      }
   }
     
 return 0;
}

double F_x_y(double x, double y)
{
    if( x >= y ) return log(fabs(x / (1 + y)));
     else
    if( y > x ) return exp(-fabs(x + y)) * ((1 + x) / (1 + y));
}

double F_x_y(double x, double y)
{
    if( x >= y ) return log(fabs(x / (1. + y)));    
    return (1. + x) / (1. + y) * exp(-fabs(x + y));
}

#include <iostream>
#include <iomanip>
#include <cmath>
int main()
{
    const double eps = 1.0e-10;
    double x, y;
    std::cout << "x\ty\tf(x,y)" << std::endl;
    for (double x = 0.2; x <= 0.6 + eps; x += 0.1)
    {
        std::cout << std::endl;
        for (double y = 0.0; y <= 0.4 + eps; y += 0.05)
        {
            double f;
            if (x >= y)
                f = log(x / (1 + y));
            else
                f = ((1 + x) * exp(-(x + y))) / (1 + y);
            std::cout << std::fixed << std::setprecision(2) << x << "\t" << y
                << "\t" << std::setprecision(3) << f << std::endl;
        }
    }
    return 0;
}