3kn/mm 1 kn/mm 5 mm 3 mm 10 mm consider the above 2d plane strain linear elastic structural problem. the left side has roller boundary conditions allowing the structure to move only in the x2 direction, while the bottom side has roller boundary conditions allowing the structure to move only in the xi direction. assuming a thickness of lmm, an elastic modulus of 100 gpa and poisson's ratio of 0.3, write an fem code in matlab to solve this problem. your code should: (a) use four-node quadrilateral plane-strain elements with the mesh as outlined below, (b) read in the nodal coordinates and element connectivity (you should create an input file that contains this information), (c) formulate the shape functions, stiffness matrix, and local force vectors accounting for local-global transformation and invoking gauss quadrature for the integration, (d) solve for the nodal displacements, (e) plot the deformed structure clearly, indicating the locations of the nodes and element edges (scale all displacements so that the maximum x displacement is lmm in order to clearly illustrate the shape of the deformed structure), and compute the σ3x and gy stresses at the four integration points of element a (f)