keisan.casio.jpに悪魔の階段をUP!

悪魔の階段の図示

です。プログラムはこんな感じで簡単。

for (j=0;j<=100;j=j+1) {
   omega = j/100;
   theta0 = 0.1;
   n=100;
   for (p=1;p<=100;p=p+1) {
      theta0 = theta0+omega-K*sin(2*pi*theta0)/(2*pi);
   }
   theta = theta0;
   for (m=1;m<=n;m=m+1) {
         theta = theta+omega-K*sin(2*pi*theta)/(2*pi);
   }
   W=(theta-theta0)/n;
   println(omega,W);
}

しかし、反復回数とかΩの分割とか、ちょっとでも大きくするとすぐタイムアウトになってしまう...もう少しタイムアウト時間を増やしてもらうともっと複雑な式がUPできるのだけれど。