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

int main(void)
{
	int i,n,nd;
	double l,g,theta,omega,domega;
	double dt,end,iv;

	theta=88.0/180.0*M_PI;
	omega=0;
	l=5.0;
	g=9.80;
	
	dt=0.001;
	end=10;
	iv=0.1;

	n=(int)(end/dt);

	nd=(int)(iv/dt);

	for(i=0;i<n;i++)
	{
		domega=-g/l*sin(theta);
		theta=theta+omega*dt;
		omega=omega+domega*dt;

		if(i%nd==0)
		{
			printf("time:%6.2f[s] θ:%6.3f[rad] ω:%6.3f[rad/s]\n",i*dt,theta,omega);
		}
	}
	return 0;
}