首先要写出最大似然函数，以及可能的约束条件，然后根据对非线性规划求出参数值。
问题在于极大似然函数出现gamma函数，如何用R编程实现