def dble_pwr_law(time_bins, tau, alpha, beta, factor=1):
    """
    Double Power Law star formation history

    Parameters
    ----------
    time_bins: list or numpy.ndarray
        Time bins
    tau: float or int
    alpha: float or int
    beta: float or int
    factor: float or int, optional
        Default = 1

    Returns
    -------
    Array containing the star formation history corresponding to the given time_bins
    """
    return factor/((time_bins/tau)**alpha + (time_bins/tau)**(-beta))