For the Euler beam element shown in the figure, derive the interpolation functions, NI(ξ), stiffness matrix, k, and nodal force vector, f. Assume that the uniformly distributed load is f(x) = f.
Note that the reference element is defined in the domain, ξ = [1, 1].
