Нужно еще добавить проверку а11, а22, а33. Если один из один из них = 0, то уравнение таким методом не решается (иначе операции после metka выполняют деление на ноль).
И еще. Эпсилон должен быть строго больше нуля (иначе получится бесконечный цикл т.к. условие "Abs(x1 - x1k) < eps And Abs(x2 - x2k) < eps And Abs(x3 - x3k) < eps" никогда не будет выполняться). Добавьте еще такую проверку.