function [u, modelForResults] = solveHarmonic(self,coefstruct,frequency)
% Internal only.

%   Copyright 2022 The MathWorks, Inc.

self.PDESystemSize = 1;
for freq = 1:length(frequency)
    u(:,freq) = solveNedelecHarmonic(self,coefstruct,frequency(freq)); %#ok<AGROW>
end
if self.IsTwoD
    self.PDESystemSize = 2;
else
    self.PDESystemSize = 3;
end
modelForResults = self;
end