Вычисления значений функции на отрезке [a; b] с шагом h - C++

#include <iostream>
#include <math.h>
using namespace std;
#define a 0
#define b 4
#define h 1
void main(){
    float x = a;
    int i = 0;
    float X[int((b-a)/h) + 1], Y[int((b-a)/h) + 1];
    while (x <= b){
        X[i] = x;
        Y[i] = sin(2*x) - cos(2*x);
        x += h;
        ++i;
    }
    for (int j = 0; j < int((b-a)/h) + 1; j++)
        cout << X[j] << '|' << Y[j] << '\n';
}