Вычислить значение функции на интервале

#include<iostream>
#include<cmath>
using namespace std;
int main()
{
    const double pi = 3.14159265;
    double a = 1;
    for (double x = -9; x <= 21; x+=3)
        if (x != 0)
        {
            double f = (a*x*x + 3) / sin(pi*x / 6);
            cout << x << " " << f << "\n";
        }
    return(0);
}

#include <iostream>
#define _USE_MATH_DEFINES
#include <math.h>
 
/******************************************
Примечания:
 
* Подразумевается что PI*x/6 - значение в радианах
* В данном конкретном случае достаточно использовать 
целочисленные типы переменных, что исключает необходимость 
учитывать погрешность вычислений.
 
******************************************/
 
int main(void)
{
    int a = 0;
 
    std::cout << "Input a: \n";
    std::cin >> a;
 
    std::cout << "F(x) = (" << a << "*x^2+3)/(sin(PI*x/6)); [-9;21] \n\n";
 
    for(int x=-9; x <=21; x+=3)
    {
        std::cout << "F(" << x << ") = ";
        int temp_zn = sin( ( M_PI*x/6.0 ) );
        if(!temp_zn) 
        {
            std::cout << "Impossible \n";
        } else { 
            std::cout << (a*x*x+3)/temp_zn << "\n";
        }
    }
 
    
    system("pause");
    return(0);
}

#include<cmath>

#include <math.h>