I think a really nice adition to Julia (either in Base or as a separate package) would be the possibility of adaptive sampling. I am thinking in terms of a stand-alone routine for n-dimensions, which could then be called for plotting, numerical integration, etc.
I am currently doing research in theoretical physics and in my day to day work I often encounter functions which are very expensive to compute and plotting them would go so much faster with the help of a good sampling scheme.
As a starting point, see this question on the Mathematica(R) stack exchange community:
http://mathematica.stackexchange.com/questions/216/adaptive-sampling-for-slow-to-compute-functions-in-2dI think this should not be hard to implement in Julia, but we could easily do better, for instance with support for parallel evaluations.
So, first I would like to know who would be interested in having such functionality in Julia, and then we decide how to proceed.