# How to solve this complicated double integral problem?

I have written this in way to make it as much as possible non-confusing. I will start describing my problem and I will walk you through my question,

I have a double integration which I am trying to solve of the following form,

$$F=\int\limits_{0<\gamma_1<\gamma_2<+\infty} \mathcal{L}\bigl(\sum_{i=1}^2\gamma_i^{-1} \bigl)\ g(\gamma_1,\gamma_2)\ d\gamma_1 d\gamma_2$$

where $$\mathcal{L}(t) = \text{exp}\left(- \int_{\gamma_2}^\infty (1- \frac{1}{1+t \ x^{-1}})\ h(x)\ dx\right)$$ Now assume that $h(x)$ is a very complex expression, I tell you that $\mathcal{L(t)}$ can not be integrated in closed form analytically and should be solved numerically.

So to solve my problem I have to take numerical values of $t$ in order to integrate numerically. But notice that to solve my main integration $F$ defined above, I need to set $$t=\sum_{i=1}^2\gamma_i^{-1}$$ and then integrate numerically to obtain $L\bigl(\sum_{i=1}^2\gamma_i^{-1} \bigl)$ by taking numerical values of $\gamma_1$ and $\gamma_2$.

Now that I took numerical values of $\gamma_1$ and $\gamma_2$ to solve my $\mathcal {L}(t)$. My inner integral is then numerical values for certain values of $\gamma_1$ and $\gamma_2$ that I picked.

My question is, since I was obliged to take certain values of $\gamma_1$ and $\gamma_2$ to solve for my non-integrable function $\mathcal{L}(t)$ - since it happened to be that by $t$ is a function of $\gamma_1$ and $\gamma_2$. How would I proceed to solve my function $F$? Am I now obliged to integrate numerically over the same values of $\gamma_1$ and $\gamma_2$ I picked to solve for $\mathcal{L}(t)$ ?

In general do you think this way of analyzing such a problem is correct. Any better ideas on tackling such a problem? for example approximations/

Thanks

• There are some inconsistencies with your notation. Please look it over and make edits as needed. It'll help. – Alexander Vlasev Nov 3 '14 at 9:08
• DO you mean the index $k$ ?I dont see any other – Tyrone Nov 3 '14 at 14:34

Since the inner integral included in $\mathcal{L(t)}$ cannot be determined in closed form and must be solved by numerical methods, the double outer integral included in $F$ has to be solved numerically as well. Calculating the inner integral allows you to obtain a tabulation of the values of $\mathcal{L(t)}$ for several values of $\gamma_1$ and $\gamma_2$. From this tabulation, a numerical integration technique can be used to calculate $F$ as well.
Methods commonly used in this context (e.g. the Newton-Cotes quadrature formulas) typically approximate the integral of a function $f(x)$ by estimating various degree polynomials, starting from numerical values of $x$ and $f(x)$ tabulated at regularly spaced intervals. As correctly noted in the OP, it is therefore necessary to integrate numerically over the same values of $\gamma_1$ and $\gamma_2$ picked to solve for $\mathcal{L}(t)$. In this view, it is convenient to pick $\gamma_1$ and $\gamma_2$ values by choosing regular intervals, whose dimension should be predefined according to the desired precision level.
To solve the problem of integrating on an infinite domain, you could consider the possibility of mapping the infinite interval to a finite interval using an appropriate substitution, i.e. choosing a function of $\gamma_1$ and $\gamma_2$ that goes to infinity in a finite range (for example $\tan\gamma_i$, $\cot \gamma_i$, or others).
• Thank you. The problem is the limit of my integration which are over $\gamma_1$ and $\gamma_2$ are from o to $\infty$. So I was trying to avoid numerical. Do you have an idea how one can approximate an integral that has limit of $\infty$ ? – Tyrone Nov 6 '14 at 22:49
• thank you ! Do you mean for example $\gamma_2$ would run from 0 to 90 if I take the example of $tan(\gamma_2)$? – Tyrone Nov 6 '14 at 23:09
• Yes, you are right. It would run from $0$ to $\displaystyle \frac{\pi}{2}$. – Anatoly Nov 6 '14 at 23:10
• OK so my strategy is the following, take numerical values of $\gamma_1$ and $\gamma_2$ and numerically integrate to find $\mathcal{L}$, this can easily be done using for example MATLAB function "integral" then I would have for each $\gamma_1$ and $\gamma_2$ a corresponding $\mathcal{L}$. However now that I have picked $\gamma_1$ and $\gamma_2$ I can't numericaly integrate to find $F$, I would have to do something like Reimann sum over $\gamma_1$ and $\gamma_2$, i.e I can't even use built in numerical integration methods I would have to do the Rieman sum, am I right? – Tyrone Nov 6 '14 at 23:17