shannon<-function(x)
{
an.exp<-function(x,y)
{
  par(mfrow=c(2,2))
b<-boxplot(x)
p<-plot(y)
bar.x<-barplot(x)
bar.y<-barplot(y)
sumario.x<-summary(x)
sumario.y<-summary(y)
par(mfrow=c(1,1))
return(c(sumario.x,sumario.y,b,p,bar.x,bar.y))
}