function y=mf(x)
y=x.^4-0.5*x.^3-28*x.^2+140;
end
%x=-5:0.1:6;
%y=x.^4-0.5*x.^3-28*x.^2+140;
%plot(x,y,'-k'), grid
%Максимум функции на интервале [-2 2]
%y=-mf(x);
%[x,y]=fminbnd(@mf,-2,2)