>
plot([sqrt(x), 4, 4+(x-16)/8, 4+(x-16)/8-(x-16)^2/512],x=1..30,color=[red,yellow,green,blue]);
>