微分方程数值解法：PCL

微分方程是数学中常见的一类方程，描述了变量之间的变化率。求解微分方程在许多科学和工程领域都是至关重要的任务。然而，许多微分方程往往没有解析解，因此需要使用数值方法来求解。本文将介绍一种常用的微分方程数值解法——PCL（Predictor-Corrector-Loop）方法，并提供相应的源代码。

PCL方法是一种迭代方法，通过不断迭代逼近微分方程的数值解。该方法包含以下步骤：

1. 预测（Predictor）步骤：根据当前的解近似值，使用数值积分方法（如欧拉法或龙格-库塔法）进行预测。假设我们要求解的微分方程为dy/dx = f(x, y)，当前的解近似值为y_i，步长为h，则预测步骤的公式为：
   y_pred = y_i + h * f(x_i, y_i)

2. 修正（Corrector）步骤：使用预测步骤得到的近似值作为初始值，再次应用数值积分方法进行修正。修正步骤的公式为：
   y_i+1 = y_i + h/2 * [f(x_i, y_i) + f(x_i+1, y_pred)]

3. 循环（Loop）步骤：重复进行预测和修正步骤，直到达到所需的精度或满足停止准则为止。停止准则可以是达到指定的迭代次数或两次迭代之