% Below are the 
%  Deviations from 1.00, in percent, of the
%  ratios of strength value (to run 1 measurement)
%  for currents of Imag = [1,10,20,40, 50,100]
%
%  Pct= (Ratio_SL2_SL1 -1.0 ).*100.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
Imag = [    1      10         20       40           50       100 ];
%
% run,  pct1,    pct10,    pct20,     pct40,      pct50,     pct100
%
A = [...   
 6  -0.13115  -0.0045659  0.0040725  0.0066228   0.006628    0.006631;
 7  -6.2954   -0.73569   -0.32524   -0.15721    -0.13353    -0.063887;
 8  -4.138    -0.4835    -0.21038   -0.095358   -0.078736   -0.04024 ;
 9  -3.4213   -0.40286   -0.17583   -0.078822   -0.064608   -0.035667;
10  -3.002    -0.35548   -0.1549    -0.067595   -0.054493   -0.030919;
11  -2.423    -0.29635   -0.13485   -0.064347   -0.05381    -0.0365;
12  -2.3883   -0.28435   -0.12382   -0.052861   -0.042018   -0.024738;
13  -2.0525   -0.24251   -0.1034    -0.040994   -0.031363   -0.01883;
14  -1.8769   -0.2268    -0.10061   -0.044686   -0.036212   -0.024469;
15  -1.4699   -0.12821   -0.034596  -0.0052227  -0.0049988  -0.014278;
20  -34.005   -5.0694    -2.7339    -1.5318      -1.2905    -0.81138;
21  -32.409   -4.7328    -2.5022    -1.3585      -1.1303    -0.68068;...
];
[n,p]=size(A)
t=1:n;
run=A(t,1);
% pct1=A(t,2);  
figure(1)
subplot(2,1,1)
plot(run,A(t,2))
legend('1A')
xlabel('Run Number');
ylabel('Percent Dev from 1.00')
title('QE-N13 Studies Mar-Apr 2006 - AWW')
grid on
%
subplot(2,1,2)
hold all
for i=3:4
plot(run,A(t,i));
end
legend('10A','20A')
xlabel('Run Number');
ylabel('Percent Dev from 1.00')
grid on
%%%%%%%%%%%5
figure(2)
hold all
for i=5:7
plot(run,A(t,i));
end
legend('40A','50A','100A')
xlabel('Run Number');
ylabel('Percent Dev from 1.00')
title('QE-N13 Studies Mar-Apr 2006 - AWW')
grid on
hold off
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%markertype =cellstr(['-ok';'-+b';'-+g';'-*k';'-or';'-+y';'-+c'])
ii=2:p;
figure(3)
subplot(2,1,1)
axis([0. 100. -7. 1.])
hold all
for j=1:10 
plot(Imag,A(j,ii)) %,char(markertype(ii)),'LineWidth',2)
end
legend(int2str(A(1,1)), int2str(A(2,1)), int2str(A(3,1)),...
int2str(A(4,1)), int2str(A(5,1)), int2str(A(6,1)),...
int2str(A(7,1)), int2str(A(8,1)), int2str(A(9,1)),...
int2str(A(10,1)))
xlabel('Current/A');ylabel('Percent Dev from 1.00')
title('QE-N13 Studies Mar-Apr 2006 - AWW')
grid on
subplot(2,1,2)
axis([0. 100. -0.5 0.1])
hold all
for j=1:10 
plot(Imag,A(j,ii)) %'-s','MarkerSize',10) %,char(markertype(ii)),'LineWidth',2)
end
xlabel('Current/A');ylabel('Percent Dev from 1.00')
grid on
%
figure(4)
subplot(2,1,1)
axis([0. 100. -35. 0.1])
hold all
for j=10:n
plot(Imag,A(j,ii))
end
legend(int2str(A(10,1)), int2str(A(11,1)), int2str(A(12,1)))
xlabel('Current/A');ylabel('Percent Dev from 1.00')
title('QE-N13 Studies Mar-Apr 2006 - AWW')
grid on
subplot(2,1,2)
axis([0. 100. -2. 0.1])
hold all
for j=10:n
plot(Imag,A(j,ii))
end
%legend(int2str(A(10,1)), int2str(A(11,1)), int2str(A(12,1)))
xlabel('Current/A');ylabel('Percent Dev from 1.00')
grid on
%%%%%%%%%%%%%%%%%%%%%%

