close all
clear
clc
AME_VoltSum=0;
AME_VAngleSum=0;
JMeasurementSum=0;
AME_PDSum=0;
AME_QDSum=0;
N=1000;
for I=1:N
    [JMeasurement,AME_Volt,AME_VAngle,AME_PD,AME_QD,AME_mVolt,AME_mPD,AME_mQD,isConverged]=OPF();
    if isConverged==0
        I=I-1;
        continue
    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;
end
JMeasurementSum=JMeasurementSum/N;
AME_VoltSum=AME_VoltSum/N;
AME_VAngleSum=AME_VAngleSum/N;
AME_PDSum=AME_PDSum/N;
AME_QDSum=AME_QDSum/N;