def METHOD_NAME(dist, interval_min, interval_max):
    """
    Compute the standard deviation of a distribution within an interval.
    Return the standard deviation of the array elements which are interpreted
    as being taken from a circular interval with endpoints given by
    interval_min and interval_max.
    Parameters
    ----------
    dist : array_like
        Distribution of values within an interval.
    interval_min, interval_max : float
        The endpoints of the interval.
    Returns
    -------
    std : float
        The standard deviation of the distribution.
    """
    # transform distribution from original interval to [-pi, pi]
    half_width = (interval_max - interval_min) / 2.0
    center = interval_min + half_width
    a = (np.asarray(dist) - center) / (half_width) * np.pi
    # compute the angular standard dev. and convert back to original interval
    a_std = angular_std(a)
    return a_std * half_width / np.pi