Listing 5

// Averaging
const int N(16);
UncertainPtr accum( constant(0.0) );
int i;
for(i=0;i<N;++i) {
   accum = accum + vm.read();
}
accum = accum / constant(N);
P = power( accum , constant(100) );
print(P);