Physics307L F09:Labs/Balmer/balmerlab1.m

%Using the values taken in runs 1 - 4 % Not using runs 5 & 6 for this as that was taken to show % variablilty of data with single calibration and weights data unfairly % towards last calibration

V1Raw = [410.8 410 412 410]; %Violet 1 data V2Raw = [435.5 434 436 434]; %Violet 2 data GBRaw = [488.8 486 489 486]; %Green Blue data RRaw = [666 657 666 656];   %Red data

V1 = mean(V1Raw); V2 = mean(V2Raw); GB = mean(GBRaw); R = mean(RRaw);

for n=3:10 %Candidates for R from Violet 1 data RV1(n) = V1 * (1/4 - 1/n.^2); RV1(n) = 1/RV1(n); %Candidates for R from Violet 2 data RV2(n) = V2 * (1/4 - 1/n.^2); RV2(n) = 1/RV2(n);

%Candidates for R from green blue data RGB(n) = GB * (1/4 - 1/n.^2); RGB(n) = 1/RGB(n); %Candidates for R from Red data RR(n) = R * (1/4 - 1/n.^2); RR(n) = 1/RR(n); end RV1 = RV1.* 1E9; RV2 = RV2.* 1E9; RGB = RGB.* 1E9; RR = RR.* 1E9;

n = 1:10; figure(1); plot(n(3:10), RV1(3:10), 'b'); hold on; plot(n(3:10), RV2(3:10),'b--'); plot(n(3:10), RGB(3:10),'g'); plot(n(3:10), RR(3:10),'r'); hold off; grid on;

Rydberg = RR(3) + RGB(4) + RV2(5) + RV1(6); Rydberg = Rydberg/4 R(n) = Rydberg;

figure(1); hold on; plot(n(3:10), R(3:10), 'k');hold off; title(sprintf('Averaged Data from Balmer lab: R = %.3E', Rydberg)); legend('Violet 1', 'Violet 2', 'Blue-Green', 'Red', 'Calculated R', -1); xlabel('Value of n2'); ylabel('Possible Value of R');