x=[2 4 7 10]

x =

     2     4     7    10

a=x(1)-((x(2)-x(1))/2)

a =

     1

n=length(x)

n =

     4

b=x(n)+((x(n)-x(n-1))/2)

b =

   11.5000

% x-ova os

u=[a x b]

u =

    1.0000    2.0000    4.0000    7.0000   10.0000   11.5000

% y-nova os

for i=1:n,
F(i)=[(2*i-1)/(2*n)];
end
F

F =

    0.1250    0.3750    0.6250    0.8750

% y-nove hodnoty

F=[0 F 1]

F =

         0    0.1250    0.3750    0.6250    0.8750    1.0000

n1=leghth(F);

%vykreslenie pcldf funkcie

for i=1:n1-1,
plot([u(i), u(i+1)],[F(i),F(i+1)]),hold on, shg
end

F=[0 F]

F =

         0    0.1250    0.3750    0.6250    0.8750    1.0000

u

u =

    1.0000    2.0000    4.0000    7.0000   10.0000   11.5000

hold on
grid on
plot(u,F,'+r')
shg

% pomocou interpolaci sa daju urcit suradnice 

%urcenie=doratanie y-ovej suradnice 

>> interp1(u,F,[4.7])

ans =

    0.4333

>> interp1(u,F,[7])

ans =

    0.6250

% doratanie x-ovej suradnice co znamena ze ked
% mame kvantil 0.5 to je vysledok 5.5 

>> interp1(F,u,[0.5])

ans =

    5.5000
%skuska median by mal vyst rovnaky

>> median(u)

ans =

    5.5000