#include <boost/units/systems/si.hpp>
#include <boost/units/systems/si/codata_constants.hpp>
#include <boost/units/io.hpp>

//kompilacja z flaga -std=c++0x tj.:
//g++ -std=c++0x boost_2.cpp

using namespace boost::units;
using namespace boost::units::si;
using namespace std;

// p = rho R_N2 T
quantity<mass_density> rho(quantity<temperature> T, quantity<pressure> p){
	auto M_N2 = 28 * 1e-3 * kilograms / moles;
	auto R_N2 = constants::codata::R / M_N2;
	return p / R_N2 / T;
}

int main()
{
	cerr << "rho = " << rho(300 * kelvins,1000e2 * pascals)<< endl;

	return 0;
}