updated barycentric lagrange interpolation, now works if evaluation points do not have one dimension cooordinate equal to an abscissa