%Returns area under h-function over interval [k*step,k*step+step]:
%KVK: This assgnment syntax, "retv=calculat.." drives us crazy: it seems upward: in reality function takes value retv.
function retv=betaPortion(k, step)
t=k*step;
power=t+1;
p2=power*power;
p4=p2*p2;
value1=1/(p2*p4); %6
t=t+step;
power=t+1;
p2=power*power;
p4=p2*p2;
retv=value1-1/(p2*p4);
end
Copyright (C) 2009 Konstantin Kirillov