Computation of the three-loop β-function of the O(N) Gross-Neveu model in minimal subtraction