Past exam of the year 2 of the earth sciences (geology) course of the University of Oxford / 2024 / Paper 3 / 5 / a

Put $h=\Delta x$ and suppose $L=(n+1)h$, so there are $n$ interior grid points $x_j=jh$, $j=1,\dots ,n$. Write $u_j(t)\simeq u(x_j,t)$. The fixed ends impose homogeneous Dirichlet boundary conditions,For a sufficiently smooth displacement field, the two Taylor series at fixed time givewhere all displayed derivatives on the right are evaluated at $(x_j,t)$. Adding cancels the odd terms, leaving the second-order central differenceSubstitute this finite difference into the wave equation and retain continuous time. The resulting method of lines gives coupled second-order ordinary differential equations,In particular, the first and last equations are ${\ddot{u}}_1=(c^2/h^2)(-2u_1+u_2)$ and ${\ddot{u}}_n=(c^2/h^2)(u_{n-1}-2u_n)$ for $n\ge 2$. With the interior displacement vector only,Equivalently, $A_{jj}=-2c^2/h^2$, $A_{j,j+1}=A_{j+1,j}=c^2/h^2$, and all other entries vanish. This entrywise definition also covers small $n$, including $A=(-2c^2/h^2)$ when $n=1$. The matrix is a symmetric matrix and a tridiagonal matrix. The spatial approximation has error $O(h^2)$ in Big O notation; time has not yet been discretized.