clear all
load RanDat.txt
whos
gnumbers=RanDat(:,1);
plot(gnumbers)
N=length(gnumbers)
avg=mean(gnumbers)
sigma=std(gnumbers)
ibig=find(gnumbers>6);
length(ibig)
iplot=1:N;
plot(iplot(ibig),gnumbers(ibig),'or')
plot(gnumbers)
hold on
plot(iplot(ibig),gnumbers(ibig),'or')
xlabel('List number in file')
ylabel('Number in file')
title('Some Gaussian-Distributed Random Numbers')
close(gcf)

subplot(2,2,1)
plot(gnumbers)
axis([0 N 1 9])
bins=2:0.25:8;
freq=hist(gnumbers,bins);
subplot(2,2,2)
plot(bins,freq,'o')
subplot(2,2,3)
stairs(bins,freq)
close(gcf)

plot(gnumbers)
title(['Mean=',num2str(avg)])
title(['Mean=',num2str(avg),'\pm',num2str(sigma)])
close(gcf)

datenum('7-Sept-2010')
datenum('8-Sept-2010')