started to implement multivariate barycentric lagrange interpolation