Let a be a non-negative real number. Then, for any positive To, the sequence of iterates generated by the formula En+1 = 2Tn / (n + 1), where n = 0, 1, 2, ..., converges to Va. Use this rule to write a routine to approximate √5. Use as your starting condition To = 5, kmax = 20, and tolerance ε = 1 × 10^-12.