warning off MATLAB:nchoosek:LargeCoefficient
[N,t]=meshgrid(72:1000,0:.01:1);
z=ones(size(N));
y=[53 57 66 67 72];
for i=1:size(N,1),
for j=1:size(N,2),
for k=1:5,

z(i,j)=z(i,j)*nchoosek(N(i,j),y(k))*t(i,j)^y(k)*(1-t(i,j))^(N(i,j)-y(k));

end 
end 
i 
end


cz=cumsum(z(:))/sum(z(:));

for i=1:1000,

ind=min(find(cz>rand));
samp(i,:)=[t(ind) N(ind)];
end

length(find(samp(:,2)>100))/1000