diff --git a/+SimplusGT/+Modal/AppModalAnalysis.m b/+SimplusGT/+Modal/AppModalAnalysis.m index d45a23c7..e7789063 100644 --- a/+SimplusGT/+Modal/AppModalAnalysis.m +++ b/+SimplusGT/+Modal/AppModalAnalysis.m @@ -56,15 +56,15 @@ for k=1:length(ApparatusType) if ApparatusType{k} <= 89 % Ac apparatus IMR.Type(count) = ApparatusType{k}; - IMR.IMRVal(count) = SigmaMag/norm(Residue{k},"fro")*norm(ZmVal{k},"fro"); + IMR.IMRVal(count) = SigmaMag/(norm(Residue{k},"fro")*norm(ZmVal{k},"fro")); count = count+1; elseif ApparatusType{k} >= 1000 && ApparatusType{k} <= 1089 % Dc apparatus IMR.Type(count) = ApparatusType{k}; - IMR.IMRVal(count) = SigmaMag/norm(Residue{k},"fro")*norm(ZmVal{k},"fro"); + IMR.IMRVal(count) = SigmaMag/(norm(Residue{k},"fro")*norm(ZmVal{k},"fro")); count = count+1; elseif ApparatusType{k} >= 2000 && ApparatusType{k} <= 2009 % Interlink apparatus IMR.Type(count) = ApparatusType{k}; - IMR.IMRVal(count) = SigmaMag/norm(Residue{k},"fro")*norm(ZmVal{k},"fro"); + IMR.IMRVal(count) = SigmaMag/(norm(Residue{k},"fro")*norm(ZmVal{k},"fro")); count = count+1; elseif ApparatusType{k} == 90 || ApparatusType{k} == 1090 % infinite bus: let IMR=inf IMR.Type(count) = ApparatusType{k};