Advertisement
Guest User

Untitled

a guest
Feb 19th, 2019
85
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.81 KB | None | 0 0
  1. Assuming[{a > 0, b > 0, n > 0, sse > 0, b > a},
  2. Integrate[(1/(Sqrt[2 Pi sigma^2]))^n
  3. Exp[-(x - mu)^2/(2 sigma^2)] PDF[UniformDistribution[{a, b}], sigma], {sigma, a, b}]]
  4.  
  5. ([Pi]^(-n/2) (1/(mu - x)^2)^(-(1/2) + n/2) (Gamma[1/2 (-1 + n), (mu - x)^2/(2 a^2)] -
  6. Gamma[1/2 (-1 + n), (mu - x)^2/(2 b^2)]))/(2 Sqrt[2] (a - b))
  7.  
  8. FullSimplify@D[g, mu]
  9.  
  10. fDeriv2[x_, mu_, n_, a_, b_] := (-2^(3/2 - n/2) E^(-((mu - x)^2/(2 a^2))) ((mu - x)^2/a^2)^(1/2 (-1 + n)) +
  11. 2^(3/2 - n/2) E^(-((mu - x)^2/(2 b^2))) ((mu - x)^2/b^2)^(
  12. 1/2 (-1 + n)) - (-1 + n) (Gamma[1/2 (-1 + n), (mu - x)^2/(2 a^2)] -
  13. Gamma[1/2 (-1 + n), (mu - x)^2/(2 b^2)]))/((mu - x) (Gamma[1/2 (-1 + n), (mu - x)^2/(2 a^2)] -
  14. Gamma[1/2 (-1 + n), (mu - x)^2/(2 b^2)]))
  15.  
  16. fDeriv2[0.1, 0.2, 100, 2, 4]
  17.  
  18. ComplexInfinity
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement