function bifurcation1
h=0.1;
p=2;
x=0:0.01:4;
y=(1-x).*(h.^p+x.^p)./x.^(p-1);
plot(y,x,'-b','LineWidth',2)
xlabel('parameter c')
ylabel('x')
xlim([0 0.5])
ylim([0 1])
axis on

