Normal distribution curve and random number in LUA

So say I have a range of numbers 1-1000. math.random(1,1000) would give me an equal chance of getting each number. Instead of this, I want to make a distribution curve so that the chance of getting 1 equals the chance of getting 1000 but the chance of getting 500 for example is much more common. How would I go about making this?