FORTRAN code for minimizing a function whose uation is expensive. At each iteration, a Bayesian posterior mean for the surface shape conditional on points already sampled is constructed and the minimum of this is found. This minimum is then used as