Finished implementing number-of-gaussians iterative algorithm found in the paper. Now I just need to debug it and integrate it with fitting functional mean dependence