I want to calculate the derivative of Bessel, I do
hh := proc (`ϖ`) options operator, arrow; (`ϖ`-V*m*Pi/LL)^2/C[f]^2-m^2*Pi^2/LL^2 end proc; H := proc (`ϖ`) options operator, arrow; piecewise(0 <= hh(`ϖ`), sqrt(hh(`ϖ`)), hh(`ϖ`) < 0, sqrt(-hh(`ϖ`))) end proc; Y := proc (`ϖ`) options operator, arrow; piecewise(0 <= hh(`ϖ`), BesselJ(n, R*sqrt(hh(`ϖ`))), hh(`ϖ`) < 0, BesselI(n, R*sqrt(-hh(`ϖ`)))) end proc; YY := proc (`ϖ`, z) options operator, arrow; piecewise(0 <= hh(`ϖ`), BesselJ(n, z), hh(`ϖ`) < 0, BesselI(n, z)) end proc; www := proc (`ϖ`) options operator, arrow; (D[2](YY))(`ϖ`, R*sqrt(hh(`ϖ`))) end proc; PF := proc (`ϖ`) options operator, arrow; rho[f]*(`ϖ`-V*m*Pi/LL)^2*Y(`ϖ`)/(H(`ϖ`)*www(`ϖ`)) end proc
Is it right or not?
Thank you very much!
