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