Circle Map

/*
/*
Circle Map (円写像)
theta_{n+1}=theta_n+Omega-K/2pi sin(2pi theta_n)

Omega:=a
K:=b
*/

#include 
#include 

main()
{
  double a=0.25,c1=0.2,c2=0.0; /*Initial condition*/
  double b=0.0;
  
  for(b=0.18; b<=0.25; b+=0.000001)
    {
      c2=fmod(c1+a+b*sin(2.*c1*M_PI), 1.);
      
      /*printf("%f %f\n",c1,c2);*/ /*Figure 1*/
      printf("%f %f\n",b,c2); /*Figure 2 */
      c1=c2;
    }
  return 0;
}

出力結果

Figure 1:

Figure 2: Bifurcation Diagram

Back to C Language