We consider the Lidstone–Euler interpolation problem and the associated Lidstone–Euler boundary value problem, in both theoretical and computational aspects. After a theorem of existence and uniqueness of the solution to the Lidstone–Euler boundary value problem, we present a numerical method for solving it. This method uses the extrapolated Bernstein polynomials and produces an approximating convergent polynomial sequence. Particularly, we consider the fourth-order case, arising in various physical models. Finally, we present some numerical examples and we compare the proposed method with a modified decomposition method for a tenth-order problem. The numerical results confirm the theoretical and computational ones.