function [u, modelForResults] = solveMSNedelec(self,coefstruct)
% Internal only.
% Create a temporary model and set its system size to 3, so we can
% process the three components of the magnetic potential vector, the
% primary variable of the PDE. This allows us to compute gradients (9
% components), which in turn provides us with magnetic flux density
% vector B = curl(Avec).

%   Copyright 2022-2023 The MathWorks, Inc.

modelForResults = createpde('electromagnetic','magnetostatic');
modelForResults.Geometry = self.Geometry;
modelForResults.MaterialProperties = self.MaterialProperties;
modelForResults.Sources = self.Sources;
modelForResults.VacuumPermeability = self.VacuumPermeability;
modelForResults.Mesh = self.Mesh;
modelForResults.PDESystemSize = 3;
u = solveNedelecMS3D(self,coefstruct);
end