function [m,c,k,b,f,g,d,mStored] = getMCKBFGD(D)
% Returns matrices of zero-order Pade approximation.
% These matrices are computed by setting all internal delays 
% to zero and adding the delay signals to the state vector
% (no elimination, no singularity of algebraic loops).

%   Copyright 2020 The MathWorks, Inc.

% Zero internal delays (error free)
D = zeroInternalDelay(D);

% Read data
c = D.c;
k = D.k;
mStored = D.m;
if isempty(mStored)
   m = speye(size(k));
else 
   m = mStored;
end
b = D.b;  
f = D.f;
g = D.g;
d = D.d;  
