using Plots

plot(x->sin(x)-(x-x^3/3))