if exist ('OCTAVE_VERSION', 'builtin') pkg load statistics; graphics_toolkit('qt'); end nrep = 10; for i = 1 : nrep tmp = rand(1,3); [A, ia] = min(tmp); [B, ib] = max(tmp); tmp([ia,ib]) = []; C = tmp; x = linspace(0, 1, 100); o_x = ones(size(x)); xi = rand(1); % distributions with C=A args = {A, B, A}; n_ABA = tst_cdf_inv(x, args{:}) n_ABA_p = tst_pdf(xi, A, args{:}) % distributions with C=B args = {A, B, B}; n_ABB = tst_cdf_inv(x, args{:}) n_ABB_p = tst_pdf(xi, A, args{:}) % distributions with A < C < B args = {A, B, C}; n_ABC = tst_cdf_inv(x, args{:}) n_ABC_p = tst_pdf(xi, A, args{:}) % distributions with A < A+eps < B args = {A, B, A + (B-A)*1.e-10}; n_ABaE = tst_cdf_inv(x, args{:}) n_ABaE_p = tst_pdf(xi, A, args{:}) % distributions with A