Используйте Matlab, чтобы получить максимальное значение функции плотности обратного гауссовского распределения

Question

Новый для Matlab, и вот мой вопрос.

Вопрос: Необходимо использовать метод отклонения, чтобы найти максимальное значение функции плотности обратного распределения Гаусса, как показано ниже.

% Using Rejection Method to generate standard normal N(0,1). 
% Check how many realizations have been rejected 

rng(123);
N = 100000;
Count = 0;
Rejection = 0;
X = [];
tic %Timer
while Count<= N
   U1 = -4+(4-(-4))*rand; % U1: Uniform(-4, 4)
   U2 = 0.5*rand; % U2: Uniform(0, 0.5)
   if U2<= 1/(sqrt(2*pi))*exp(-U1^2/2)
       Count = Count + 1;
       X = [X,U1];
   else
      Rejection = Rejection + 1; 
   end

Пожалуйста, дайте мне знать, если вам нужна какая-либо другая информация. Спасибо!