working on optimize-theta
construction of the B matrix, beta, the trace_RDR, inversion of the B matrix completed
added get_sigma2.f90
refactored rxy to get_rxy, so that scalar function are prefixed by get
finished optimize_mle