clear all;
% accepted value is 1.76 E11 C/Kg

% Data entry for constant I, Day 1

% Constant magnet current
I1=1.5;
V1 = [339.8 350 359.8 369.7 440];
R1 = [3.4 3.0 2.7 2.9
      3.0 3.5 3.0 3.3
      3.0 3.5 3.4 4.4
      3.7 3.2 3.7 3.2
      4.8 3.2 3.5 4.8];

B1 = 7.8E-4 * I1;
R1 = R1 .* .01;
R1 = R1';
for i=1:5
    R_1(i)=mean(R1(:,i));
    EM1(i) = (2*V1(i))/(R_1(i)*B1)^2;
end
display(sprintf('Day 1, Constant magnet current: E/m = %.3E with std dev of %.3E', mean(EM1), std(EM1)));  

%Data entry for constant V Day 1
V2=375;
I2=[1.2,1.3,1.4,1.5,1.6];
B2 = 7.8E-4 .* I2;
R2=[3.8 4.8 4.0 4.5; 3.8 4.2 3.9 4.2; 3.8 4.6 3.9 4.4; 3.8 4.0 3.5 4.1; 3.4 3.2 3.2 3.5];
R2 = R2 .* .01;
R2 = R2';
for i=1:5
    R_2(i)=mean(R2(:,i));
    EM2(i) = (2*V2)/(R_2(i)*B2(i))^2;
end

display(sprintf('Day 1, Constant accelerating voltage: E/m = %.3E with std dev of %.3E', mean(EM2), std(EM2)));

% Day 2 data processing

% Constant magnet current
I3=1.35;
V3 = [441.2 430.5 420.0 409.8 399.9 390.1 380.0 369.8 360.0 350.0];
R3 = [4.0 6.0 4.0 5.5;
      4.0 5.5 4.0 4.5;
      3.9 5.4 3.9 5.3;
      3.9 5.3 3.8 5.3;
      3.8 5.2 3.8 5.1;
      3.7 5.1 3.7 5.0;
      3.5 5.0 3.6 4.8;
      3.4 4.9 3.5 4.8;
      3.3 4.8 3.5 4.8;
      3.3 4.8 3.4 4.5];


B3 = 7.8E-4 * I3;
R3 = R3 .* .01;
R3 = R3';
for i=1:10
    R_3(i)=mean(R3(:,i));
    EM3(i) = (2*V3(i))/(R_3(i)*B3)^2;
end
display(sprintf('Day 2, Constant magnet current: E/m = %.3E with std dev of %.3E', mean(EM3), std(EM3)));  

% Constant magnet current data set
V4=360;
I4=[1.22 1.32 1.42 1.52 1.62 1.72 1.27 1.37 1.47 1.57];
R4 = [3.4 5.2 3.8 5.2;
      3.4 5.1 3.6 5.0; 
      3.0 4.7 3.3 4.7; 
      2.9 4.5 3.5 4.4; 
      2.9 3.9 3.0 4.0; 
      2.4 3.5 3.0 3.7; 
      3.5 5.2 3.5 5.3; 
      3.4 5.0 3.4 4.8; 
      3.2 4.8 3.5 4.5; 
      3.0 4.3 3.3 4.4];

B4 = 7.8E-4 .* I4;
R4 = R4 .* .01;
R4 = R4';
for i=1:10
    R_4(i)=mean(R4(:,i));
    EM4(i) = (2*V4)/(R_4(i)*B4(i))^2;
end
display(sprintf('Day 2, Constant accelerating voltage: E/m = %.3E with std dev of %.3E', mean(EM4), std(EM4)));