¹This code was written by Herbert Melenk and Francis J. Wright.