format long
[x fval]=fsolve(@(h) (1-h).^0.461.*(0.47880*h+0.32256)-0.28385,    [-0.5 0.5])
h=linspace(-1,1,500);
y=(1-h).^0.461.*(0.47880*h+0.32256)-0.28385;
figure(1)
plot(h,y,'b',[min(h) max(h)],[0 0],'r')