using Plots gr(); 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")