Create BLR_fit.R