function res = pr1zad1func(x, y)
res = [ y(2);
        y(3);
        x.^2-8*y(3)-2*y(2)-5*y(1) ];
end
