We present the definitions, derive the relevant Euler-Lagrange equations, and establish various properties concerning biharmonic maps. We investigate several classes of examples exhibiting singular behavior. Existence of a weak solution to the associated evolution equation is proved using a penalization argument and the Galerkin method. We prove higher integrability of order greater than two for derivatives of Laplacian energy minimizers contingent upon certain energy constraints. We initiate a numerical analysis of biharmonic maps using a discrete Laplacian energy, a finite difference scheme, and involving spherical coordinates in a variety of dimensions in order to understand isolated singularities. Since singularities of Laplacian energy...