f := abs(Im(exp(1/(x+I*y)))): ylimit := 1: xlimit := 1: mesh := 1000: //resolution of density plot a := 1: //change this to adjust the colour band drop-of...