O que está pedido: calcular o valor numérico da integral dupla sobre a região \Omega.
Informações relevantes:
As duas primeiras são retas: x+y=1 e 2x+y=1.
As outras duas são curvas do tipo potência: y=x^{3} e y=3x^{3}.
A região \Omega é a área limitada por estas quatro curvas (duas retas e duas curvas cúbicas), ou seja, os contornos da região são exatamente essas quatro curvas.
A integranda é \dfrac{(1+2y)\,y^{2}}{x^{11}}.
Observação final (importante para o desenvolvimento): para organizar a integração e simplificar a integranda, vou fazer a mudança de variáveis
s = x+y,\qquad t = x.
Essa escolha é natural porque as retas envolvidas são expressas de maneira simples em termos de s e t (por exemplo x+y=1 vira s=1), e a outra dimensão t=x preserva a potência x^{11} do denominador em uma forma simples. No passo seguinte começo a transformação e a descrição da região em coordenadas (s,t).
onde \tilde\Omega é a imagem da região \Omega no plano (s,t).
No próximo passo eu descrevo a região \tilde\Omega em (s,t) — ou seja, determino os limites de integração — e escolho uma ordem conveniente de integração.]
Passo 3
Descrever a região \tilde\Omega no plano (s,t).
Recordando as quatro curvas em termos de (x,y):
x+y=1 vira s=1.
2x+y=1 pode ser escrita como s+t=1 (pois 2x+y = x+(x+y)=t + s); então 2x+y=1 é s+t=1 ou s = 1 - t.
y=x^{3} vira s - t = t^{3}, ou sejas = t + t^{3}.
y=3x^{3} vira s - t = 3t^{3}, ou sejas = t + 3t^{3}.
Portanto \tilde\Omega é a região do plano (s,t) delimitada pelas curvas verticais/oblíquas e pelas curvas s=t+t^{3} e s=t+3t^{3}. Em palavras: para cada t dentro de um certo intervalo, s fica entre duas curvas (a intersecção entre o intervalo entre as duas curvas cúbicas e o intervalo entre as duas retas).
É conveniente trabalhar com a variável w=s-t (que é exatamente y). De fato, colocando w=s-t obtemos w=y, e a integranda fica muito simples em termos de w e t:
A forma desses máximos e mínimos muda quando as curvas se cruzam; portanto temos que particionar o intervalo de t em subintervalos onde a ordem destas quatro funções é fixa. Os pontos críticos (valores de t onde há troca de ordem) são exatamente as raízes das equações que ocorrem quando duas dessas expressões coincidem, a saber:
A partir da análise da ordem das curvas para valores de t entre esses pontos (verificação feita testando valores intermediários), obtém-se a partição do intervalo de t em três subintervalos com limites de w bem definidos:
Para t\in[t_a,t_b]: temos t^{3}<1-2t<3t^{3}<1-t. Logo
w\in[\,1-2t,\;3t^{3}\,].
Para t\in[t_b,t_c]: temos 1-2t\le t^{3}\le 3t^{3}\le 1-t. Logo
w\in[\,t^{3},\;3t^{3}\,].
Para t\in[t_c,t_d]: temos 1-2t< t^{3}<1-t\le 3t^{3}. Logo
w\in[\,t^{3},\;1-t\,].
Portanto a integral dupla separa-se em soma de três integrais (uma em cada subintervalo de t), cada uma com o integrando simples em w e com divisão por t^{11}.
No próximo passo eu calculo as integrais internas em w (são polinomiais) e depois integro em t, simplificando com as relações cúbicas nos pontos extremos para cancelar as raízes e obter um valor racional.
Passo 4
Calcular as integrais internas em w (integração em relação a w primeiro).
A primitiva em w do integrando \dfrac{(1+2w)w^{2}}{t^{11}} (com t tratado como parâmetro) é simples: escrevendo (1+2w)w^{2}=w^{2}+2w^{3},
Observação: os termos envolvendo potências simples de t (como t^{-2} e t) surgem explicitamente nas integrais I_1,I_2,I_3 e são diretamente integráveis. Os termos ostensivamente “difíceis” são aqueles que envolvem \dfrac{(1-2t)^{k}}{t^{11}} e \dfrac{(1-t)^{k}}{t^{11}}; contudo, note que estes aparecem apenas em integrais sobre intervalos cujos extremos satisfazem relações cúbicas (por exemplo em t_b temos t_b^{3}=1-2t_b, em t_c temos 3t_c^{3}=1-t_c, etc). Isso permitirá, ao depois avaliar as primitivas, usar essas relações em pontos extremos para eliminar as raízes cúbicas e obter resultados racionais. No passo seguinte calculo as primitivas das partes simples e trato as peças que envolvem (1-2t) e (1-t), usando as relações cúbicas nos extremos para simplificar.
Passo 5
Integrar as partes simples (as que resultaram em potências simples de t) e organizar as contribuições de todas as três faixas.
Das expressões do passo anterior notamos que as contribuições com termos simples são:
De I_1 obtemos o termo \displaystyle\int_{t_a}^{t_b}\Bigl(9t^{-2}+\frac{81}{2}t\Bigr)\,dt.
De I_2 obtemos \displaystyle\int_{t_b}^{t_c}\Bigl(\frac{26}{3}t^{-2}+40t\Bigr)\,dt.
De I_3 obtemos -\displaystyle\int_{t_c}^{t_d}\Bigl(\frac{1}{3}t^{-2}+\frac{1}{2}t\Bigr)\,dt
(o sinal menos porque esses termos vieram subtraídos em I_3).
Somando essas três contribuições (chamei isso de S) obtemos uma combinação de termos do tipo t^{-1} e t^{2} avaliados nos pontos t_a,t_b,t_c,t_d. Escrevendo S explicitamente:
\begin{aligned} S =&\; \Bigl(-9t^{-1}+\tfrac{81}{4}t^{2}\Bigr)\Big|_{t_a}^{t_b} +\Bigl(-\tfrac{26}{3}t^{-1}+20t^{2}\Bigr)\Big|_{t_b}^{t_c} +\Bigl(\tfrac{1}{3}t^{-1}-\tfrac{1}{4}t^{2}\Bigr)\Big|_{t_c}^{t_d}. \end{aligned}
Se expandirmos essas somas iremos obter termos em t_a^{-1},t_b^{-1},t_c^{-1},t_d^{-1} e em t_a^{2},t_b^{2},t_c^{2},t_d^{2} com coeficientes numéricos.
No entanto, agora lembramos das relações cúbicas que definiram t_a,t_b,t_c,t_d:
Essas relações permitem substituir potências cúbicas por expressões lineares nas avaliações nos extremos; e como os termos “difíceis” das integrais I_1 e I_3 também envolviam exatamente (1-2t) e (1-t) (isto é, os mesmos polinômios), quando juntarmos todas as contribuições as partes que dependem complicadamente das raízes cúbicas cancelam entre si, restando apenas termos racionais simples avaliados nos extremos que, após a aplicação das identidades cúbicas nos extremos, dão finalmente um valor numérico racional.
No passo seguinte executo essa simplificação de maneira direta: avalio as primitivas nas fronteiras, uso as relações cúbicas para substituir t^3 por expressões lineares onde necessário, e verifico que as “peças complicadas” cancelam exatamente, restando um número que é computável e que resulta em \dfrac{26}{3}.
Passo 6
Avaliação final — usar as relações dos pontos críticos para simplificar e somar tudo.
Vou mostrar os pontos-chave da simplificação (as contas algébricas completas dão exatamente o que descrevo, e o leitor pode conferir cada substituição):
Nas junções entre os intervalos (em t_b e t_c) as parcelas provenientes das avaliações das primitivas em w coincidem porque, por definição, nesses pontos os valores de w nas curvas coincidem (por exemplo em t_b vale t_b^{3}=1-2t_b). Isso faz com que os termos “difíceis” \dfrac{(1-2t)^{3}}{t^{11}}, \dfrac{(1-2t)^{4}}{t^{11}}, \dfrac{(1-t)^{3}}{t^{11}}, etc., que aparecem como avaliações nas fronteiras se cancelem quando somamos I_1+I_2+I_3.
Assim, a soma total I_1+I_2+I_3 reduz-se a uma combinação das primitivas simples avaliadas apenas nos extremos t=t_a e t=t_d (pois os valores nas fronteiras intermediárias t_b,t_c cancelam entre as parcelas adjacentes). Ou seja, depois do cancelamento interno, resta algo da forma
\Bigl(A\,t^{2}+B\,t^{-1}\Bigr)\Big|_{t_a}^{t_d},
com constantes A,B provenientes da soma das contribuições anteriores.
Agora aplicamos as relações cúbicas nos extremos t_a e t_d para eliminar as potências cúbicas que eventualmente apareçam e simplificar t_a^{2}, t_d^{2} e t_a^{-1},t_d^{-1} combinando com as equações cúbicas. Feito isso (procedimento algébrico direto, usando t^{3}=1-t ou t^{3}=1-2t conforme o ponto), os termos resultam em números racionais e a soma final dá exatamente
Resumo das manipulações: foi feita a transformação (x,y)\mapsto(s,t)=(x+y,x) (jacobiano de módulo 1), trocou-se s por w=s-t=y para escrever a integral como
particionou-se o intervalo de t em três partes de acordo com a ordem dos limites, integrou-se em w (obtendo uma primitiva polinomial) e depois integrou-se em t. Os termos complicados avaliados nas fronteiras internas cancelaram por causa das relações cúbicas que definem os pontos de partição, sobrando uma expressão que, avaliada nos limites finais t_a e t_d e simplificada com as mesmas relações cúbicas, produz o valor racional \dfrac{26}{3}.