An adaptive numerical integration routine for the IBM 1620 II