x1=[0.1:0.01:0.4];
x2=[0.4:0.01:0.8];
n1=size(x1,2);
n2=size(x2,2);
y1=zeros(n1,1);
y2=zeros(n2,1);
for i=1:n1
  y1(i)=-0.1*x1(i)^5-x1(i)^4-8*x1(i)^3+12*x1(i)^2-3*x1(i);
end;
for i=1:n2
  y2(i)=-0.1*x2(i)^5-x2(i)^4-8*x2(i)^3+12*x2(i)^2-3*x2(i)+0.4;
end;
y3=[y1(n1):0.01:y2(1)];
n3=size(y3,2);
x3=zeros(n3,1);
for i=1:n3
  x3(i)=0.4;
end;
hold off;
plot(x1,y1,'k',x2,y2,'k',x3,y3,':k');
