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

int main(int argc,char **argv)
{
  int i,n,nd;
  double theta,omega,domega,l,g;
  double dt,end,iv,lth;
  int sw=0;

  l=3.123; // シミュレーション対象パラメータ
  g=9.80;
  theta=90.0/180.0*M_PI; // θの初期値θ0

  if(argc>1) 
    {
      int id=atoi(argv[1]);
      l=3+(id/100)%10;
      theta=(id%10+80)/180.0*M_PI;
      printf("id=%d,  l=%f,  theta=%f(%f)\n",id,l,theta,theta/M_PI*180);
    }

  // /180*M_PIで度→ラジアン
  omega=0.0; // ωの初期値

  dt=0.005;   // 計算区分
  end=10;       // 終了時間
  iv=0.1;      // 表示間隔
 

  for(i=2;i<argc;i++)
    {
      if(strncmp(argv[i],"dt=",3)==0) { dt =atof(argv[i]+3); }
      if(strncmp(argv[i],"en=",3)==0) { end=atof(argv[i]+3); }
      if(strncmp(argv[i],"iv=",3)==0) { iv =atof(argv[i]+3); }
      if(strncmp(argv[i],"th=",3)==0) { theta=(atof(argv[i]+3))/180.0*M_PI; }
      if(strncmp(argv[i],"l=",2)==0) { l=(atof(argv[i]+2)); }
      if(strncmp(argv[i],"sw",2)==0) { sw=1; }
    }

  n=(int)(end/dt);  // 計算繰返数
  // "(int)"は「整数にする(切捨)」という指示。
  nd=(int)(iv/dt+0.5);   // 表示周期

  // printf("%d %d \n",n,nd);

  for(i=0;i<n;i++)  // 繰り返し
    {
      domega=(-g/l)*sin(theta);

      if(sw)
	{
	  omega=omega+domega*dt;
	  theta=theta+omega*dt;
	}
      else
	{
	  theta=theta+omega*dt;
	  omega=omega+domega*dt;
	}
      if((lth*theta<0)||(i%nd==0))    // nd回に１回表示
	{
	  if((i<n/20)||(i>n*19/20)||(lth*theta<0))
	  
	  printf("time:%5.2lf[s]  θ:%6.3lf[rad](%7.3lf[deg])  ω:%6.3lf[rad/s](%8.2lf[deg/s])\n",
		 i*dt,theta,theta/M_PI*180,omega,omega/M_PI*180);
	}
      lth=theta;
    }
  return 0;
}