% parametri
b = 5;
t = 3;
L = -5;
U = 5;

% mantise
M = zeros(b^t,1);
i = 1;
for c1 = 0:b-1
    for c2 = 0:b-1
        for c3 = 0:b-1
            m = c1*b^-1 + c2*b^-2 + c3*b^-3;
            M(i) = m;
            i = i+1;
        end
    end
end

% normalizirana stevila
bm = b^(t-1);
d = U-L+1;
Xpn = zeros((b-1)*bm,d);
for j = 0:d-1
    Xpn(:,j+1) = M(bm+1:end)*b^(L+j);
end
Xpn = Xpn(:);
Xn = [-Xpn(end:-1:1); Xpn];

% denormalizirana stevila
Xpd = M(2:bm)*b^L;
Xd = [-Xpd(end:-1:1); Xpd];

% predstavljiva stevila (brez 0, Inf, -Inf, NaN)
X = [Xn(1:end/2); Xd(1:end/2); Xpd; Xpn];

%1) Delez denormaliziranih stevil
disp('Delez denormaliziranih stevil')
disp(length(Xd)/length(X))

%2) Koliko predstaviljivih normaliziranih stevil je manjsih od pi?
%>>moznost 1
disp('Stevilo normaliziranih stevil manjsih od pi') 
disp(sum(Xn<pi));

%>>moznost 2
a = Xn < pi;
ai = find(a);
manjsaOdPi=ai(end);

%3) Povprecna razdalja med predstavljivimi stevili, ki so za manj kot 1
%oddaljena od pi
disp('Povprecna razdalja')
disp(mean(diff(X(abs(X-pi)<1))))

%4) fl(pi)
a = X < pi;
ai = find(a);
i = ai(end);
xm = X(i);
xv = X(i+1);

if pi-xm < xv-pi
    flpi = xm;
elseif pi-xm > xv-pi
    flpi = xv;
else
    disp('obe stevili sta enako oddaljeni od pi')
end

disp('fl(pi)')
disp(flpi)