Matlab code for okada85:
Copyright (c) 1997-2009, Francois Beauducel