We construct Lagrange interpolating polynomials for a set of points and values belonging to the algebra of real quaternions ℍ ≃ ℝ0,2, or to the real Clifford algebra ℝ0,3. In the quaternionic case, the approach by means of Lagrange polynomials is new, and gives a complete solution of the interpolation problem. In the case of ℝ0,3, such a problem is dealt with here for the first time. Elements of the recent theory of slice regular functions are used. Leaving apart the classical cases ℝ0,0 ≃ ℝ, ℝ0,1 ≃ ℂ and the trivial case ℝ1,0 ≃ ℝ⊕ℝ, the interpolation problem on Clifford algebras ℝp,q with (p,q) ≠ (0,2), (0,3) seems to have some intrinsic difficulties.