Dear SciPy developers, as requested I created an enhancement ticket (http://projects.scipy.org/scipy/ticket/1352). I attached FORTRAN 95, FORTRAN 77 and ANSI C versions of the subroutine, so you may pick your favorite language :) Just let me know what the next step (if any) should be. Merry Christmas! Regards, Pim Schellart