我有3个变量b1，b2和b3，每次都可以有不同的值。我编写了函数f，我希望使用以下约束来最大化f。
不平等约束
b1，b2> = 0且b3 <= 1000
b1 = c(1, 2)
b2 = 3
b3 = c(1, 1)
theta = c(b1, b2, b3)
f = function(theta){
sum(-exp(theta[1:length(b1)]^2
+ theta[(length(b1)+1):(length(b1)+length(b2))]^2)
+ theta[(length(b1)+length(b2)+1):length(theta)]^2)
}
我可以在R中使用哪些包？