Mandelbrotov set

#Mandelbrotov set crtanje
from turtle import*
print('crtanje je kvadratno ')
#pocetak_x=float(input('Unesi početak intervala za realni dio '))
#kraj_x=float(input('Unesi kraj intervala za realni dio '))
#pocetak_y=float(input('Unesi početak intervala za imaginarni dio '))
#kraj_y=float(input('Unesi kraj intervala za imaginarni dio '))

dimenzija=float(input('Unesi dimenziju kvadratnog crtanja seta '))
gustoca=int(input('Unesi gustocu crtanja '))
provjera=int(input('Unesi broj provjera za svaki broj '))
def funkcija(x,y):
    realni=x
    imaginarni=y
    c_realni=x
    c_imaginarni=y
    broj_provjera=1
    rezultat=(realni**2 + imaginarni**2)**(1/2)
    da_il_ne=0
    if rezultat<2:
        while broj_provjera<provjera and rezultat<2:
            privremena=realni
            privremena2=imaginarni
            realni=privremena**2 - privremena2**2 + c_realni
            imaginarni= 2*privremena*privremena2 + c_imaginarni
            rezultat = ((realni)**2 + (imaginarni)**2)**(1/2)
            broj_provjera= broj_provjera+1
    else:
        da_il_ne=0
    if rezultat<2:
        da_il_ne=1
    else:
        da_il_ne=0
    c=[broj_provjera, da_il_ne]
    return(c)
#neke provjere

#print('provjera za 0,i ',funkcija(0,1)[1])
#print('provjera za 0,0.3i ',funkcija(0,0.3)[1])


#crtanje
#ht()
lt(90)
pu()
fd (100)
lt(90)
fd(100)
rt(180)
pd()

#razmak_x = (kraj_x-pocetak_x)/gustoca
#razmak_y = (kraj_y-pocetak_y)/gustoca
#neki_y=pocetak_y
#gustoca_y=(kraj_y-pocetak_y)*200
#gustoca_x=2

razmak_x = 2*dimenzija/gustoca
razmak_y= 2*dimenzija/gustoca
neki_x=-dimenzija
neki_y=dimenzija

for i in range (2*gustoca):
    neki_x = -dimenzija
    for j in range(2*gustoca):
        if funkcija(neki_x,neki_y)[1]==1:
            dot(2)
        neki_x=neki_x+razmak_x
        pu()
        fd(1)
        pd()
    neki_y=neki_y-razmak_y
    pu()
    back(gustoca*2)
    rt(90)
    fd(1)
    lt(90)
    pd()