% C. Talarico
% exampleIdivider.m

clear all; close all; clc;

filename = 'idivider.txt';
[nvars blocks points pData] = ltsExp2Matlab(filename);

mR = blocks*nvars; % number of data sets (rows of pData)
nR = points;       % number of points in each data set (columns of pData)

figure(1);
hold on;
for i=1:nvars:mR % steps nvars latge
  Iin = pData(i,:);
  IR1_IR2 = pData(i+1,:);
  IR2 = pData(i+2,:);
  IR1 = pData(i+3,:);
  plot(Iin*1e6,IR1*1e6,'linewidth',2,'color','b');
  plot(Iin*1e6,IR2*1e6,'linewidth',2,'color','m');
end
ax = gca;
set(ax, 'fontName', 'Arial Narrow', 'fontsize', 16);
xlabel('I_{in} [\muA]','fontsize',16, 'fontName', 'Arial Narrow');
ylabel('I_{R1} and I_{R2} [\muA]','fontsize',16, 'fontName', ...
    'Arial Narrow');
title('I_{R1} and I_{R2} vs. I_{in}','fontsize',18,'fontName', ...
    'Arial Narrow','fontWeight','bold');
legend('I_{R1}','I_{R2}', 'location','best');

figure(2);
plot(Iin*1e6,IR1_IR2,'linewidth',2,'color','b');

ax = gca;
set(ax, 'fontName', 'Arial Narrow', 'fontsize', 16);
xlabel('I_{in} [\muA]','fontsize',16, 'fontName', 'Arial Narrow');
ylabel('I_{R1}/I_{R2}','fontsize',16, 'fontName', 'Arial Narrow');
title('I_{R1}/I_{R2} vs. I_{in}','fontsize',18,'fontName', ...
    'Arial Narrow','fontWeight','bold');