I read this thread, but it didnt really help.
test.mw （i noticed an error in the expression pj, now fixed!）
See this file above.
The function needs to be maximized is "loglik", with constraints,
f0>0, sigma>0 and mu is FREE. Only three variables.
Using some calculus and by changing variables, I was able to use Gauss–Hermite quadrature, and do some approximation, to get the values:
I am wondering if Maple can do this straight away. By that, I mean, without doing much further, how could I get the values? (most efficiently and quickly)