def ip(T,x,c):
    """
    Performs intrinsic plasticity
    
    Parameters:
        T: array
            The current thresholds
        x: array
            The state of the network
        c: Bunch
            The parameter bunch
    """
    # T += c.eta_ip*(x.sum()/float(c.N_e)-c.h_ip)
    T += c.eta_ip*(x-c.h_ip)
    return T