Find *your* computer's machine epsilon (smallest available precision) for both single and double precision floating point variables so you can use it as a stopping criterion in future iterative computations. To do this, you should write 2 procedures (functions or subroutines) to return the machine epsilon, one for single precision variables, and one for double precision variables. You can refer to the pseudo-code presented in class as a starting point:

epsilon = 1
DO
IF(epsilon + 1 <= 1) EXIT
epsilon = epsilon/2
LOOP
epsilon = 2 * epsilon