y=@(x) (-2*x.^2-5*x-3).*(x>=-4&x<-1) + (x.^2-1).*(x>=-1&x<1) + (-x.^2+4*x-3).*(x>=1&x<=4); x = -4:0.01:4; plot(x,y(x))
