ABD_matrix.m

function ABD = ABD_matrix(C_r, z)
    % ABD_matrix Returns ABD matrix.
    %
    % Arguments:
    %   C_r : Cell array of size n with stiffness matrices in ply CS.
    %   z   : Array of size n+1 with the locations of the ply edges.
    %
    % Returns:
    %   ABD : ABD matrix.
    A = zeros(3,3);
    B = zeros(3,3);
    D = zeros(3,3);
    if length(z) - 1 == length(C_r)
        for i=1:length(C_r)
            A = A + % COMPLETE THESE STATEMENTS
            B = B +
            D = D +
        end
    end
    ABD = [A, B; B, D];
end