用C++解决下面问题：一天文学家要确定一颗小行星绕太阳运行的轨道，他在轨道平面内建立以太阳为原点的直角坐标系，其单位为天文测量单位，在5个不同的时间对小行星作了5次观察，测得轨道上的5个点的坐标数据为(5.764，0.648)，(6.826，1.202)，(6.759，1.823)，(7.168，2.526)，(7.408，3.360)。试确立小行星的轨道方程并画出小行星的运动轨线图形。

题目中给出了小行星轨道上的5个点坐标，我们可以使用最小二乘法拟合出小行星的轨道方程。

假设小行星的轨道方程为 y = a*x^2 + b*x + c，其中a、b、c为待求系数，x、y为点的坐标。

将5个点的坐标代入轨道方程，得到如下方程组：

a*(5.764)^2 + b*(5.764) + c = 0.648
a*(6.826)^2 + b*(6.826) + c = 1.202
a*(6.759)^2 + b*(6.759) + c = 1.823
a*(7.168)^2 + b*(7.168) + c = 2.526
a*(7.408)^2 + b*(7.408) + c = 3.360

将方程组写成矩阵形式，即 Y = XA，其中：

Y = [0.648, 1.202, 1.823, 2.526, 3.360]
X = [[5.764^2, 5.764, 1], [6.826^2, 6.826, 1], [6.759^2, 6.759, 1], [7.168^2, 7.168, 1], [7.408^2, 7.408, 1]]
A = [a, b, c]

使用最小二乘法求解得到 A = [-0.00627, 0.192, -1.276]，因此小行星的轨道方程为 y = -0.00627*x^2 + 0.192*x - 1.276。

最后，我们可以使用绘图库（如matplotlib）绘制小行星的运动轨迹。具体代码如下（假设轨道方程已经求解出来）：

#include <iostream>
#include <cmath>
#include <vector>
#include "matplotlibcpp.h"

namespace plt = matplotlibcpp;

int main() {
    // 生成 x 坐标轴数据
    std::vector<double> x(1000);
    for (int i = 0; i < 1000; i++) {
        x[i] = i / 100.0; // x 范围在 0~10 之间
    }

    // 计算 y 坐标轴数据
    std::vector<double> y(x.size());
    for (int i = 0; i < x.size(); i++) {
        y[i] = -0.00627 * std::pow(x[i], 2) + 0.192 * x[i] - 1.276;
    }

    // 绘制轨迹图形
    plt::plot(x, y);
    plt::show();

    return 0;
}

运行程序后，会弹出一个窗口显示小行星的轨迹图形。