The one-dimensional derivation of the Meissner effect is only meant to give a qualitative behaviour, and of course is not "the" solution to the 3D case.
The point is the 3D equation is a PDE, which is therefore complicated, and doesn't really have a general solution. If you assume spherical symmetry, then the solution which doesn't diverge at $+\infty$ is given by $B\sim e^{-r/\lambda}/r$, which roughly gives the same kind of behaviour (the exponent is indeed $r/\lambda$, which is consistent with the OP's reference, but not with the equation written down by the OP as of when this answer was written). The blow up at $r=0$ should perhaps worry us.
But why should we assume spherical symmetry? There's actually no reason why we should. In fact it is often useful to write the equation in cylindrical coordinates, for instance in the study of (London) vortices:
$$\left(\frac{d^2}{dr^2}+\frac{1}{r}\frac{d}{dr}-\frac{1}{\lambda^2} \right) B=0$$
under some appropriate conditions. This is a Bessel-type equation, and can be solved as $B\sim K_0(r/\lambda)$, where $K_0$ is one of the modified Bessel functions. Asymptotically,
$$B \sim \begin{cases}\log{(\lambda/r)} & r \ll \lambda \\e^{-r/\lambda}/\sqrt{r} & r \gg \lambda\end{cases}$$
which gives a similar behaviour (this is a vortex: constant field inside the vortex, logarithmic decrease up to $\approx \lambda$ and then exponential decrease).
It is best to think of the Meissner effect in terms of "expulsion of magnetic fields" rather than in terms of a specific formula.