p=-0.5;
h=[];b=[];
syms  x
for a=0:0.1:5
    y1=[(-log(x)).^(2*a-1)].*[(x.^(-p)-1)/p];
    z1=int(y1,x,0,1);
    y2=(-log(x)).^a;
    z2=int(y2,x,0,1);
    y3=(1-x).*(x.^(-p)-1)/p;
    z3=int(y3,x,0,1);
    y4=(1-x).*(-log(x)).^a;
    z4=int(y4,x,0,1);
    z=[a*z1/z2-(a-1)*z3/z4]^2;
    h=[h;z];
    b=[b;a];
end
plot(b,h)