¹This code was written by Herbert Melenk.