编程模拟自然(一)：如何画一颗静态树

    Dim ph As Graphics
    Dim PaintBoard As Bitmap
    Dim centerX, centerY As Long
    Dim modeSequence() As Integer = {4, 5, 4, 6, 5, 6, 1, 2, 1, 3, 2, 3}
    Dim rnd As New Random
    Public Function Calculator(ByVal x As Double, ByVal y As Double, ByVal mode As Integer, ByVal r As Double, ByVal depth As Integer, ByVal radio As Double) As Double()
        Dim R2 As Double
        r = r * radio
        R2 = Math.Sqrt(3) * (r / 2)
        Dim Position(2) As Double
        Dim tempRnd As Double = 1 ‘* rnd.Next(depth, 20) / 10
        If mode = 1 Then
            x = 0 : y = -r * tempRnd
        ElseIf mode = 2 Then
            x = -R2 * tempRnd : y = (r / 2) * tempRnd
        ElseIf mode = 3 Then
            x = R2 * tempRnd : y = (r / 2) * tempRnd
        ElseIf mode = 4 Then
            x = -R2 * tempRnd : y = -(r / 2) * tempRnd
        ElseIf mode = 5 Then
            x = R2 * tempRnd : y = -(r / 2) * tempRnd
        ElseIf mode = 6 Then
            x = 0 : y = r * tempRnd
        End If
        Position(0) = x
        Position(1) = y
        Position(2) = r
        Return Position
    End Function
    Public Sub PaintLine(ByVal x As Double, ByVal y As Double, ByVal mode As Integer, ByVal r As Double, ByVal depth As Integer, ByVal radio As Double)
        Dim position() As Double
        Dim tempR As Double
        Dim tempWidth As Double
        Dim tempColor As Color
        Dim tempRnd As Integer = rnd.Next(1, 3)
        Dim tempAlpha As Integer
        For j = 0 To 1
            For pff = 0 To tempRnd
                position = Calculator(x, y, modeSequence((mode - 1) * 2 + j), r, depth, radio)
                Dim LineVector As New Vector2D(position(0), position(1))
                LineVector.x *= rnd.NextDouble
                LineVector.y *= rnd.NextDouble
                LineVector.SetMag(IIf(depth <= 1, 20, position(2)))
                position(0) = LineVector.x
                position(1) = LineVector.y
                Dim tempMag As Double = LineVector.Magnitude
                Dim PointList As New List(Of PointF)
                For m = 1 To tempMag
                    LineVector.SetMag(m)
                    PointList.Add(New PointF(x + LineVector.x, y + LineVector.y))
                Next
                If PointList.Count >= 2 Then
                    Dim tempNext As Double = rnd.NextDouble
                    tempAlpha = rnd.Next(20, 100)
                    For p = 0 To PointList.Count - 2
                        tempR = r * 0.015 * depth
                        tempWidth = IIf(depth <= 1, Math.Sin(p / 6) * 10, tempR - p / PointList.Count * (tempR - tempR * radio + r * 0.015 * radio))
                        tempColor = IIf(depth <= 1, Color.FromArgb(tempAlpha, Color.Green), Color.Black)
                        If depth <= 1 Then
                            If tempNext < 0.01 Then
                                tempWidth = Math.Sin(p / 6) * 20
                                tempColor = Color.FromArgb(rnd.Next(200, 255), Color.Red)
                            End If
                        End If
                        ph.DrawLine(New Pen(tempColor, tempWidth), PointList(p), PointList(p + 1))
                    Next
                End If
                PaintHexgon(position(0) + x, position(1) + y, modeSequence((mode - 1) * 2 + j), tempMag, depth - 1, radio)
            Next
        Next
    End Sub
    Public Sub PaintHexgon(ByVal x As Double, ByVal y As Double, ByVal mode As Integer, ByVal r As Double, ByVal depth As Integer, ByVal radio As Double)
        If depth > 0 Then
            If mode = 0 Then
                ph.DrawLine(New Pen(Brushes.Black, r * 0.02 * depth), New PointF(x, y), New PointF(x, y + r))
                PaintLine(x, y, 1, r, depth, radio)
            Else
                PaintLine(x, y, mode, r, depth, radio)
            End If
        End If
    End Sub
    ‘画树
    Private Sub Button2_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles Button2.Click
        PaintBoard = New Bitmap(PictureBox1.Width, PictureBox1.Height)
        ph = Graphics.FromImage(PaintBoard)
        ph.SmoothingMode = Drawing2D.SmoothingMode.AntiAlias
        centerX = PictureBox1.Width / 2 : centerY = PictureBox1.Height / 2
        PaintHexgon(centerX, centerY, 0, 200, 5, 0.6180339887)
        ‘ PaintHexgon(centerX, centerY, 0, 20, 5, 1.61812)
        PictureBox1.Image = PaintBoard
    End Sub