#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();
}
评论0