Incoherent estimates and bug for truncated normal functions #2

Hello,

Trying to run the examples for rNormal_trunc_ab and eNormal_trunc_ab that you provide in the documentation at https://search.r-project.org/CRAN/refmans/ExtDist/html/Normal_trunc_ab.html, I obtain incoherent estimates :

Console :

X <- rNormal_trunc_ab(n= 500, mu= 2, sigma = 5, a = 1, b = 2)
est.par <- eNormal_trunc_ab(X); est.par

Parameters for the Normal_trunc_ab distribution.
(found using the numerical.MLE method.)

Parameter Type Estimate
mu shape -93.752527
sigma shape 20.881979
a boundary 1.001197
b boundary 1.996527

...clearly far off for the shape parameters !?

Then other tests end up in errors, for example setting the bound "b" to 3 :

Console :

X <- rNormal_trunc_ab(n= 500, mu= 2, sigma = 5, a = 1, b = 3)
est.par <- eNormal_trunc_ab(X); est.par

Error in trans.est.par[[k]] : subscript out of bounds

Thanks for looking into this,
Benoît