Dopasowanie wykresu funkcji do histogramu

Hej!
Usiłuję napisać projekt, w którym mam dopasować do histogramu wykres funkcji, jednak kompletnie nie wiem, jak ją wpisać. Konkretnie chodzi o funkcję μ(A) z poniższego linka. Mogłabym prosić o jakieś wskazówki?
https://en.wikipedia.org/wiki/Marchenko%E2%80%93Pastur_distribution

Dodam, że do rysowania histogramu użyłam ggplot2, jednak stat_function zawodzi

Pokaż kod, choć w R nie piszę to spróbujemy poszukać błędu ;)

Jednak całkiem źle używałam stat_function, ale dalej problemem jest wpisanie jakoś poniższej funkcji (załącznik) w kod - przecież sigma nie jest w żaden sposób zdefiniowana, nawet gdy mamy lambdę. Niżej wysyłam co mam, ale oczywiście gdy chcemy wykonać zwraca object 'sgm' not found

ni <- function(x, lam){
  lamplus = sgm^2(1+sqrt(lam))
  lamminus = sgm^2(1-sqrt(lam))
  if (x > lamminus && x < lamplus){
    result = 1/(2*pi*sgm^2) * sqrt((lamplus - x)*(x - lamminus))/(lam*x)
  }else{
    result = 0
    }
}