∫d0− segundo x + unC−−−−−−−√k1(− segundo x + unC−−−−−−−√)d x
sustituyendo:
− segundo x + unC= t
−Cb∫- segundo re+ unCaCt√k1(t√)dt _
sustituyendo:
t =X2
−Cb∫- segundo re+ unC√aC√2X2k1( X )d x
Por partes:
− 2k0(−bd _C+aC−−−−−−−−√) re+ 2abk0(−bd _C+aC−−−−−−−−√) −2abk0(aC−−√) −4Cb∫- segundo re+ unC√aC√k0( x ) xd x
y por partes:
∫- segundo re+ unC√aC√k0( x ) xre x=k1(aC−−√)aC−−√−k1(- segundo re+ unC−−−−−−−√)- segundo re+ unC−−−−−−−√
entonces nosotros tenemos:
∫d0− segundo x + unC−−−−−√k1(− segundo x + unC−−−−−√)re x=1b( ( - 2bd _+ 2un )k0(- segundo re+ unC−−−−−√) +4k1(- segundo re+ unC−−−−−√)- segundo re+ unC−−−−−√do - 4k1(aC−−√)aC−−√do - 2k0(aC−−√) un)
Código de arce :
int(sqrt((-b*x + a)/c)*BesselK(1, sqrt((-b*x + a)/c)), x = 0 .. d) = ((-2*b*d + 2*a)*BesselK(0, sqrt((-b*d + a)/c)) + 4*BesselK(1, sqrt((-b*d + a)/c))*sqrt((-b*d + a)/c)*c - 4*BesselK(1, sqrt(a/c))*sqrt(a/c)*c - 2*BesselK(0, sqrt(a/c))*a)/b
Juan María
Apacs