"matern" <- function(u, kappa){ out <- .C("cormatern", as.integer(length(u)), as.double(u), as.double(kappa), res = as.double(rep(0,length(u))))$res return(out) }