fit2tsmodel_ucminf function

Numerical optimization of the 2ts model