7/10 課題

問題

class Mandel{
  // 実数を引数(パラメタ)に持つメソッド mandel の定義
  static int mandel(double x, double y){
    double zr=0, zi=0,new_zr;
    int i;
    for(i=0;i<50;i=i+1){
      if(zr*zr+zi*zi>4.0) { // |Zi|>2.0 の時は発散と判断
        return i;
      }
       // 複素数の計算
       // (zr+zi*i)*(zr+zi*i)+x+y*i=(zr*zr-zi*zi+x)+(2*zr*zi+y)*i
      new_zr=zr*zr-zi*zi+x;
      zi=2*zr*zi+y;
      zr=new_zr;
    }
    return 50;
  }
  public static void main(String[] args){
    int width=256, height=256;
    System.out.println("P3 "+width+" "+height+" 255");
    int i,j,r;
    double x, y;
    for(i=0;i<height;i=i+1){
      for(j=0;j<width;j=j+1){
           // -1.7 <= x < 0.7, -1.2 <= y < 1.2の範囲で計算
           // 実数(double)から整数(int)への変換は「(double)整数式」
           // のようにおこなう．
        x= -1.7+2.4*((double)j/(double)width);
        y= -1.2+2.4*((double)i/(double)height);
           // 0 <= r <= 50の値が返る
        r=mandel(x,y);
           // r=50の時に黒(0,0,0), r=0の時にほぼ白(250,250,250)とする．
        System.out.println((250-r*5)+" "+(250-r*5)+" "+(250-r*5));
      }
    }
  }
}

1. 上のプログラムを, ~/jousho00 以下に Mandel.java という名前で保存して, コンパイル，実行してみなさい．
2. 上のプログラムは，-1.7 <= x < 0.7, -1.2 <= y < 1.2の範囲でマンデ ルブロー集合を計算しているが, より狭い範囲(元のプログラムよりも縦横そ れぞれ 1/4 以下)を計算するように Mandel.java を書き換えなさい(範囲の選 び方を間違えると，境界部分が含まれないつまらない画像になってしまうので， 注意すること，元のプログラムの画像で左上の点が(-1.7,-1.2), 右下の点が(0.7, 1.2)なので，複雑そうな所を見当をつけてそこを拡大すると良い)．
3. 上のカラー画像は色使いにセンスがないので，Mandel.java

           System.out.println(((50-r)*71)%256+" "+((50-r)*111)%256+" "+((50-r)*97)%256);
   

   の部分を書き換えて(71, 111, 97あたりの数字)，より綺麗な画像を作成するように書き変えなさい．

結果

• 1点1名 元のまま提出
• 2点22名 拡大率不足(縦横それぞれ4倍が目安)，縦横比の変更の結果自 己相似性が崩れている(効果を上げている場合は減点せず)，マンデルブロー集 合の境界が現れない面白くない範囲を指定．
• 3点149名

優秀作品

  System.out.println(((50-r)*0)%256+" "+((50-r)*7)%256+" "+((50-r)*7)%256);

     System.out.println(((100+r)*5)%256+" "+((200+r)*5)%256+" "+(200+r)%256);

        System.out.println(((50-r)*11)%256+" "+((50-r)/2)+" "+((70-r)*2));

        System.out.println((0+r*17%256)+" "+r*23%256+" "+(250-r*37%256));

          System.out.println("255"+" "+((50-r)*43289324)%256+" "+"0");

        System.out.println(((50-r)*2)%256+" "+((50-r)*4)%256+" "+((50-r)*16)%256);

        System.out.println((250+r*3)%256+" "+" 0 " +""+(r*5)%256);

        System.out.println((255-r*51/10)+" "+(255-r*51/10)+" "+255);

        System.out.println(((5-r)*23)%256+" "+((128-r)*23)%256+" "+((39-r)*23)%256);

        System.out.println(((239*r)%256)+" "+((r*r*37)%256)+" "+((123*r*r*r)%256));

System.out.println(((50-r)*0)%256+" "+((50-r)*0)%256+" "+((50-r)*500)%256);  

 System.out.println(((50-r)*255)%256+" "+((50-r)*25)%256+" "+((50-r)*25)%256);

System.out.println(((50-r)*10)%256+" "+((50-r)*20)%256+" "+((50-r)*5)%256);

        System.out.println(((50-r)*3)%256+" "+((50-r)*11)%256+" "+((50-r)*7)%256);

        System.out.println(((50-r)*1)%256+" "+((50-r)*22)%256+" "+((55-r)*241)%256);

        System.out.println((255-r)+" "+(255-r)+" 245");

        System.out.println(((50-r)*22)%256+" "+((50-r)*0)%256+" "+((50-r)*240)%256);

        System.out.println(((50-r)*r+(i+j))%256+" "+((50-r)*111*((j+i)/r))%256+" "+((50-r)*57)%256);

        System.out.println(((50-r)*20)%256+" "+((50-r)*7)%256+" "+((50-r)*10)%256);

         System.out.println(((50-r)*256)%256+" "+((50-r)*0)%256+" "+((50-r)*244)%256);

        System.out.println(((50-r)*0)%256+" "+((50-r)*250)%256+" "+((50-r)*250)%256);

        System.out.println((50-r)*5+" "+(50-r)*5+" "+(220-r));

        System.out.println(((50-r)*6-120)%256+" "+((70-r)*13-80)%256+" "+(67-r*9-30)%256);

        System.out.println((r*2)%256+" "+(r*102)%256+" "+(r*194)%256);