# 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')

```