% Return the initial condition at the given x.
function initC = initialCondition(x)
initC = 400 * 4 * (x - 0.5) * (x - 0.5);