I need help with MATLAB on this assignment, Please provide code and graphs

Write a MATLAB program to plot the exact values from part a) using 100 points. On the same graph, plot the three FEA nodal solution values from part c). In your program, assign non-zero values of your choice for A, f, and L, and an appropriate value for E based on your chosen material type. You will need to assign values such that the engineering strain at the free end is no more than 2%. The purpose of this is to keep the problem linear, for better comparison with ABAQUS in part g). (10 points). e) Next, create a new MATLAB program to repeat part d), but for a Finite Element discretizaton using a user-specified number of nodes, N. You can solve the global matrix system, [K]{u}={f}, for the nodal displacement vector {u}, using MATLAB as follows (30 points):

By using the inverse of K: u = inv(K) * f Or by Gaussian elimination: u = K\f However, before you do this, you will need to make [K] non-singular (to prevent rigid body motion of the rod). Since the displacement for node 1 is known, with value of zero (fixed to wall), you will make [K] non-singular by eliminating the first row and first column, as follows:

K = K(2:N, 2:N) where N is the original number of rows and columns in [K]. You will also need to eliminate the 1st row of load vector {f} to make sure the size of each matrix or vector in the linear system is compatible. This above procedure effectively throws out the first equation in the original matrix system since the displacement u1=0 is known, and we only need the number of unknowns to match the number of equations, so the node 1 column is also deleted in the above operation.

NOTE: A, E(x) E3 forxL2 Force per unit length, icx EE for x -L2 f(x) = f for x -L2 f(x) = f for x