%% Continuous-Time Sparse Second-Order Model
% For this example, consider the sparse matrices for the 3-D beam model subjected to an impulsive point load at its tip in the file sparseBeam.mat.
% Extract the sparse matrices from sparseBeam.mat.

load('sparseBeam.mat','M','K','B','F','G','D');

% Create the mechss model object by specifying [] for matrix C, since there is no damping.

sys = mechss(M,[],K,B,F,G,D)

% The output sys is a mechss model object containing a 3-by-1 array of sparse models with 3408 degrees of freedom, 1 input, and 3 outputs.
% You can use the spy command to visualize the sparsity of the mechss model object.

spy(sys)

% Copyright 2020 The MathWorks, Inc.