I personally think the mere idea that you can map on abstract functions that have infinities and boundary conditions, on to a representation that extremely precisely models our physical world, is more astonishing than having a low level understanding of the equations. The final image of the probability density summed over the region around the nucleus, really drove that home for me