本程序采用改进欧拉法,即梯形法用C++编写。