main=mapM putStr[" Z | 0.0",['0'..'8']>>=(:" 0.0"),"9\n----+",1#69>>"-",do a<-10#30;'\n':shows a" |"++do b<-100#9;" 0."++show(round$5e3+sum[sqrt$exp(-x^2)*8/pi|x<-[2e-4,6e-4..a+b-2e-4]])] a#b=map(/a)[0..b]
Note that non-ascii characters in the above source code will be escaped (such as \x9f).