Convex envelope of bivariate cubic functions over rectangular regions