#include <iostream>
#include <complex>
using namespace std;
//use units
#include <boost/typeof/std/complex.hpp>
#include <boost/units/cmath.hpp>
#include <boost/units/systems/si.hpp>
#include <boost/units/systems/si/prefixes.hpp>
#include <boost/units/systems/si/codata/universal_constants.hpp>
using namespace boost::units;
using namespace boost::units::si;
using namespace boost::units::si::constants::codata;
//use units end

double psi(
      int n                             // n - order of wavefunction
    , const quantity<mass> m            // m - mass of particle
    , const quantity<frequency>  omega  // oscillator frequency
)
{
    typedef derived_dimension<length_base_dimension,-2>::type  reciprocal_area;

    quantity<energy>            E      ( hbar*omega*(n+0.5)       );        // energy of oscillator
    quantity<reciprocal_area>   lambda ( 2.0*m*E / (hbar*hbar)    );        // 
    quantity<reciprocal_area>   alpha  ( m*omega/root<2>(hbar)    );
    double                      lambda_per_alpha ( lambda / alpha );
    complex<double> i(0,1);
    complex<double> a(0);
}

int main()
{
    return false;
}