h = 0:.01:1.;
n = length(h);
r = zeros(n,1);

for i=1:n-1
    [r(i),~,~] = mybis(@per, h(i), 0, 1);
end

plot(h, r);
