"produce the same graph" The code below produces the same graph.
However
- it may be disputed whether it's simpler
- less intermediate results are saved
%%
c = 5;
retint = 0;
%%
cac = { '-go','Displayname','gap = 12(8,6,4)1'
'-ro','Displayname','gap = 14(9,6,3)1'
'-bo','Displayname','gap = 16(10,6,2)1' };
num = [ [12:-1:8,6,4:-1:1]+retint
[14:-1:9 6 3:-1:1]+retint
[16:-1:10 6 2:-1:1]+retint ];
%%
figure(11);
hold on
%%
for jj = 1 : 3
dist = log(num(jj,:));
eta = exp(-c*abs(dist-reshape(dist,[],1)));
discrim = 1./sum(eta,1);
plot( discrim, cac{jj,:} )
end
%%
axis([0,length(dist(1,:)),0,1.02*max(max(discrim))])
xlabel('Serial Position','fontsize',14)
ylabel('Discriminability','fontsize',14)
title({'SIMPLE Memory Model'...
'\fontsize{9}\color{blue}Length = 10; Various Gaps'})
lgd = legend ('Location','northwest');
lgd.NumColumns = 1;
lgd.FontSize = 9;
set(gcf, 'color', 'white');
hold off
