is.diffusion <- function(x){ return(inherits(x,"diffusion")) }
is.diffusion <- function(x){
return(inherits(x,"diffusion"))
}