clear all;
data = load('noise1');
y = data(:,2);
x = data(:,1);
m = mean(y);
variance  = var(y);
standard = std(y);
rms = rms(y);
bins =standard*[-4.25 -3.75 -3.25 -2.75 -2.25 -1.75 -1.25 -0.75 -0.25,0.25 0.75 1.25 1.75 2.25 2.75 3.25 3.75 4.25];
b=hist(data,bins);
one_std = b(8)+b(9)+b(10)+b(11);
two_std = b(12)+b(13)+b(14)+b(15);
percent_one_std = 100.0 * (one_std/5000);
percent_two_std = 100.0 * (two_std/5000);
fit_to_gauss_cdf(y,100)
fit_to_gauss_pdf(y,100)
% hist(y,bins); grid;
% title([num2str(percent_one_std),'qqqq = ',num2str(percent_two_std)]);
% hold on;
% plot(x,y);
% title(['\sigma = ', num2str(standard,4),'\sigma^2 = ', num2str(variance,4),'RMS = ', num2str(rms,4),'m = ',num2str(m,4)]);
% ylabel('Amplitude (V)')
% xlabel('Time (Sec)');