* change argument optarg to _optarg so that it doesn't colide with the global
  variable optarg
* crop sum over m when value gets very small
* change default LMAX_FACTOR to 5
* use other scaling for XiRL
* very simple wrapper for python
change order to n=80
solve integrals analytically
changed style of integration
add simple suite for unittests
add simple unittests
add unittests