Abstract We discuss a pointwise numerical differentiation formula on multivariate scattered data, based the coefficients of local polynomial interpolation at Discrete Leja Points, written in Taylor’s monomial basis. Error bounds for approximation partial derivatives any order compatible with function regularity are provided, as well sensitivity estimates to functional perturbations, terms inver...