Как рисовать графики в С++?

#include <iostream>
#include<cmath>
using namespace std;
 
int main()
{
    int G = 73000, Iz = 660000, TF = 20;
    double Xt = 0.24, V = 97.2, p = 0.1190, g = 9.81, S = 201.45, ba = 5.285,
        cy0 = -0.255, cya = 5.78, cy_dv = 0.2865, cx = 0.046, mz0 = 0.2, mzw_ = -13.0,
        mza_ = -3.8,
        mza = -1.83, mz_dv = -0.96, DD = 0.5, DT = 0.01, TD = 0;
    double C1, C2, C3, C4, C5, C9, C16, DV_bal, a_bal, cy_bal, m, T = 0, Vpr;
    double XV = -17.86, Kw0 = 0.112, DVD, KWZ = 1, TWZ = 0.7, DV, DVS, NY;
    
    m = G / g;
    C1 = (-mzw_ / Iz) * S * pow(ba, 2) * (p * V / 2);
    C2 = (-mza / Iz) * S * ba * (p * pow(V, 2) / 2);
    C3 = (-mz_dv / Iz) * S * ba * (p * pow(V, 2) / 2);
    C4 = ((cya + cx) / m) * S * (p * V / 2);
    C5 = (-mza_ / Iz) * S * pow(ba, 2) * (p * V / 2);
    C9 = (cy_dv / m) * S * (p * V / 2);
    C16 = V / (57.3 * g);
    cy_bal = 2 * G / (S * p * pow(V, 2));
    a_bal = 57.3 * ((cy_bal - cy0) / cya);
    DV_bal = (-57.3 * (mz0 + mza * a_bal / 57.3 + cy_bal * (Xt - 0.24))) / mz_dv;
    
    cout << "C1 = " << C1 << endl;
    cout << "C2 = " << C2 << endl;
    cout << "C3 = " << C3 << endl;
    cout << "C4 = " << C5 << endl;
    cout << "C9 = " << C9 << endl;
    cout << "C16 = " << C16 << endl;
    cout << "cy_bal = " << cy_bal << endl;
    cout << "a_bal = " << a_bal << endl;
    cout << "DV_bal = " << DV_bal << endl;
 
    double* X = new double [5];
    double* Y = new double [5];
 
    int max = 2 * TF + 1;
    
    double o = TF / DD;
    
    // Массивы значений в таблице
    double* Time = new double[max];
    double* massNY = new double[max];
    double* massY1 = new double[max];
    double* massY3 = new double[max];
    double* massY4 = new double[max];
    double* massDV = new double[max];
    double* massXV = new double[max];
 
    int count = 0;
 
    while (T <= TF) {
        DVS = Kw0 * XV;
        DVD = 0;
        DV = DVS + DVD;
        X[0] = Y[1];
        X[1] = -C1 * Y[1] - C2 * Y[3] - C5 * X[3] - C3 * DV;
        X[2] = C4 * Y[3] + C9 * DV;
        X[3] = X[0] - X[2];
        X[4] = DVD;
        
 
        NY = C16 * X[2];
        
        for (int i = 0; i < 5; i++) {
            Y[i] = Y[i] + X[i] * DT;
        }
        if (T >= TD) {
            massNY[count] = NY;
            massY1[count] = Y[0];
            massY4[count] = Y[3];
            massY3[count] = Y[2];
            Time[count] = T;
            massXV[count] = XV;
            massDV[count] = DV;
            TD += DD;
            count++;
        }
        T += DT;
    }