I have us being different by a factor of 10^40, but yeah, that's a bit surprising. Maybe we're far enough out in the tails that the normal approximation is breaking down?
Oh, I misbracketed your formula. Yes, 10^40.
I don't offhand have a model for why we expect your method to work, so I don't know why it fails. But another approach using the normal approximation gets within a factor of 10, so that shouldn't be it.
Um, I think you're just counting standard deviations in the wrong direction? You're counting standard deviations from 500,000 and doubling them, but the relevant distribution means are 510,000 and 490,000.
But no, those should be equivalent.
Oh! You're squaring a sum, not summing a square. You're counting the correct...
If it's worth saying, but not worth its own post (even in Discussion), then it goes here.
Notes for future OT posters:
1. Please add the 'open_thread' tag.
2. Check if there is an active Open Thread before posting a new one.
3. Open Threads should be posted in Discussion, and not Main.
4. Open Threads should start on Monday, and end on Sunday.