clear all
close all
%% 
load sim_1.mat 
%% 
subplot(3,3,1:2)
plot(time,sup,'k','LineWidth',3), hold on 
plot(time,storexl,'r','LineWidth',1), hold on 
xlabel('Time')
ylabel('Supervisor')
title('Supervisor (black) vs. Trained Network (red)')
xlim([0,300])
patch([imin*dt,imax*dt,imax*dt,imin*dt],[-1,-1,10,10],'b','edgealpha',0,'facealpha',0.1)
text(imax*dt-50,9,'RLS ON')
text(imax*dt+1,9,'RLS OFF')

subplot(3,3,3)
plot(time,sup,'k','LineWidth',3), hold on 
plot(time,storexl,'r','LineWidth',1), hold on 
xlim([200,300])
title('Zoom')
patch([imin*dt,imax*dt,imax*dt,imin*dt],[-1,-1,10,10],'b','edgealpha',0,'facealpha',0.1)
text(imax*dt-30,9,'RLS ON')
text(imax*dt+1,9,'RLS OFF')


%%
subplot(3,3,4:5)
for j = 1:5
    plot(time,storerl(:,j)*0.5+j), hold on 
end
title('Firing Rates')
ylabel('Neuron Index')
xlabel('Time')
ylim([0,6])
subplot(3,3,6)
for j = 1:5
    plot(time,storerl(:,j)*0.5+j), hold on 
end
xlim([200,300])
xlabel('Time')
ylabel('Neuron Index')
title('Zoom')
ylim([0,6])



%% 
subplot(3,3,7)
for j = 1:5
plot(time,storet(:,j)*0.4+j), hold on 
end
ylabel('Neuron Index')
xlabel('Time')
title('Firing Rates, g=1.5xN^{-1/2}')
ylim([0,6])
subplot(3,3,8)
for j = 1:5
plot(time,store2(:,j)*0.4+j), hold on 
end
title('Firing Rates, g=1.9xN^{-1/2}')
ylim([0,6])
ylabel('Neuron Index')
xlabel('Time')
%% 
subplot(3,3,9)
plot(storext2(:,1),storext2(:,2),'g','LineWidth',2), hold on 
plot(storext(:,1),storext(:,2),'k'), hold on 
legend('g=1.9xN^{-1/2}Reservoir','g=1.5xN^{-1/2}Reservoir')
xlim([-3,3])
ylim([-3,3])
f1 = figure(1)

%%
set(f1,'position',[0,0,1000,1000])
print(f1,'f1.svg','-dsvg','-painters')