"krige" <-
function (formula, locations, data, newdata, model = NULL, beta = NULL, 
	nmax = Inf, block = numeric(0), nsim = 0, ...) 
{
    g = gstat(formula = formula, locations = locations, model = model, 
        data = data, beta = beta, nmax = nmax, ...)
    predict.gstat(g, newdata = newdata, block = block, nsim = nsim)
}