遠心力の微分方程式の数値解　前進・後退オイラー法

circle.cpp

#include <stdio.h>

int main() {
	double xn = 1.0;
	double yn = 0.0;
	double t = 0.0;
	double dt = 0.01;
	double k = 1.0;
	double x = xn;
	double y = yn;
	for (int i = 0; i < 5000; ++i) {
		//x = xn - yn * dt;
		//y = yn + xn * dt;
		x = (xn - yn * dt) / (1.0 + dt * dt);
		y = (yn + xn * dt) / (1.0 + dt * dt);
		xn = x;
		yn = y;
		printf("%lf %lf %lf\n", t, x, y);
		t = t + dt;
	}
	return 0;
}