在RStudio下，代码如下：
x0<-rep(0,23)
fn1<-function(x){
x1<-as.matrix(x)
t(x1)%*%D%*%x1-sum(mm*x1) ###D是23*23矩阵，mm是长度为23的向量
}
res <-solnp(x0,fun = fn1,eqfun = sum(x),eqB = 1)
Error in fun(pars, ...) : 没有"fun"这个函数
怎么试都不行，solnp函数明明有这个fun啊
http://huangfujie1987.blog.163.com/blog/static/212049011201531351014534/