我需要利用c++绘制二次B样条基函数曲线的代码

绘制二次B样条基函数曲线的代码如下：

#include <iostream>
#include <graphics.h>
#include <math.h>

const int MAXN = 1000;

int n; // 控制点数
double u[MAXN], x[MAXN], y[MAXN]; // 控制点

double B(int i, int k, double t) {
    if (k == 0) {
        if (u[i] <= t && t < u[i + 1]) return 1.0;
        else return 0.0;
    }
    double b1 = B(i, k - 1, t), b2 = B(i + 1, k - 1, t);
    double w1 = (t - u[i]) / (u[i + k] - u[i]), w2 = (u[i + k + 1] - t) / (u[i + k + 1] - u[i + 1]);
    return w1 * b1 + w2 * b2;
}

void drawBspline() {
    double t;
    int step = 100; // 绘制精度
    double x1, y1, x2, y2;
    for (int i = 0; i < n - 2; i++) {
        for (int j = 0; j <= step; j++) {
            t = u[i + 1] + (u[i + 2] - u[i + 1]) / step * j;
            x1 = B(i, 2, t) * x[i] + B(i + 1, 2, t) * x[i + 1] + B(i + 2, 2, t) * x[i + 2];
            y1 = B(i, 2, t) * y[i] + B(i + 1, 2, t) * y[i + 1] + B(i + 2, 2, t) * y[i + 2];
            t = u[i + 1] + (u[i + 2] - u[i + 1]) / step * (j + 1);
            x2 = B(i, 2, t) * x[i] + B(i + 1, 2, t) * x[i + 1] + B(i + 2, 2, t) * x[i + 2];
            y2 = B(i, 2, t) * y[i] + B(i + 1, 2, t) * y[i + 1] + B(i + 2, 2, t) * y[i + 2];
            line(int(x1), int(y1), int(x2), int(y2));
        }
    }
}

int main() {
    initgraph(800, 600);

    // 读入控制点
    std::cin >> n;
    for (int i = 0; i < n; i++) {
        std::cin >> u[i] >> x[i] >> y[i];
    }

    // 绘制控制多边形
    for (int i = 0; i < n; i++) {
        circle(int(x[i]), int(y[i]), 3);
    }
    for (int i = 0; i < n - 1; i++) {
        line(int(x[i]), int(y[i]), int(x[i + 1]), int(y[i + 1]));
    }

    // 绘制B样条曲线
    setcolor(RED);
    drawBspline();

    getch();
    closegraph();
    return 0;
}

其中，`n` 表示控制点数，`u`、`x`、`y` 分别为节点向量、控制点的横坐标和纵坐标。`B(i, k, t)` 表示第 `i` 个控制点关于第 `k` 次B样条基函数在参数 `t` 的值。`drawBspline` 函数用于绘制二次B样条基函数曲线。