jakauma2.m

% Tuotetaan tuhat satunnaislukua normaalijakautuneesti keskiarvolla x0
% ja keskihajonnalla sigma.
x0 = 2.5
sigma = 2.5
x = normrnd(x0,sigma,1000,1);

% Esitetään satunnaisluvut graafisesti
figure(1)
plot(x,'o')

% Tuotetaan histogrammi 
N = 30;
rajat = linspace(-5,10,N+1);
osumat = histc( x, rajat );

% Esitetään histogrammi graafisesti
figure(2)
bar( rajat, osumat, 'histc' )

% Esitetään samassa kuvassa miten keskimäärin normaalijakautuneet
% luvut jakautuvat
hold on
f = @(x) exp(-(x-x0).^2/(2.*sigma^2))./sqrt(2*sigma^2*pi)
xx = linspace(-20,20);
h = rajat(2)-rajat(1);
plot(xx,1000*f(xx)*h)
hold off