# The following code is one method for simulating the average score from 3-pt shots that the basketball player 
# from our classroom example earns. This method simulates the individual shots.    
number_simulations := 100000:   
scores := array(1..number_simulations):                                             
for i from 1 to number_simulations do   
scores[i] := 0;
u := rand()/10^12;  
if (u < .25) then shots := 0 elif (u<.5) then shots := 1 elif (u<.75) then shots := 2 else shots := 3 end if;  
for j from 1 to shots do;
u := rand()/10^12;
if (u < .1) then scores[i] := scores[i] + 3 end if;
end do:
end do:  
evalf(add(scores[i], i=1..number_simulations)/number_simulations);