This version seems to be bug free in terms of compiling and running

Need to fix the 'na' and all '0' parameters returned by 'lm_fit'