Чтобы максимизировать сумму переменных Ci с учетом ограничения, заданного в вопросе и предполагающего, что они неотрицательны:
library(lpSolve)
cc <- c(0.35, 0.15, 0.15, 0.2, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15)
result <- lp(direction = "max",
objective = rep(1, length(cc)),
const.mat = t(cc),
const.dir = "<=",
const.rhs = 3750 - sum(cc))
result
## Success: the objective function is 24988.33
result$solution
## [1] 0.00 24988.33 0.00 0.00 0.00 0.00 0.00 0.00
## [9] 0.00 0.00