x=[-10,-9,-8.5,-8,-7,-6,-5,-4,-3,-2,-1,0,0.5,0.8,1,2,3,4,5,6,7,8,9,10,10.1];
y0=x.^3-6*x.^2+5*x-3
y=[-1653,-1263,0.485375648722841,-939,-675,-465,-303,-183,-99,-45,-15,-3,0.957166948242946,0.800280468888800,-3,-9,-15,-15,-3,27,81,165,285,447,0.141886338627215];
y3 =0.7453*x.^3 - 5.406*x.^2 + 7.774*x + 3.825
y2 =- 4.833*x.^2 + 59.84*x - 11.47
y4 =- 0.04962*x.^4 + 0.7783*x.^3 - 0.7526*x.^2 + 6.738*x - 40.94
plot(x,y0,'.k',x,y,'-.g',x,y3,'--y',x,y2,'.m',x,y4,'-c')