目标函数是copula函数密度、两个边缘分布密度函数的乘积，三个密度函数的表达式都比较复杂,求其最大值。。求指点，用什么方法做比较好？代码如下，fn是目标函数，heq是等式约束。。u、v是边缘分布概率，parq、parw是边缘分布参数。。

自己试验过一般的非线性规划求极值，包括alabama::auglag函数、Rsolnp::solnp函数。但效果不理想。。