% file orderstat.m
% Version of 2/3/95
% Calculates order statistics from dbn function
n  = input('Enter the sample size  ');
t  = input('Enter a row matrix of times  ');
F  = input('Distribution function F as a function of t  ');
m = length(t);
FK = zeros(m,n);
for i = 1:m
  FK(i,:) = cbinom(n,F(i),1:n);
end
disp('      t        F(t)     Values of FK(k) for k = 1:n ')
disp([t' F' FK])