Multistable dynamics mediated by tubuloglomerular feedback in a model of coupled nephrons