import math
tear=10
angles=(tear-180)/2
def caesaro(t,order,size):
    if order==0:
        t.forward(size)
    else:
        
        for angle in [angles, 180-tear, angles, 0]:
           x=size/(2*(1+math.cos(math.radians(angles))))#some trignometry to make the size of largest square almost constant irrespective of order
          # x=size/2
           caesaro(t, order-1, x)
           t.left(angle)
if __name__=='__main__':
    import turtle
    wn=turtle.Screen()
    nivi=turtle.Turtle()
    nivi.hideturtle()
    nivi.penup()
    
    nivi.pendown()
    nivi.color('blue')
    wn.bgcolor('lightgreen')
    nivi.speed(0)

sz=100
   
    caesaro(nivi,3,sz)

tear1=10
angles1=(tear1-180)/2    
def caesaro1(t,order1,size1):
    if order1==0:
        t.forward(size1)
    else:
        
        for angle in [angles1, 180-tear1, angles1, 0]:
           x=size1/(2*(1+math.cos(math.radians(angles1))))#some trignometry to make the size of largest square almost constant irrespective of order
          # x=size/2
           caesaro1(t, order1-1, x)
           t.left(angle)
if __name__=='__main__':

nivi1=turtle.Turtle()
    nivi1.hideturtle()
    nivi1.penup()
    nivi1.goto(20,0)
    nivi1.left(-90)
    nivi1.pendown()
    nivi1.color('blue')
    nivi1.speed(0)

sz1=100
   
    caesaro1(nivi1,3,sz1)
    
tear2=10
angles2=(tear2-180)/2    
def caesaro2(t,order2,size2):
    if order2==0:
        t.forward(size2)
    else:
        
        for angle in [angles2, 180-tear2, angles2, 0]:
           x=size2/(2*(1+math.cos(math.radians(angles1))))#some trignometry to make the size of largest square almost constant irrespective of order
          # x=size/2
           caesaro2(t, order2-1, x)
           t.left(angle)
if __name__=='__main__':

nivi2=turtle.Turtle()
    nivi2.hideturtle()
    nivi2.penup()
    nivi2.goto(20,-20)
    nivi2.left(-180)
    nivi2.pendown()
    nivi2.color('blue')

nivi2.speed(0)

sz2=100
   
    caesaro2(nivi2,3,sz2)

tear3=10
angles3=(tear3-180)/2    
def caesaro3(t,order3,size3):
    if order3==0:
        t.forward(size3)
    else:
        
        for angle in [angles3, 180-tear3, angles3, 0]:
           x=size3/(2*(1+math.cos(math.radians(angles1))))#some trignometry to make the size of largest square almost constant irrespective of order
          # x=size/2
           caesaro3(t, order3-1, x)
           t.left(angle)
if __name__=='__main__':
    nivi3=turtle.Turtle()
    nivi3.hideturtle()
    nivi3.penup()
    nivi3.goto(0,-20)
    nivi3.left(-270)
    nivi3.pendown()
    nivi3.color('blue')

nivi3.speed(0)

sz3=100
   
    caesaro3(nivi3,3,sz3)

WindMill from Caesaro