Sierpinski Triangle - Chaos Game¶

P = [0. -1. 1.; sqrt(3) 0. 0.]
scatter(P[1, :], P[2, :], markersize = 8, ratio=:equal)#, dpi=15)

s = [0.5, 0.]
result = Array{Float64}(2, 0)

for count in 1:1000
    i = rand(1:3)
    @. s = 0.5 * (P[:, i] + s)
    result = hcat(result, s)
end

scatter!(result[1, :], result[2, :],
        markeralpha = 0.3,
        markersize = 2,
        markershape = :utriangle)

#savefig("sier-chaos.png")