>> x=[2 4 7 10]

x =

     2     4     7    10

>> n=length(x)

n =

     4
%vytvorenie hodnot pre y-ovu osu

>> for i=1:n
y(i)=(i/n);
end
>> y

y =

    0.2500    0.5000    0.7500    1.0000

>> y=[0 y]

y =

         0    0.2500    0.5000    0.7500    1.0000

>> x

x =

     2     4     7    10

%vyuzity vzorec pre vypocet a

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

a =

     1

%vyuzity vzorec pre vypocet b

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

b =

   11.5000

>> u=[a x b]


%intervaly pre x-ovu os

u =

    1.0000    2.0000    4.0000    7.0000   10.0000   11.5000

>> y

y =

         0    0.2500    0.5000    0.7500    1.0000

>> n1=length(u)

n1 =

     6

>> axis([0 12 -0.1 1.1])

%vykreslenie

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