reset; model hospital.mod; data hospital.dat; load amplgsl.dll; option solver cplex; param B := 2000; var count integer := 0; for {b in 1..B} { set S := {rand(1,n) | i in 1..726}; let count_b := card{i in S : sum{c in C} Total_Time_PSC[c,i] * X[c] + sum{h in H} Total_Time_CSC[h,i] * X[h] < Threshold}; let count := count + count_b; } let avg_count := count / B; let std_count := sqrt(sum{b in 1..B} (count[b]/726 - avg_count)^2 / (B-1)); let lower_bound := avg_count - 1.96 * std_count / sqrt(B); let upper_bound := avg_count + 1.96 * std_count / sqrt(B); printf "count/726 mean = %f\n", avg_count; printf "count/726 standard error = %f\n", std_count; printf "count/726 95%% confidence interval = (%f, %f)\n", lower_bound, upper_bound;