First version
Working, but not beautiful.

TODO
1) Replace estimation procedure with one function
2) Add monte-carlo testing scripts (clear and good one)
3) Check reasons for bias (long story, you know it)
Dev.branch of functions
Added separate bandwidth options, rlm options list.
precise.h turned off, lm() replaced by rlm() for initial residuals
First try in pushing
Merge branch 'master' of github.com:Loiisso/modal.param.0.8