Dutor

//! 作者：Hou Fenglin
//! 程序名：euler.cpp
//! euler预测－校正法解微分方程
#include <iostream>
#include <cmath>
#include <iomanip>
 
using namespace std;
 
class euler
{
    private:
        int n;
        double f, h, x0, x_last, yc, yp;
    public:
        double func(double z, double t)
        {
            f = (1+z) * t * t / 2;
            return f;
        }
        void pc();
};
 
//主函数
int
main()
{
    euler predictor_corrector;
    predictor_corrector.pc();
}
 
//Euler 预测－校正法
void euler::pc()
{
    cout<<"输入初始条件："<<endl;
    cout<<"\n输入 x0 ";
    cin>>x0;
    cout<<"\n输入y0 ";
    cin>>yc;
    cout<<"\n输入y需要的x值 ";
    cin>>x_last;
    cout<<"\n输入等分数";
    cin>>n;
    h = (x_last - x) / n;
    for(int i = 0; i < n; i++)
    {
        yp = yc + h * func(x0, yc);
        yc += 0.5 * h * (func(x0, yc) + func((x0+h), yp));
        x0 += h;
        cout.precision(10);
        cout<<x0<<setw(15)<<yc<<endl;
    }
}