def NR_scalar(p):
    """Not vectorizable.
    pow is for scalars only, and if statement"""
    if p > 0:
        return M*pow(p,N)/(pow(s,N)+pow(p,N))
    else:
        return 0
dataset/ETHPy150Open robclewley/compneuro/NR_functions.py/NR_scalar/original