We propose efficient methods for the numerical approximation of the field of values of the linear pencil A − λ B, when one of the matrix coefficients A or B is Hermitian and λ ∈ . Our approach builds on the fact that the field of values can be reduced under compressions to the bidimensional case, for which these sets can be exactly determined. The presented algorithms hold for matrices both of small and large size. Furthermore, we investigate spectral inclusion regions for the pencil based on certain fields of values. The results are illustrated by numerical examples. We point out that the given procedures complement the known ones in the literature.