import numpy as np
def LoG(sigma):
#window size
n = np.ceil(sigma*6)
n = 2
y,x = np.ogrid[-n//2:n//2+1,-n//2:n//2+1]
y_filter = np.exp(-(y*y/(2.*sigma*sigma)))
x_filter = np.exp(-(x*x/(2.*sigma*sigma)))
print(x + y)
final_filter = (-(2*sigma**2) + (x*x + y*y) ) * (x_filter*y_filter) * (1/(2*np.pi*sigma**4))
return final_filter
LoG(1)
aW1wb3J0IG51bXB5IGFzIG5wCgpkZWYgTG9HKHNpZ21hKToKICAgICN3aW5kb3cgc2l6ZSAKICAgIG4gPSBucC5jZWlsKHNpZ21hKjYpCiAgICBuID0gMgogICAgeSx4ID0gbnAub2dyaWRbLW4vLzI6bi8vMisxLC1uLy8yOm4vLzIrMV0KICAgIHlfZmlsdGVyID0gbnAuZXhwKC0oeSp5LygyLipzaWdtYSpzaWdtYSkpKQogICAgeF9maWx0ZXIgPSBucC5leHAoLSh4KngvKDIuKnNpZ21hKnNpZ21hKSkpCiAgICBwcmludCh4ICsgeSkKICAgIGZpbmFsX2ZpbHRlciA9ICgtKDIqc2lnbWEqKjIpICsgKHgqeCArIHkqeSkgKSAqICAoeF9maWx0ZXIqeV9maWx0ZXIpICogKDEvKDIqbnAucGkqc2lnbWEqKjQpKQogICAgcmV0dXJuIGZpbmFsX2ZpbHRlcgogICAgCkxvRygxKQ==