In [KP16] (arXiv:1605.07880) the authors introduced a second-order variational problem in L∞. The associated equation, coined the ∞-Bilaplacian, is a \emph{third order} fully nonlinear PDE given by Δ2∞u:=(Δu)3|D(Δu)|2=0. In this work we build a numerical method aimed at quantifying the nature of solutions to this problem which we call ∞-Biharmonic functions. For fixed p we design a mixed finite element scheme for the pre-limiting equation, the p-Bilaplacian Δ2pu:=Δ(|Δu|p−2Δu)=0. We prove convergence of the numerical solution to the weak solution of Δ2pu=0 and show that we are able to pass to the limit p→∞. We perform various tests aimed at understanding the nature of solutions of Δ2∞u and in 1-d we prove convergence of our discretisation to an appropriate weak solution concept of this problem, that of -solutions
