Question: The command "singular" returns a wrong result


s,k,mu,sigma are parameters.k is real number

f := unapply(k^k*(int(exp(-x^2)/((s-chi)*exp(sqrt(2)*sigma*x+mu)+k)^k, x = -infinity .. infinity))/(sqrt(Pi)*chi^(k+1)), chi, s, k, mu, sigma)

proc (chi, s, k, mu, sigma) options operator, arrow; k^k*(int(exp(-x^2)/((s-chi)*exp(2^(1/2)*sigma*x+mu)+k)^k, x = -infinity .. infinity))/(Pi^(1/2)*chi^(k+1)) end proc

(1)

singular(f(chi, s, k, mu, sigma), chi)

{chi = 0}, {chi = (exp(2^(1/2)*sigma*x+mu)*s+k)/exp(2^(1/2)*sigma*x+mu)}

(2)

``

The result includes x,which is integration variable

Download singular.mw

Please Wait...