Уравнение Эрмита

#include <iostream>
 
void Hermite(int *H, int n);
void differential(int *polynomial, int n);
void printHermite(int *H, int n);
 
int main()
{
    int *H;
    int n;
 
    std::cout << "Enter n: ";
    std::cin >> n;
 
    H = new int [n + 1];
 
    Hermite(H, n);
 
    std::cout << "H" << n << "(x) = ";
 
    printHermite(H, n);
 
    delete [] H;
 
    std::cin.get();
    return 0;
}
 
void Hermite(int *H, int n)
{
    for (int i = 0; i <= n; i++)
        H[i] = 0;
 
    H[0] = 1;
 
    for (int i = 1; i <= n; i++)
    {
        int *Hn = new int [i];
 
        for (int j = 0; j < i; j++)
            Hn[j] = H[j];
        
        for (int j = i - 1; j >= 0; j--)
            H[j + 1] = H[j];
 
        H[0] = 0;
 
        differential(Hn, i);
 
        for (int j = 0; j < i; j++)
            H[j] -= Hn[j];
 
        delete [] Hn;
    }
}
 
void differential(int *polynomial, int n)
{
    polynomial[0] = 0;
 
    for (int i = 0; i < n - 1; i++)
        polynomial[i] = (i + 1) * polynomial[i + 1];
 
    polynomial[n - 1] = 0;
}
 
void printHermite(int *H, int n)
{
    if (n != 0)
    {
        if (H[n] != 1)
            std::cout << H[n];
        
        std::cout << "x";
 
        if (n > 1)
            std::cout << "^" << n;
    }
    else
        std::cout << "1";
 
    for (int i = n - 1; i > 0; i--)
    {
        if (H[i] == 0)
            continue;
        else
        {
            if (H[i] == 1)
                std::cout << " + ";
            else
            {
                if (H[i] == -1)
                    std::cout << " - ";
                else
                {
                    if (H[i] < 0)
                        std::cout << " - " << -H[i];
                    else
                        std::cout << " + " << H[i];
                }
            }
 
            std::cout << "x";
 
            if (i > 1)
                std::cout << "^" << i;
        }
    }
 
    if (n != 0)
    {
        if (H[0] != 0 && H[0] != 1)
        {
            if (H[0] < 0)
                std::cout << " - " << -H[0];
            else
                std::cout << " + " << H[0];
        }
    }
}