% life = wblrnd(0.5,2,100,1);  % Simulate lifetimes.
% [p,ci] = wblfit(life)
% 
% t = 0:0.1:4.5;
% h1 = exppdf(t,0.6267) ./ (1-expcdf(t,0.6267));
% h2 = wblpdf(t,2,2) ./ (1-wblcdf(t,2,2));
% plot(t,h1,'--',t,h2,'-')


