changed to numerical solution, from analytical. Using BFGS algorithm to optimize
plotting bootstraps works, unfound scaling issue with MDS