I found this function in 2012 and was amazed by the beauty! It still catches my eye after so many years. Here is the MATLAB code!
Department of Mathematics
University of Oklahoma
601 Elm Avenue
Norman, Oklahoma 73019-3103