Моделирование (маятник)

[alexmikt]

подскажите пожалуйста как в данном коде, добавить время движения маятника с начала его движения до его полной остановки:

#include "pch.h"
#pragma hdrstop
#include "graphics.h"


double dt = 0.001;
int dt_display = 16, k_display = 40;

double x[2], v[2];
double g=9.8, v0=1, alpha=0.01, l=5, m=1, c=-0.5;


void init()
{
	x[0] = 5;
	x[1] = 2;
	v[0] = v0*alpha;
	v[1] = v0*alpha;
}

bool tick( double t )
{

	x[0] += v[0]*dt;
	v[0] -=((g*x[0]/l)-(c*v[0]/m))*dt;

	if( fabs(v[0]) < 0.001 && fabs(x[0]) < 0.1 )

		return false;
	return true;
}

void paint()
{
	double k=60;
	text(0,500) << " x=" << setw(10) << x[0] << "     y=" << setw(10) << x[1];
	text(0,480) << "vx=" << setw(10) << v[0] << "    vy=" << setw(10) << v[1];
	text(500,500) << "x1=" << setw(10) << v0*v0*sin(2*alpha)/g;
	particle( x[0]*k+300, x[1]*k );
}