% PR = (Md_MSF_seclatent_steam)./(Ncqsum); % Performance Ratio
PR = rand(365,1);
PR_ave=sum(PR)/365; % Average Performance Ratio
hold on
plot(PR)
plot(PR_ave*ones(1,365));
ylabel("PR")
xlabel("Days")
title("Performance Ratio profile")
legend('PR','PR-ave')