#include <stdio.h>
#include <math.h>

#define OMEGA sqrt(4/((1/3.14159)*(1/3.14159)))
#define A 1

double x(double i)
{
	return A*sin(OMEGA*i);
}

double v(double i)
{
	return A*OMEGA*cos(OMEGA*i);
}

int main()
{
	double i = 0;
	FILE *analyticp, *analyticv;
	analyticp = fopen("analyticp.dat", "w");
	analyticv = fopen("analyticv.dat", "w");

	for(i = 0; i < 10; i += 0.01)
	{
		fprintf(analyticp, "%f\t%f\n", i, x(i));
		fprintf(analyticv, "%f\t%f\n", i, v(i));
	}
	return 0;
}