%Function that plots and saves the figures of the conditional pdfs of parameters that are
%used for computing MIRI. It takes in input the struct of the model, the
%folder where the results are saved, the output node of interest and the
%number of independent realizations


function plotpdf_param(folder,variable_name,model,Nr)

nominal_parameters=model.nominal_parameters;
parameters_name=model.parameters_name;


    try
        load(fullfile(folder,variable_name,'pdf_param.mat'));
        for i=1:length(nominal_parameters)
            figure
            for k=1:Nr
                plot(xbin_param{1,k}(i,:),ks_param{1,k}(i,:),'b',xbin_param{2,k}(i,:),ks_param{2,k}(i,:),'r');
                hold on
            end
            xlabel(parameters_name{1,i},'Interpreter','none');
            ylabel('pdf');
            legend('Conditional upper density','Conditional lower density');
            savefig(fullfile(folder,variable_name,parameters_name{1,i}));
        end
    catch ME
        errordlg(ME.message)
    end
end