double T = 0.1; double XE[NN]; double YN[NN] = { 0 }; XE[0]=0; for (i = 1; i < NN; i++) { XE[i] = XE[i - 1] + (DVe[i - 1] + DVe[i]) * 0.05; YN[i] = YN[i - 1] + (DVn[i - 1] + DVn[i]) * 0.05; }写成子函数形式 并且在主函数中调用

好的，那么你可以将这段代码封装在一个函数中，函数的参数为DVe、DVn、XE、YN、NN和T。具体代码可以像下面这样：

void calculateXY(double DVe[], double DVn[], double XE[], double YN[], int NN, double T)
{
    int i;
    XE[0] = 0;
    for (i = 1; i < NN; i++)
    {
        XE[i] = XE[i - 1] + (DVe[i - 1] + DVe[i]) * T / 2;
        YN[i] = YN[i - 1] + (DVn[i - 1] + DVn[i]) * T / 2;
    }
}

然后在主函数中调用这个函数即可：

int main()
{
    double DVe[200], DVn[200], XE[NN], YN[NN];
    int NN = 200;
    double T = 0.1;
    // 读取数据到DVe和DVn数组中
    readDataFromFile("shuju3.txt", data);
    // 计算XE和YN数组中的值
    calculateXY(DVe, DVn, XE, YN, NN, T);
    // 其他操作
    return 0;
}