Finite Element Analysis (FEA) Truss Analysis

I have a 29x4 matrix named [E] which represents the element properties of a truss structure with 29 elements and 15 nodes. Each element has a correspnding row in matrix [E] with four columns [A, L, Theta, E] which represent cross sectional area, length, angle, and Young Modulus values. How can I write a function in MATLAB that uses [E] as an input to return an elemental stiffness matrix [K] which is 4 x 4 x N in size (N=#of elements=29)?

% Angle
Theta1 = 45;
Theta2 = 90;
Theta3 = 180;
Theta4 = 135;

% Length
L = 4;

% Cross-sectional Area
% A_all = 40(20) + 80(20) = 2400mm
% A_inside = 70(10) + 10(35) = 1050mm
A = (1350*10^-3);

% Young's Modulus
Ealuminum = (8*10^10);
Esteel = (2.1*10^11) ;

% [E] = [A, L, Theta, E]
[E] = [ A, L, Theta2, Esteel;
        A, L, Theta2, Esteel;
        A, L, Theta2, Esteel; 
        A, L, Theta3, Esteel;
        A, L, Theta4, Esteel;
        A, L, Theta3, Esteel;
        A, L, Theta4, Esteel;
        A, L, Theta3, Esteel;
        A, L, Theta4, Ealuminum;
        A, L, Theta3, Ealuminum;
        A, L, Theta3, Ealuminum;
        A, L, Theta2, Ealuminum;
        A, L, Theta4, Ealuminum;
        A, L, Theta3, Ealuminum;
        A, L, Theta2, Ealuminum;
        A, L, Theta4, Ealuminum;
        A, L, Theta3, Ealuminum;
        A, L, Theta2, Ealuminum;
        A, L, Theta4, Ealuminum;
        A, L, Theta3, Ealuminum;
        A, L, Theta2, Ealuminum;
        A, L, Theta4, Ealuminum;
        A, L, Theta3, Ealuminum;
        A, L, Theta2, Ealuminum;
        A, L, Theta3, Ealuminum;
        A, L, Theta4, Ealuminum;
        A, L, Theta3, Ealuminum;
        A, L, Theta2, Ealuminum;
        A, L, Theta2, Ealuminum;];

Answer 1

[KSSV]

Useful file exchange function: https://in.mathworks.com/matlabcentral/fileexchange/14313-truss-analysis