% Run the solver u = poisson1d(f, nx);
A stepped bar fixed at both ends, with axial load at mid-node.