################### Representación gráfica del plano  ############
################### No precisa la librería espacial   ############


dibujaplano<-function(x1,x2,plano){
r1<-seq(min(x1),max(x1),length=30)
r2<-seq(min(x2),max(x2),length=30)
r3<-c(plano$coefficients)
fplano<-function(x1,x2,r3){r3[1]+r3[2]*x1+r3[3]*x2}
persp(r1,r2,outer(r1,r2, r3,FUN = fplano),xlab="x1",ylab="x2",zlab="y")
}


###################   En R ################################## 

dibujaplano<-function(x1,x2,plano){
r1<-seq(min(x1),max(x1),length=30)
r2<-seq(min(x2),max(x2),length=30)
r3<-c(plano$coefficients)
fplano<-function(x1,x2,r3){r3[1]+r3[2]*x1+r3[3]*x2}
persp(r1,r2,outer(r1,r2, r3,FUN = fplano),theta=275,phi=25,xlab="x1",ylab="x2",zlab="y")
}