Вывести таблицу значений функции двух переменных

#include <iostream>
#include <cmath>
using namespace std;
 
double function()
{
    double a, b, f, amax, bmax, amin, bmin, da, db;
    int i,n;
    for (i = 0; i <= n; i = i + 1) {
        for (a < amax; a = a + da;) {
            for (b < bmax; b = b + db;) {
                if (b > 0)
                {
                    f = (pow((a + pow(-1, i) * b * i), 2)) / (i * sqrt(i + b * b));
                }
                else
                {
                    f = (pow((a + b * i), 2)) / (i * sqrt(i + b * b));
                }
            }
        }
    }
    return(0);
}
 
int main()
{
    double a, b, f, amax,bmax,amin,bmin,da,db;
    int i, n;
    cin >> amax >> amin >> bmin >> bmax >> da >> db >> n;
    a = amin;
    b = bmin;
    function();
}

#include <iostream>
#include <cmath>
using namespace std;
 
double function()
{
    double a, b, f, amax, bmax, amin, bmin, da, db;
    int i,n;
    for (i = 0; i <= n; i = i + 1) {
        for (a < amax; a = a + da;) {
            for (b < bmax; b = b + db;) {
                if (b > 0)
                {
                    f = (pow((a + pow(-1, i) * b * i), 2)) / (i * sqrt(i + b * b));
                }
                else
                {
                    f = (pow((a + b * i), 2)) / (i * sqrt(i + b * b));
                }
            }
        }
    }
    return(0);
}
 
int main()
{
    double a, b, f, amax,bmax,amin,bmin,da,db;
    int i, n;
    cin >> amax >> amin >> bmin >> bmax >> da >> db >> n;
    a = amin;
    b = bmin;
    function();
}

//объявляем
double function(double amax, double bmax, double amin, double bmin, double da, double db, int n)
{
//исправьте заголовки 2 последних циклов. Они напоминают эротические дневники
//и уберите объявления переменных объявленных параметрами
 
}

//вызываем
double a, b, f, amax,bmax,amin,bmin,da,db;
int i, n;//i уберите - она не нужна
cin >> amax >> amin >> bmin >> bmax >> da >> db >> n;
 function(amax, bmax, amin, bmin, da, db, n);

#include <iostream>
#include <cmath>
using namespace std;
 
double function(double a,double b,int n)
{
    double f=0;
    double sign=1.;
    for (int i = 0; i <= n; i = i + 1) {
        double t1=i * sqrt(i + b * b);
        if (b > 0)
        {
            double t2=a+sign*b*i;
            sign=-sign;
            f += t2*t2 / t1;
        }
        else
        {
            f += (a+b*i)*(a+b*i) / t1;
        }
 
    }
    return f;
}
int main()
{
    int n;
    cin>>n;
    double da, db,amin,amax,bmin,bmax;
    cin>>da>>db>>amin>>amax>>bmin>>bmax;
    cin.get();
    for (double a=amin;a < amax+da; a += da) 
        for (double b=bmin;b < bmax+db; b += db)
            cout<<a<<'\t'<<b<<'\t'<<function(a,b,n)<<endl;
    cin.get();
    return 0;
}