Generator N(0,1)f := x->1/sqrt(2*Pi) * exp(-x^2 / 2); <0, sup(sqrt(f(g))> x <inf(z-a)/b * f(z), sup(zamietacia3 := proc (dlzka) local yH, zH, zD, f,i,v,zoznam,Y,Z,W,efektivita;
yH := (2*Pi)^(-1/4);
zD := -(2/Pi)^(1/4) * exp(-1/2);
zH := (2/Pi)^(1/4) * exp(-1/2);
f := x->1/sqrt(2*Pi) * exp(-x^2 / 2);
i := 0;
v := 0;
zoznam := [];
while i<dlzka do
v := v + 1;
Y:=stats[random, uniform[0, yH]](1);
Z:=stats[random, uniform[zD, zH]](1);
W := evalf(Z/Y);
if(evalf(Y^2) <= evalf(f(W))) then
i := i+1;
zoznam := [op(zoznam), W];
# print("Prijimame" || W);
else
# print("Zamietame");
end
end;
efektivita := dlzka / v;
print("Efektivita", evalf(efektivita));
return(zoznam);
end:zamietacia3(150);stats[statplots, histogram](%);zamietacia4 := proc (dlzka) local f,i,v,zoznam,R,X,y1,y2,Z,W,efektivita;
f := x->1/sqrt(2*Pi) * exp(-x^2 / 2);
i := 0;
v := 0;
zoznam := [];
while i<dlzka do
v := v + 1;
R:=stats[random, uniform[-1, 1]](1);
X:=(R+1)/(2*root(2 * Pi, 4));
y1 :=-root(1-R^2, 2) * (1/(root(2 * Pi, 4)));
y2 :=+root(1-R^2, 2) * (1/(root(2 * Pi, 4)));;
Z:=stats[random, uniform[y1, y2]](1);
W := evalf(Z/X);
if(evalf(X^2) <= evalf(f(W))) then
i := i+1;
zoznam := [op(zoznam), W];
# print("Prijimame" || W);
else
# print("Zamietame");
end
end;
efektivita := dlzka / v;
print("Efektivita", evalf(efektivita));
return(zoznam);
end:zamietacia4(150);stats[statplots, histogram](%);