Построить таблицу значений функции z=sh*(x+y)-sh*(x-2y)/sh*sh*(x-y)

for(double x=1.;x<=2.;x+=0.2)
{
    for(double y=2.;y<=3.;y+=0.1)
    {
       //z=sh*(x+y)-sh*(x-2y)/sh*sh*(x-y); странная очень формула
       double z=sinh(x+y)-sinh(x-2.*y)/sinh(x-y); // не это ли имеется ввиду?
       printf("%lf %lf %lf\n",x,y,z);
    }
}

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
int main()
{
    double xn, xk, hx, yn, yk, hy, x, y, z;//Áåç ïîñëåäíèõ òð¸õ ìîæíî áûëî áû îáîéòèñü
    cout<<"x = [xn,xn + hx...xk]\n";
    cout<<"Enter xn = ";cin>>xn;
    cout<<"Enter xk = ";cin>>xk;
    cout<<"Enter hx = ";cin>>hx;
 
    cout<<"y = [yn,yn + hy...yk]\n";
    cout<<"Enter yn = ";cin>>yn;
    cout<<"Enter yk = ";cin>>yk;
    cout<<"Enter hy = ";cin>>hy;
 
    cout<<setw(3)<<"x"<<setw(3)<<"|"
        <<setw(3)<<"y"<<setw(3)<<"|"
        <<setw(3)<<"z"<<setw(3)<<"|\n";
    for(x = xn; x <= xk; x += hx)
    {
        for(y = yn; y <= yk; y += hy)
        {
            z = sinh(x+y) - sinh(x-2*y)/sinh(x-y)*sinh(x-y);
            cout<<setw(3)<<x<<setw(3)<<"|"
                <<setw(3)<<y<<setw(3)<<"|"
                <<setw(3)<<x<<setw(3)<<"|\n";
        }
    }
    system("pause");
    return 0;
}