differential.cpp

常微分方程数值解法

#include"BetterEuler.h"
#include"Adamsin.h"
#include"RungeKutta.h"
double f(double t,double u){
	return (-5)*u;
}
double u(double t){
	return exp(-5*t);
}
double f1(double t,double u){
	return exp(10*(t-u));
}
double u1(double t){
	return log(exp(10*t)+exp(1)-1)/10;
}
void main(){
    //BetterEuler<double> c(f1);
	//c.BetterEulermeth();
	Adamsin<double> a(u1,f1);
	a.Calculate();
	a.Error();
	RungeKutta<double> b(u1,f1);
	b.Calculate();
	b.Error();
}