Fixed a bug in the model
Check on the variance of the estimations of b and epsilon added