efficient_bfgs function

BFGS Implementation for REML Parameter Estimation