stem <- function(x,scale=1, width=80, atom=0.00000001) {
	if( !is.numeric(x) )
		stop("stem: x must be numeric")
	x <- x[!is.na(x)]
	if(length(x)==0) stop("no non-missing values")
	.C("stemleaf", as.double(x), length(x), as.double(scale), as.integer(width), as.double(atom))
	invisible(NULL)
}
