oula数值积分.cpp

来自「欧拉法」· C++ 代码 · 共 35 行

#include <iostream.h>
class OULA
{
  private:
  double h,b,u,a;
  public:
  void seth(double l=0){h=l;} //设步长
  void setf(double xa=0,double xb=0,double y=0) //设初值和范围（xa,xb）
{
  b=xb;
  a=xa;
  u=y;
}
 double f(double t,double u) //函数值，修改它以适应各自需要
{
//函数设定
  double f=u-2*t/u; 
  return f;
}
  void oula() //R-K 主函数
{
   for(int count=0;count<(b-a)/h;count++)
   {
    u=u+h*f(a+count*h,u);
    cout<<u<<endl;
   }
  }
}; 
void main()
{
  OULA my;
  my.seth(0.1);
  my.setf(0,1,1);
  my.oula();
}