User:Ginevra Cochran/Formal Report/MATLAB

X = 4;
 hydrogen = [650,653,653.2, 653, 652.9;486.05,486.05,485.9,485.95,486;434.4,434.1,434,434.4,433.8;403.9,403.95,404,404,404.05]
 deuterium = [652.9,652,652.1,651.9,651.9;485.9,485.95,486,486,485.9;433.9,433.9,433.9,433.7,433.8]
 %Y1 = mean(hydrogen(1,:));
 Y1 = [mean(hydrogen(2,:))-std(hydrogen(2,:))/sqrt(5),mean(hydrogen(2,:)),mean(hydrogen(2,:))+std(hydrogen(2,:))/sqrt(5)]
 Y2 = [mean(deuterium(2,:))-std(deuterium(2,:))/sqrt(5),mean(deuterium(2,:)),mean(deuterium(2,:))+std(deuterium(2,:))/sqrt(5)]
 %Y3 = mean(hydrogen(3,:));
 %Y4 = mean(hydrogen(4,:));
 %Y = [Y1 Y2 Y3 Y4]
 %Z1 = mean(deuterium(1,:));
 % Z2 = mean(deuterium(2,:))
 %Z3 = mean(deuterium(3,:));
 %Z = [Z1 Z2 Z3 0]
 %E1 = std(hydrogen(1,:))/sqrt(5);
 E2 = std(hydrogen(2,:))/sqrt(5)
 %E3 = std(hydrogen(3,:))/sqrt(5);
 %E4 = std(hydrogen(4,:))/sqrt(5);
 %E = [E1 E2 E3 E4]
 %F1 = std(deuterium(1,:))/sqrt(5);
 F2 = std(deuterium(2,:))/sqrt(5)
% F3 = std(deuterium(3,:))/sqrt(5);
 %F = [F1 F2 F3]
 plot(X,Y1(1),'rx',X,Y1(2),'rx',X,Y1(3),'rx',X,Y2(1),'bo',X,Y2(2),'bo',X,Y2(3),'bo')