  A numerical scheme for computing arc-length parametrized curves of low
bending energy that are confined to convex domains is devised. The convergence
of the discrete formulations to a continuous model and the unconditional
stability of an iterative scheme are addressed. Numerical simulations confirm
the theoretical results and lead to a classification of observed optimal curves
within spheres.
