t

theta = rand(1000, 1);
a = rand(1, 2000);
b = rand(1, 2000);
c = rand(1, 2000);

n = size(a,2);
N = size(theta, 1);

s = bsxfun(@minus, bsxfun(@times, a, theta), b);
Information = 1e5 * repmat(a .^ 2, N, 1) .* ((exp(s) ./ (1-exp(s))) .^ 2) .* ...
  bsxfun(@times, 1-c, 1 ./ bsxfun(@plus, c, exp(s)));