Develop a 'MATLAB program' for $1$D bar problem as given in the figure above and carry out the following tasks:

$($a$)$ Global stiffness matrix by assembling the element stiffness matrices

$($b$)$ Formation of global load vector

$($c$)$ Application of boundary conditions

$($d$)$ Solution for displacement at every node and reaction forces