function [M,K,F,H,R,faceWithBCs] = discretizedModelHarmonicNedelec(self,coefstruct,edges,edge_elem,numEdges,frequency)

    [p,~,t] = self.Mesh.meshToPet();

    % pass dummy val to solution to fill in state vector appropriately
    % (needed for axisymmetric version)
    femat = formGlobalMatrix(self,p,t,coefstruct,zeros(size(p,2),1),[],'kfm',edge_elem,numEdges,frequency);
    K = femat.K;
    M = femat.M;
    F = femat.F;

    [H,R,faceWithBCs] = assembleDirichletBCNedelec(self,edges);

end