Summary: Given \(\epsilon > 0\), we compute a function taking prescribed values at \(N\) given points in \(\mathbb{R}^2\), whose \(C^2\)-norm is within a factor \((1 + \epsilon)\) of least possible. The computation takes \(C(\epsilon)N\log N\) computer operations.