Error in solve.default(H) :
  Lapack routine dgesv: system is exactly singular: U[1,1] = 0
> H
     [,1] [,2]
[1,]    0    0
[2,]    0    0
我觉得可能是H值得问题，H不能够全是0
然后，我有看下
> f1
[1] 423.304
> f2
[1] 423.304
> f3
[1] 423.304
> f4
[1] 423.304
> f5
[1] 423.304
> f6
[1] 423.304

这几个值都相等，符合你的要求吗，如果不符合，可能是你自定义的LL_func函数有问题