close all
clear
clc
AME_VoltSum=0;
AME_VAngleSum=0;
JMeasurementSum=0;
AME_PDSum=0;
AME_QDSum=0;
AME_mVoltSum=0;
AME_mPDSum=0;
AME_mQDSum=0;
totalTimeSum=0;
N=200;
loopN=1;
while 1
    [JMeasurement,AME_Volt,AME_VAngle,AME_PD,AME_QD,AME_mVolt,AME_mPD,AME_mQD,isConverged,totalTime]=OPF();
    if isConverged==0
        %I=I-1;
        continue
    end
    loopN=loopN+1;
    if loopN>N
        break;
    end
% [JMeasurement]=run();
    AME_VoltSum=AME_VoltSum+AME_Volt;
    AME_VAngleSum=AME_VAngleSum+AME_VAngle;
    AME_PDSum=AME_PDSum+AME_PD;
    AME_QDSum=AME_QDSum+AME_QD;
    AME_mVoltSum=AME_mVoltSum+AME_mVolt;
    AME_mPDSum=AME_mPDSum+AME_mPD;
    AME_mQDSum=AME_mQDSum+AME_mQD;
    totalTimeSum=totalTimeSum+totalTime;
end
JMeasurementSum=JMeasurementSum/N;
AME_VoltSum=AME_VoltSum/N;
AME_VAngleSum=AME_VAngleSum/N;
AME_PDSum=AME_PDSum/N;
AME_QDSum=AME_QDSum/N;
AME_mVoltSum=AME_mVoltSum/N;
AME_mPDSum=AME_mPDSum/N;
AME_mQDSum=AME_mQDSum/N;
totalTimeSum=totalTimeSum/N;