function [B R] = bootstrap(Obbl,rho,T)

[m n] = size(rho);

for i = 1:m
  for j = 1:i
    A(i,j) = rho(i)*(T(j+1)-T(j));
  end
  A(i,i) = A(i,i) + 1;
end 

B = A\Obbl;

R = [0;-log(B/100)]./(T-T(1));

endfunction