Thanks, IRstuff. Your solution is right. However, I am not going to use it after all. The plot should produce an ellipse (joint confidence region of the two parameters), but due to the fact that MC returns one only root from the quadratic equation, the corresponding plot shows one only branch of...