Nonlinear Aspects of the BOLD Response in Functional MRI