function .SEheavyModel works now. I have two concerns:
1. solve(Jt): there is a caution that it returns a number
2. time to run: about 2-3 minutes.
Hope for your feedback.
I have 2 problems to ask for your advice:
1. solve(): in case the inverse of a matrix is error, how should we deal with it?
2. .heavy_likelihood_c does not run properly.
Thank you