函数如下：

fun1=function(data,S0,alpha)
{
  data=as.matrix(data)
  n=nrow(data)
  p=ncol(data)
  X.bar=apply(data,2,mean)
  T1.square=((p*(n-1))/(n-p))*qf(1-alpha,n,n-p)
  x1<-seq(3,6,0.01)
  x2<-seq(35,55,0.01)
  f<-function(x1,x2,x3){t(X.bar-matrix(c(x1,x2,x3),3,1,T))%*%solve(S0)%*%(X.bar-matrix(c(x1,x2,x3),3,1,T))-T1.square/n}
  x3<-outer(x1,x2,f)
  contour(x1,x2,x3)
}
报错： Error in matrix(c(x1, x2, x3), 3, 1, T) :
  argument "x3" is missing, with no default
求问怎么做出三个自变量在一个坐标系内的图形啊，拜托了