如题，u是前面步骤产生的矩阵2列，不知道是什么算法问题吗不太懂
L<-function(pa){
  a<-pa
  norm.cop<-normalCopula(a,dim=2)
  l<-1/nT*sum(log(dCopula(u,norm.cop)))
  return(l)
  }
est2<-maxLik(L,start=0.5)

报错：Error in maxNRCompute(fn = function (theta, fnOrig, gradOrig = NULL, hessOrig = NULL,  :
           NA in the initial gradient