⁴This code was written by Herbert Melenk.