Learning BOLD Response in fMRI by Reservoir Computing