#include <algorithm>
#include <cmath>
#include <functional>
#include <iostream>
#include <list>
#include <numeric>
#include <vector>
using namespace std;
template <class T>
void print(T& c){
for( typename T::iterator i = c.begin(); i != c.end(); i++ ){
std::cout << *i << endl;
}
}
int main( )
{
const float a[] = { 1, 1.3, 1.5, 0.9, 0.1, 0.2};
vector<float> data( a,a + sizeof( a ) / sizeof( a[0] ) );
print( data );
float mean = accumulate( data.begin(), data.end(), 0.0f )/ data.size();
vector<float> zero_mean( data );
transform( zero_mean.begin(), zero_mean.end(), zero_mean.begin(),bind2nd( minus<float>(), mean ) );
// compute the sample standard deviation
float deviation = inner_product( zero_mean.begin(),zero_mean.end(), zero_mean.begin(), 0.0f );
deviation = sqrt( deviation / ( data.size() - 1 ) );
// erase all points more than three standard deviations greater than the mean
const int num_deviations = 3;
vector<float>::iterator end = remove_if( data.begin(), data.end(), bind2nd( greater<float>(),mean + num_deviations * deviation ) );
data.erase( end, data.end() );
}