function [ K,M,F ] = getKMF(self)
% getKMF Calculate the K and M matrices and the F vector for tetrahedral models.
% Models may be composed of 4-node and 10-node tetrahedron elements.
% This undocumented function may be changed or removed in a future release.

%       Copyright 2014-2021 The MathWorks, Inc.

% call builtin functions for actual computation
coefstruct.ElementsInSubdomain = 1:size(self.ElemConn,2);
coefstruct.NumElementsInSubdomain = size(self.ElemConn,2);
coefstruct.Coefficients.m = {[0]};
coefstruct.Coefficients.d = {[0]};
coefstruct.Coefficients.c = {[self.c]};
coefstruct.Coefficients.a = {[self.a]};
coefstruct.Coefficients.f = {[self.f]};


femat = formGlobalMatrix(self.thePde,self.Points,self.ElemConn,coefstruct,self.u,self.time,'akf');
M = femat.A;
K = femat.K;
F = femat.F;

end

