blue = [57 106 177]./255;
green = [62 150 81]./255;
brown = [146 36 40]./255;
purple = [107 76 154]./255;
cl_colors = {blue, red, black, ...
cl_legend = {'For Borr', 'Inf Borr', 'For+Inf Br', 'For+Br+Save', 'Bridge Loan', 'For Save'};
cl_scatter_shapes = {'s','x','o','d','p','*'};
cl_linestyle = {'--','-',':','-.','--','-'};
cl_scatter_csizes = {10*it_sca_bs, 20*it_sca_bs, 10*it_sca_bs, 10*it_sca_bs, 5*it_sca_bs, 8*it_sca_bs};
cl_line_csizes = {1*it_line_bs, 2*it_line_bs, 1*it_line_bs, 1*it_line_bs, 1*it_line_bs, 2*it_line_bs};
it_x_groups_n = length(cl_scatter_csizes);
mat_y = rand([it_x_n, it_x_groups_n]);
mat_y = mat_y + sqrt(1:it_x_groups_n);
mat_y = mat_y + log(1:it_x_n)';
for it_fig = ar_it_graphs_run
it_graph_counter = it_graph_counter + 1;
ar_y = mat_y(:, it_fig)';
it_csize = cl_scatter_csizes{it_fig};
ar_color = cl_colors{it_fig};
st_shape = cl_scatter_shapes{it_fig};
st_lnsty = cl_linestyle{it_fig};
st_lnwth = cl_line_csizes{it_fig};
% plot scatter and include in legend
ls_chart(it_graph_counter) = scatter(ar_x, ar_y, it_csize, ar_color, st_shape);
% plot line do not include in legend
line.HandleVisibility = 'off';
line.LineStyle = st_lnsty;
line.HandleVisibility = 'off';
line.LineWidth = st_lnwth;
cl_legend{it_graph_counter} = cl_legend{it_fig};
legend(ls_chart, cl_legend, 'Location', 'southeast');
title('Optimal Savings');
ylabel('Savings Levels');
xlabel('Cash-on-Hand Today');