On the numerical solution of a class of problems in a linear first strain–gradient theory of elasticity