Bivariate interpolation with quadratic box splines