stat1_1.m

Satatistical program for variancy

clear;clc;

N=100000;
a = -3.;
b=3.;
W = 0.01;
X = randn(1,N);
%N = length(X)
prog=(a:W:b);
[nr nc] = size(prog);

for i1=1:nc
    ref(i1)=(1./sqrt(2*pi))*exp(-prog(i1)^2/2);
    K=find(X > prog(i1));
    kc = length(K);
    F(i1)=kc;
end

for i1=1:nc-1
    p(i1)=(F(i1)-F(i1 + 1))/(N*W);
end

prog1=prog(1:nc-1);
plot(prog1,p,'r*',prog,ref,'g')