Есть несколько способов вычислить этот интеграл.
пусть .
Тогда пусть и подставим :
TrigSubstitutionRule(theta=_theta, func=sin(_theta), rewritten=cos(_theta)**2, substep=RewriteRule(rewritten=cos(2*_theta)/2 + 1/2, substep=AddRule(substeps=[ConstantTimesRule(constant=1/2, other=cos(2*_theta), substep=URule(u_var=_u, u_func=2*_theta, constant=1/2, substep=ConstantTimesRule(constant=1/2, other=cos(_u), substep=TrigRule(func='cos', arg=_u, context=cos(_u), symbol=_u), context=cos(_u), symbol=_u), context=cos(2*_theta), symbol=_theta), context=cos(2*_theta)/2, symbol=_theta), ConstantRule(constant=1/2, context=1/2, symbol=_theta)], context=cos(2*_theta)/2 + 1/2, symbol=_theta), context=cos(_theta)**2, symbol=_theta), restriction=And(u < 1, u > -1), context=sqrt(-u**2 + 1), symbol=u)
Если сейчас заменить ещё в:
Перепишите подынтегральное выражение:
пусть .
Тогда пусть и подставим :
TrigSubstitutionRule(theta=_theta, func=sin(_theta), rewritten=cos(_theta)**2, substep=RewriteRule(rewritten=cos(2*_theta)/2 + 1/2, substep=AddRule(substeps=[ConstantTimesRule(constant=1/2, other=cos(2*_theta), substep=URule(u_var=_u, u_func=2*_theta, constant=1/2, substep=ConstantTimesRule(constant=1/2, other=cos(_u), substep=TrigRule(func='cos', arg=_u, context=cos(_u), symbol=_u), context=cos(_u), symbol=_u), context=cos(2*_theta), symbol=_theta), context=cos(2*_theta)/2, symbol=_theta), ConstantRule(constant=1/2, context=1/2, symbol=_theta)], context=cos(2*_theta)/2 + 1/2, symbol=_theta), context=cos(_theta)**2, symbol=_theta), restriction=And(u < 1, u > -1), context=sqrt(-u**2 + 1), symbol=u)
Если сейчас заменить ещё в:
Теперь упростить:
Добавляем постоянную интегрирования:
Ответ:
(0.0 + 2.60671648277426j)