stat2.m

Satatistical program for variancy

clear;clc;

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

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

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

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