close all;
clc;
clear;

N=10000;
x=randn(N,1);
for i=1:100


F_norm=quad(@NumF,-20,20);  %real distribution = NumF/F_norm

%Mean

MeanX(i) = 1/N*sum((x(i)).*NumF(x(i))./(F_norm*Normal((x(i)))));
VarX(i)  = 1/N*sum((((x(i))-MeanX(i)).^2).*NumF((x(i)))./(F_norm*Normal((x(i)))));
KurX(i)  = 1/N*sum((((x(i))-MeanX(i)).^4).*NumF((x(i)))./(F_norm*Normal((x(i)))))/VarX(i)^2;

end

figure(1);
plot(MeanX);
figure(2);
plot(VarX);
figure(3);
plot(KurX);