% Octave session on p. 193

format long;

function y = intfun (t, alpha)
x=2*t./(1+t); dx=2./(1+t).^2; y = dx.*x.^ alpha;
end
alfa=pi/4; [ts,ws]=q_ossinf(16,alfa/2,1/8);
[tc,wc]=q_ossinf(16,alfa/2,5/8);
wf=ws.*intfun(ts,alfa); q=sum(wf), sum(abs(wf)>eps*q)

x=2*ts./(1+ts);
wf1=wc.*intfun(tc,alfa).*(1+tc)/2; wf2=ws.*log(x).*intfun(ts,alfa);
dq=sum(wf1)+sum(wf2), sum(abs(wf1)>eps*dq), sum(abs(wf2)>eps*dq)
disp('\n')
