xg = a:0.1:b; x(end + 1)=b; plot (x,f(x),'-',xg,f(xg),'-');