import math


def returnNewDelta(delta_old):
    """returns the side of the new polygonal approximation.

    Arguments:
    delta_old -- the side of the previous approximation
    """
    return math.sqrt( 2. * (1. - math.sqrt(1. - 0.25 * delta_old**2) ) )