Inspiration for implementing the Gaussian process methods came from
the work of Carl Rasmussen and Chris Williams.  I also benefited from
looking at Carl Rasmussen's implementation.

The sampling method for top-level hyperparameters using adaptive
rejection sampling, and the features allowing cpu time to be stored,
displayed, and used to control the number of iterations, are also
based on work by Carl Edward Rasmussen.

Thanks also to Tim Howells, Dirk Husmeier, Marcus Lee, David MacKay,
Tony Plate, Jeff Rosenthal, and Chris Williams for comments on the
documentation, for testing the programs, and for finding bugs.