x=0:0.1:2.5;
y1=-x.^3+2*x.^2+2*x+1;
y2=0*ones(1,size(x,2));
X=[x,fliplr(x)];
Y=[y1,fliplr(y2)];
fill(X,Y,'b');