function plotfunc (hamso,khoangve)
    %
    x=linspace(khoangve(1),khoangve(2),1000);
    y=feval(hamso,x);
    plot(x,y,'b-','LineWidth',3);
    grid on;
    