better with alpha & beta
Modifying this case. However, it does not work.
new procedure for fitting
better
linear constraint
more