function ci=f(x,n, percent)
% binciold(x,n,percent) computes the old confidence interval for the binomial parameter p,
% at the confidence "percent" in the range 0-100, and
% x is the number of successes as a row vector, n is the number of trials.
z=norminv(1- (100-percent)/200);
u=x/n+ z*sqrt(x/n.*(1-x/n)/n);
l=x/n - z*sqrt(x/n.*(1-x/n)/n);
ci=[l; u];