Numerical computation of the optimal vector field in a fishery model