A Numerical Treatment of the rf SQUID: I. General Properties and Noise Energy