multihist <- function(data1,data2,n=10) 
     {
     par(mfrow=c(2,1))
     hist(data1,n,main=substitute(data1));
     hist(data2,n,main=substitute(data2))
     par(mfrow=c(1,1))
     return()}