In a scalar parameter estimation case, if the covariance of the observed data is involved with parameter, the distribution function (likelihood function) or the joint distribution of the observation becomes much more complicated. Set the derivative of it to zero and are we guaranteed an unique solution of the parameter? In my case, the derivative is the cubic function of the parameter and there would be three roots. Is it for sure that there would one real and two complex roots?