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

int main(void)
{
	double base_freq = 440.0;
	double low=-14, high=36, step=1;
	int N=(high - low) / step + 0.5;
	double pitch[N];
	double freq[N];
	double mass[N];
	double pos[N];
	double vel[N];
	double amp[N];
	double phase[N];
	int i;
	int sample_rate = 44100;
	double dt = 1.0 / sample_rate;
	double t=0;
	double pi=acos(-1);
	double damp=0.000001;	// amplitude should be damped to 1/1000000th after 1 second
	double dampf=pow(damp, 2*dt);	// I don't understand why it needs to be 2*
	double sample;
	double outdt = 0.05;
	double outt = 0;
	int focus = 0;
	focus = (focus-low)/step+0.5;

	for (i=0; i<N; ++i) {
		pitch[i] = low + i * step;
		freq[i] = base_freq * pow(2, pitch[i]/12.0);
		mass[i] = 1.0/(4.0*pi*pi*freq[i]*freq[i]);
		pos[i] = 0.0;
		vel[i] = 0.0;
	}

	while(1) {
		if (scanf("%lf", &sample) <= 0)
			break;
		for (i=0; i<N; ++i) {
			vel[i] += (sample-pos[i])/mass[i]*dt;
			vel[i] *= dampf;
			pos[i] += vel[i]*dt;
			amp[i] = sqrt(pos[i]*pos[i] + mass[i]*vel[i]*vel[i]);
//			phase[i] = atan2(pos[i], vel[i]*sqrt(mass[i]))/pi/2.0*360.0;
		}
//		printf("%f %f %f %f %f\n", t, pos[focus], vel[focus], amp[focus], phase[focus]);
		t+=dt;
		if (t > outt) {
			for (i=0; i<N; ++i) {
				printf("%f %f %f\n", t, pitch[i], amp[i]);
			}
			outt += outdt;
		}
	}
	exit(0);
}