I need help understanding how to do this problem through Finite Element Analysis. Using the global stiffness matrix to derive the solution of the problems below.

The beam shown in the figure below is a wide$-$flange W$16$x$31$ with a cross$-$sectional area of $9\mathrm{.}12\backslash (\backslash$mathrm$\{$in$\}^\{2\}\backslash )$ and a depth of $15\mathrm{.}88$ in $.$ The second moment of inertia is $\backslash (375\backslash$mathrm$\{$in$\}^\{4\}\backslash ).$ The beam is subjected to a uniformly distributed load of $\backslash (1000\backslash$mathrm$\{$lb$\}.\backslash$mathrm$\{$ft$\}\backslash )$ and a point load of $500$ ib $.$ The modulus of elasticity of the beam is $\backslash (\backslash$mathrm$\{$E$\}=29\backslash$mathrm$\{$x$\}106\backslash$mathrm$\{$lb$\}/\backslash$mathrm$\{$in$\}^\{2\}\backslash ).$ Using manual calculation $($FEA$),$ determine

a$)$ the vertical displacement at node $3$ and

b$)$ the rotations $($slopes$)$ at nodes $2$ and $3.$

c$)$ The vertical displacement at the point where the $500$ lb force is applied.

d$)$ Compute the reaction forces at nodes $1$ and $2$ and reaction moments at node $1.$

Take the length of element $(1)$ as $10$ ft and the length of element $(2)$ as $5$ ft $.$ Note that the $500$ Ib force is applied at the middle of element $(2).$