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

double function(double t, double k1, double phase)
	{

	double k2=10,m=1;
	double x;

	x= log(t+1)*cos(k1*k2*t+phase)*exp(-t/2);

	return x;
}

main()
	{
	srand(time(NULL));
	double time=.0001, v=1, phase,xvalue=0,yvalue=0.0001, k1;
	double loop=1;
	FILE *out;
	out=fopen("simulation.csv","w");
	double k2=.50,m=1;
	double x;


	while (loop<6)
		{
		phase= rand()*3.14;


		while (yvalue<10)
			{
			xvalue=xvalue+.1*cos(k2*yvalue+phase)*exp(-yvalue/6);

			fprintf(out,"%e,%e\n",xvalue,-yvalue);

			yvalue=yvalue+.001;
		}
		loop=loop+1;
		xvalue=0;
		yvalue=.0001;
	}
}