olive2682562
2011-03-28, 10:27
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)
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)