Star Sky Night

import turtle
import math

t = turtle.Turtle()
t.speed(0)

s = 1
for i in range(5000):
    t.color((255 * 10 ** (-i * 0.002), 255 * 10 ** (-i * 0.05), 255 * 10 ** (-i * 0.002)))
    print(s, i)
    t.forward(s / 2)
    t.left(90)
    t.forward(1)
    s = math.sqrt(s ** 2 + s)
    t.forward(s / 2)
    t.left(90 + math.sqrt(i))