Consider the following $7$ steps, write a Matlab program to solve the frame structure below.

Step $1$: Data Preparation and Input.

Step $2$: Computation of element matrices.

Step $3$: Assembly of elements' stiffness matrices and elements' force vectors in the global stiffness matrix and global force vector.

Step $4$: Imposition of boundary conditions such as prescribed displacement. This is ussaully carried out simultaneously with the assembly of global matrix in Step $3.$

Step $5$: Solution of the global system of equations for the nodal unknowns.

Step $6$: Computation of secondary variables such as stresses and strains from displacements

Step $7$: Print and$/$or plot desired results.