R中constrOptim函数的问题

我在R中进行受约束的最优化，用了constrOptim函数，但是没办法计算grad函数（似然函数的导函数？），因此用了Nelder-Mead的方法来最优化，我想问一下R里面怎么用其他的最优化方法来最大化没有导函数的似然函数？求助！紧急！

还有"BFGS", "CG", "L-BFGS-B"，grad函数可以自己写function吧
function(x) { ## Gradient of 'fr'[＃梯度“FR”]
    x1 <- x[1]
    x2 <- x[2]
    c(-400 * x1 * (x2 - x1 * x1) - 2 * (1 - x1),
       200 *      (x2 - x1 * x1))
}