For a 2 dimensional 4 noded membrane element, it is required to implement a FE code through MATLAB. I am incapable of implementing these three steps. Can you implement a MATLAB code to build up a stiffness matrix for 2D linear membrane element which has 4 nodes?

Step 2, solution: build_K_matrix (2/4) --- Integrate stiffness matrix for id = 1 : nint_p for j = 1 : nint_p Loop over x- and y-direction for numerical integration xi = xGauss (MG); n = xGauss ( iG); --- Step 1 Get kinematic gradients KINEMATICS = get_kinematic_gradients ( ELEMENT, xi, n); Returns: Jacobian of the transformation and its determinant, grad N Step 2 --- Get elasticity tensor ELEMENT = get_C_matrix( ELEMENT, MATERIAL, KINEMATICS); Returns: the elasticity tensor (ELEMENT.C) Step 3 --- Get B matrix [ Bc, ~) get_B_matrix( KINEMATICS); Returns: the Bc matrix for evaluating the material stiffness Step 2, solution: build_K_matrix (2/4) --- Integrate stiffness matrix for id = 1 : nint_p for j = 1 : nint_p Loop over x- and y-direction for numerical integration xi = xGauss (MG); n = xGauss ( iG); --- Step 1 Get kinematic gradients KINEMATICS = get_kinematic_gradients ( ELEMENT, xi, n); Returns: Jacobian of the transformation and its determinant, grad N Step 2 --- Get elasticity tensor ELEMENT = get_C_matrix( ELEMENT, MATERIAL, KINEMATICS); Returns: the elasticity tensor (ELEMENT.C) Step 3 --- Get B matrix [ Bc, ~) get_B_matrix( KINEMATICS); Returns: the Bc matrix for evaluating the material stiffness