diff(diff((1/2)*log(2*pi) + log(s) + (x-m)^2/(2*s^2),m),s)