get.ab function

Hyperparameter estimation for a and b.