function result= estimator(X, n) result= (sum(X) +1)/(n+2); endfunction