def compute_max_firing_rate(spiketrains, runtime):
    """
    Compute maximum firing rate
    :param spiketrains: spike trains of neurins
    :param runtime: simulation runtime
    :return: float
    """
    rates = list()
    for spiketrain in spiketrains:
        rates.append(len(spiketrain) * (1000 / runtime))
    return max(rates)