multivariate barycentric lagrange interpolation now works (only tested in 2d)