#include <iostream>
#include <boost/numeric/interval.hpp>
#include <cmath>

using namespace std;
using namespace boost::numeric;
using namespace boost::numeric::interval_lib;

template<class T, class Policies>
ostream& operator <<(std::ostream &os,const interval<T, Policies> &i) {
  return os << "[" << i.lower() << "," << i.upper() << "]" << endl;
}

int main(void) {
  typedef boost::numeric::interval<double, policies<save_state<rounded_transc_std<double> >, checking_base<double> > > I;
  
  I alpha;
  
  alpha.assign(-2.1, -2.1);
  cout << "sin(alpha) = " << sin(alpha) << endl;
  
  alpha.assign(-2.1 - M_PI / 2.0, -2.1 - M_PI / 2.0);
  cout << "cos(alpha) = " << cos(alpha);
  
  return 0;
};