Gibbs sampler and "plotErr2" modified to 2D-benchmarks

Gibbs sampler works in 2D, but it often can't find any points with
nonzero probability of improvement -- has to be corrected.
gpedaStep2d() and plotErr2() plots 3D plots of surfaces.
Minor changes in optimiztation_test.m