We experimentally implement an optical algorithm for integration of a real-valued bivariate func- tion. A user-defined function is encoded in the position-dependent phase of one of the polarization ...