« 鳩@豊中駅 (IXY 20IS) | トップページ | カッシーニの卵形線をkeisan.casio.jpにUP! »

2009年4月20日 (月)

マンデルブロー集合をExcelで描いてみる(ソース付き)。

一周年記念ソース公開第三弾。Excel VBAでマンデルブロー集合を描いてみる。準備として、セルを小さくしてもらって。。。(セル1つがドット1つのイメージ。設定がめんどくさいんで手動で。。。)。こんな感じで掛ける。色は適当に変えてもらったら。

Mandel

「mandel.txt」をダウンロード

中身はこんな感じ。

Private Sub CommandButton1_Click()
    Dim a As Double, b As Double
    Dim i As Integer, j As Integer, k As Integer
    Dim amax As Double, amin As Double
    Dim bmax As Double, bmin As Double
    Dim x As Double, y As Double, xdummy As Double, ydummy As Double
    Dim n As Integer, kmax As Integer
   
    Application.ScreenUpdating = False
   
    amax = 1.2
    amin = -2
    bmax = 1.2
    bmin = -1.2
    n = 200
    kmax = 100
   
    For i = 0 To n
        For j = 0 To n
            a = amin + (amax - amin) * i / n
            b = bmin + (bmax - bmin) * j / n
            x = 0
            y = 0
            For k = 1 To kmax
                xdummy = x
                ydummy = y
                x = x * x - y * y + a
                y = 2 * xdummy * ydummy + b
                If Sqr(x * x + y * y) > 2 Then
                    col = k
                    Exit For
                End If
                col = 0
            Next k
            Worksheets("Sheet1").Range(Cells(j + 2, i + 2), Cells(j + 2, i + 2)).Interior.Color = RGB((col * 20) Mod 255, (col * 20) Mod 255, (col * 20) Mod 255)
        Next j
    Next i
   
    Application.ScreenUpdating = True
   
End Sub

« 鳩@豊中駅 (IXY 20IS) | トップページ | カッシーニの卵形線をkeisan.casio.jpにUP! »

学問・資格」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

トラックバック

この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/44722889

この記事へのトラックバック一覧です: マンデルブロー集合をExcelで描いてみる(ソース付き)。:

« 鳩@豊中駅 (IXY 20IS) | トップページ | カッシーニの卵形線をkeisan.casio.jpにUP! »

最近のコメント

2018年11月
        1 2 3
4 5 6 7 8 9 10
11 12 13 14 15 16 17
18 19 20 21 22 23 24
25 26 27 28 29 30  
フォト
無料ブログはココログ