put return back in algorithm 3.1

more clear now that gamma_k only is computed once, and not for every j,
for consistency and use in next algorithm there should be a return