so=page_screen_output;
if so
   more
end
#page_screen_output=0;
x=[16 34 53 75 93 120];
wblplot(x)

# Benchmark eta = 76.318, beta = 1.4301 rho = 0.9956
       
 
x=[2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 59 61 67]';

c=[0 1 0 1 0 1 1 1 0 0 1 0 1 0 1 1 0 1 1]';

figure;

[h p]=wblplot(x,c)

# Benchmark eta = 82.02, beta = 0.895
#pause


x=[16, 34, 53, 75, 93, 120, 150, 191, 240 ,339];

# Benchmark eta = 146.2, beta = 1.20 & rho = 0.998703
# Reliasoft eta = 146.2545 beta 1.1973 rho = 0.9999

figure;

[h2 p]=wblplot(x,[],[],0.05)


x=[46 64 83 105 123 150 150];
c=[0 0 0 0 0 0 1];
f=[1 1 1 1 1 1 4];


figure;

[h3 p2]=wblplot(x,c,f,0.05)


figure;
disp("subtract 30.92 from x to simulate a 3 parameter wbl with gamma = 30.92")

wblplot(x-30.92,c,f,0.05)



if so
   more
end

