# include # include using namespace std; double euler ( double t0, double u0, double dt, double f ( double t, double u ) ); double f1 ( double t, double u ); int main ( ) { double dt, t0, t1, tmax, u0, u1; u0 = 10000; t0 = 2011; tmax = 2020; dt = 1.0; while ( true ) { cout << " " << t0 << " " << u0 << "\n"; if ( tmax <= t0 ) { break; } t1 = t0 + dt; u1 = euler ( t0, u0, dt, f1 ); t0 = t1; u0 = u1; } return 0; } double euler ( double t0, double u0, double dt, double f ( double t, double u ) ) { double u1; u1 = u0 + dt * f ( t0, u0 ); return u1; } double f1 ( double t, double u ) { double dudt; dudt = 0.05 * u; return dudt; }