/*
* File: euler-cromer.cpp
* Author: Muhammad Faruq Nuruddinsyah
*
* Created on October 27, 2015, 9:26 PM
*/


#include <iostream>
#include <fstream>

using namespace std;

int main(){
	float xn, xn_1, vn, vn_1;
	float k = 1.0, m = 1.0, dt = 0.001;
	
	ofstream file;
	file.open("euler-cromer.dat");
	
	xn_1 = 1.0;
	vn_1 = 0.0;
	xn = 1.0;
	vn = 0.0;
	
	for(int i = 0; i < 10000; i++){
		xn = xn_1 + vn * dt;
		vn = vn_1 - ((k * xn) * dt) / m;
		
		file << xn << "\t" << vn << endl;
		
		xn_1 = xn;
		vn_1 = vn;
	}
	
	file.close();
	return 0;
}