x1=linspace(-1, 0);
x2=linspace(-1, 0.25);
x3=linspace(.25, (13/20));
x4=linspace(.25, (13/20));
x5=linspace((13/20),1);
x6=linspace(1, 2);

y=2.*((x1.*(x1-1).*(x1-4))./(3.*x1-8)).^.5;
y2=.25.*((((9-4.*x2).*(13-20.*x2)).^.5)-(9-4.*x2));
y3=.25.*((9-4.*x3)-((9-4.*x3).*(13-20.*x3)).^.5);
%y4=.25.*((9-4.*x4)+((9-4.*x4).*(13-20.*x4)).^.5);
y4=.25.*((9-4.*x3)+((9-4.*x3).*(13-20.*x3)).^.5);
y5=((2.*(x5-1).*(x5-4).*(x5-9))./(x5-5)).^.5;
y6=2.*((x6.*(x6-1).*(x6-4))./(3.*x6-8)).^.5;
figure
plot(x1,y,x2,y2,x3,y3,x4,y4,x5,y5,x6,y6);
ylim([0 2])
color = [0, 0, 1];
hold on;
area(x1,y,'FaceColor','b')
hold on
a = fill([x2 x3], [y2 y3], color);
hold on
a = fill([-1 x4], [2 y4], color);
hold on
a = fill([x5 x6], [y5 y6], color);
hold on
x = [-1 -1 .65];
y = [1.92808 2 1.6];
fill(x,y,'b','EdgeColor','b')