def heron(a):
""" Approximiert die Quadratwurzel von a"""
eps = 0.000000001
x = a / 2.0
while True:
y = (x + a/x) / 2.0
if (x - y) < eps:
break
x = y
return y