无法绘制pdf

Question

我无法绘制pdf，因为包含保存概率的值的矩阵无法正常运行。以下是我写的代码，请告诉我哪里出错了。

mu=c(1,-1)
covar=matrix(c(1,-1.6,-1.6,4),2,2,byrow=TRUE)
pdff= function(z){
pdf= ((exp(-t(z-mu)%*% solve(covar) %*% (z-mu))/2)/(6.283* sqrt(det(covar))))
return(pdf)
}


x=seq(3,5,length=50)
y=seq(2,4,length=50)
prob=matrix(NA, 50, 50)
for(i in x){
  for(j in y){

    prob[i,j]= as.vector(pdff(c(i,j)))
  }
}

这里的概率矩阵没有充满概率。有很多值仍然是NA。

我需要一个有助于找到填充了所有值的prob矩阵的解决方案。