梯形法求解常微分方程（c++）

#include<iostream>

#include<iomanip>

using namespace std;

int main()

{

      

      double x,y,yn,h,temp,f;     

      x=0;      //对x赋初值

      y=1;      //对y