load Nominalmatrices.mat ;
n=17; % n is number of candidate measurments
for i = 3:n
    [B,sset,ops,ctime,flag]=pb3wc(G1,F,Wd,Wn,Juu,Jud,i,tlimit,1);
    meas{i,1} = sset;
    loss(i) = B;
end
figure(1)
plot(3:17,loss(3:end),'r-*','LineWidth',2)
xlabel('No. of Measurements');
ylabel('Loss')

figure(2)
plot(3:17,loss(3:end)*1e8/(365*24),'r-*','LineWidth',2)
xlabel('No. of Measurements');
ylabel('Loss $/hr')
gtext({'individual measurements'})