Is there an integral that proves $\pi > 333/106$?

Question

The following integral,

$$ \int_0^1 \frac{x^4(1-x)^4}{x^2 + 1} \mathrm{d}x = \frac{22}{7} - \pi $$

is clearly positive, which proves that $\pi < 22/7$.

Is there a similar integral which proves $\pi > 333/106$?

Answer 1

This integral would do the job:

$$\int_0^1 \frac{x^5(1-x)^6(197+462x^2)}{530(1+x^2)}\:dx= \pi -\frac{333}{106}$$

• Also you can refer to S.K. Lucas Integral proofs that $355/113 > \pi$, Gazette, Aust. Math. Soc. 32 (2005), 263-266.

• This is the link. (Thanks to lhf for pointing out.)

Answer 2

Although this is not exactly an answer to the question, it seems sufficiently related to mention: there are some direct generalizations, given on the Wikipedia page about this integral. For instance, $$0 < \frac14\int_0^1\frac{x^8(1-x)^8}{1+x^2}\ dx=\pi -\frac{47171}{15015}$$

In general, $$\frac1{2^{2n-1}}\int_0^1 x^{4n}(1-x)^{4n}\ dx <\frac1{2^{2n-2}}\int_0^1\frac{x^{4n}(1-x)^{4n}}{1+x^2}\ dx <\frac1{2^{2n-2}}\int_0^1 x^{4n}(1-x)^{4n}\ dx$$

which for $n=1$ (the integral in the question) gives slightly better bounds than just $\pi < 22/7$: $$ \frac{1}{1260} < \frac{22}{7} - \pi < \frac{1}{630}$$

Answer 3

In the beginning of 2009 I was posting re similar issue at several sites, namely, at sci.math.symbolic, www.math.utexas.edu, etc.

To repeat: In Paper 1 Lucas found, by brute-force search using Maple programming, several different variants of integral identities which relate each of several first Pi convergents (described in terms of OEIS <oeis.org> sequences as A002485(n)/A002486(n)) to Pi.

Further, in my above-mentioned postings, I conjectured the following identity below, which represents a generalization of Stephen Lucas' experimentally obtained identities between Pi and its convergents:

$$(-1)^n\cdot(\pi - \text{A002485}(n)/\text{A002486}(n))$$

$$=(|i|\cdot2^j)^{-1} \int_0^1 \big(x^l(1-x)^m(k+(i+k)x^2)\big)/(1+x^2)\; dx$$

where integer n = 0,1,2,3,... serves as index for terms in OEIS A002485(n) and A002486(n), and {i, j, k, l, m} are some integers (to be found experimentally or otherwise), which are probably some functions of n.

The "interesting" (I think) part of my generalization conjecture is that "i" is present in both:

denominator of the coefficient in front of the integral and in the body of the integral itself

For example for $\frac{22}{7}$

$$\frac{22}{7} - \pi = \int_{0}^{1}\frac{x^4(1-x)^4}{1+x^2}\,\mathrm{d}x$$

with $n=3, i=-1, j=0, k=1, l=4, m=4$ - with regards to my above suggested generalization.

In Maple notation

i:=-1; j:=0; k:=1; l:=4; m:=4;Int(x^l*(1-x)^m*(k+(k+i)x^2)/((1+x^2)(abs(i)*2^j)),x= 0...1)

yields 22/7 - Pi

It also works for found by Lucas

http://www.math.jmu.edu/~lucassk/Papers/more%20on%20pi.pdf

formula for $\frac{333}{106}$

$$\pi - \frac{333}{106} = \frac{1}{530}\int_{0}^{1}\frac{x^5(1-x)^6(197+462x^2)}{1+x^2}\,\mathrm{d}x$$

with $n=4, i=265, j=1, k=197, l=5, m=6$ -with regards to my above suggested generalization.

In Maple notation i:=265; j:=1; k:=197; l:=5; m:=6;Int(x^l*(1-x)^m*(k+(k+i)x^2)/((1+x^2)(abs(i)*2^j)),x= 0...1)

yields Pi - 333/106

And it works for Lucas's formula for $\frac{355}{113}$

$$\frac{355}{113} - \pi = \frac{1}{3164}\int_{0}^{1}\frac{(x^8(1-x)^8(25+816x^2)}{(1+x^2)}$$

with $n=5, i=791, j=2, k=25, l=8, m=8$ -with regards to my above suggested generalization.

In Maple notation

i:=791; j:=2; k:=25; l:=8; m:=8;Int(x^m*(1-x)^l*(k+(k+i)x^2)/((1+x^2)(abs(i)*2^j)),x= 0...1)

yields 355/113 - Pi

And it works as well for Lucas's formula for $\frac{103993}{33102}$

$$\pi - \frac{103993}{33102} = \frac{1}{755216}\int_{0}^{1}\frac{x^{14}(1-x)^{12}(124360+77159x^2)}{1+x^2}\,\mathrm{d}x$$

with $n=6, i= -47201, j=4, k=124360, l=14, m=12$ -with regards to my above suggested generalization.

In Maple notation

i:=-47201; j:=4; k:=124360; l:=14; m:=12;Int(x^l*(1-x)^m*(k+(k+i)x^2)/((1+x^2)(abs(i)*2^j)),x= 0...1)

yields Pi - 103993/33102

And also it works Lucas's formula for $\frac{104348}{33215}$

$$\frac{104348}{33215} - \pi = \frac{1}{38544}\int_{0}^{1}\frac{x^{12}(1-x)^{12}(1349-1060x^2)}{1+x^2}\,\mathrm{d}x$$

with $n=7, i= -2409, j=4, k=1349, l=12, m=12$ - with regards to my above suggested generalization.

In Maple notation

i:=-2409; j:=4; k:=1349; l:=12; m:=12;Int(x^l*(1-x)^m*(k+(k+i)x^2)/((1+x^2)(abs(i)*2^j)),x= 0...1)

yields 104348/33215 - Pi

And it works as well for $\frac{618669248999119}{196928538206400}$

which, by the way, is not part of A002485/A002486 OEIS sequences:

$$\frac{618669248999119}{196928538206400} - \pi = \frac{1}{755216}\int_{0}^{1}\frac{x^{14}(1-x)^{12}(77159+124360x^2)}{1+x^2}\,\mathrm{d}x$$

with $i= 47201, j=4, k=77159, l=14, m=12$ -with regards to my above suggested generalization.

In Maple notation

i:=47201; j:=4; k:=77159; l:=14; m:=12;Int(x^l*(1-x)^m*(k+(k+i)x^2)/((1+x^2)(abs(i)*2^j)),x= 0...1)

618669248999119/196928538206400 - Pi

I do not have computer math resources (Mathematica, Maple, etc.) to experimentally prove or disprove it for all larger n (but see my comment below).

Best Regards, Alexander R. Povolotsky

UPDATE #1:

Matt B. in his answer to my question on Mathematics Stack Exchange has provided analytical proof and improved my parametric formula by reducing the number of parameters from 5 to 4 (see Seeking proof for the formula relating Pi with its co seenvergents).

$$(-1)^n (\pi- \frac{p_n}{q_n}) = \int_0^1 \frac{x^{\epsilon+2m'}(1-x)^{2m'}(\alpha + \beta x^2) }{(\alpha - \beta) 2 ^{m'-2} (-1)^{\epsilon}(1+x^2)}dx$$.

Below is the list of parameters in Matt B.'s formula for all cases, covered in Stephen Lucas' publications - the stuff to the right of the arrow sign is the actual Maple code, which one could copy (while in the "edit" mode) and then paste into (let say) Inverse Symbolic Calculator (it accepts Maple code) and run it there.

NB I replaced for brevity some parameter names used by Matt B: "alpha" by "a", "beta" by "b", "epsilon" by "c", "m' " by "p".

104348/33215 - Pi -> a:=1349;b:=-1060;p:=6;c:=0;Int((x^(c+2p)(1-x)^(2p)(a+bx^2))/((a-b)2^(p-2)((-1)^(c)(1+x^2))),x=0...1)

Pi - 103993/33102 -> a:=124360;b:=77159;p:=6;c:=2;Int((x^(c+2p)(1-x)^(2p)(a+bx^2))/((a-b)(2^(p-2))((-1)^(c)(1+x^2))),x=0...1)

355/113 - Pi -> a:=25;b:=816;p:=4;c:=0;Int((x^(c+2p)(1-x)^(2p)(a+bx^2))/((a-b)(2^(p-2))((-1)^(c)(1+x^2))),x=0...1)

Pi - 333/106 -> a:=197;b:=462;p:=3;c:=-1;Int((x^(c+2p)(1-x)^(2p)(a+bx^2))/((a-b)(2^(p-2))((-1)^(c)*(1+x^2))),x=0...1)

22/7 - Pi -> a:=1;b:=0;p:=2;c:=0;Int((x^(c+2p)(1-x)^(2p)(a+bx^2))/((a-b)(2^(p-2))((-1)^(c)(1+x^2))),x=0...1)

Obviously parameters in the formula somehow depend on "n". The most straight forward dependency on "n" is observed in what Matt B named as "m'" (and I call it "p") : {2,3,4,4,6,6} ... Note that when "n" -> infinity, then the integral should come to 0 ...

UPDATE #2:

Analysis of Thomas Baruchel calculations results (see his answer to my question in

Seeking proof for the formula relating Pi with its co seenvergents)

led to observation that in originally supplied five parameter notation (i,j,k,l,m)

j = m/2 - 2

and correspondingly

m=2*(j+2)

This makes the original conjecture to depend on 4 parameters and to look like:

$$ (-1)^n\cdot(\pi - \text{A002485}(n)/\text{A002486}(n)) =(|i|\cdot2^j)^{-1} \int_0^1 \big(x^l(1-x)^{2(j+2)}(k+(i+k)x^2)\big)/(1+x^2)\; dx $$

This observation confirms previous Matt B result (see previous update), which also is based on 4 parameters.

Based on his calculations results, Thomas Baruchel also found that even with 4 parameters, this formula yields infinite number of solutions for each n.

Thomas shared with me his calculations results and supplied me with quite a few of valid combinations of i,j,k,l values - so now I have a lot of experimentally found five-tuples {n,i,j,k,l}, which satisfy above parametrization, where n varies in the range from 2 to 26.

Based on this data, of course, it would be nice to find how (if at all) i,j,k,l are inter-related between each other and with "n" - but such inter-relation (if exists) is not obvious and difficult to derive just by observation ... (though it is clearly seen that an absolute value of "i" is strongly increasing as "n" is growing from 2 to 26).

Looking at all available {i,j,k,l} solutions (initial ones found by me and those which were found by Thomas Baruchel program) one could observe that in all of them abs(l-j)=2*m where "m" is some positive integer.

If I didn't make a mistake RHS could be reduced (after performing integration) to:

(abs(i)2^j)^(-1)Gamma(2j+5)((k+i)Gamma(l+3)HypergeometricPFQ(1,l/2+3/2,l/2+2;j+l/2+4,j+l/2+9/2;-1)/Gamma(2j+l+8)+kGamma(l+1)HypergeometricPFQ(1,l/2+1/2,l/2+1;j+l/2+3, j+l/2+7/2;-1)/Gamma(2j+l+6))

May be from discussed parametric identity one could derive irrationality measure for pi, if to assume that RHS in this identity holds true, when the rational fraction on the LHS is equal to 0, then we have:

Pi = (abs(i)2^j)^(-1)Gamma(2j+5)((k+i)Gamma(l+3)HypergeometricPFQ(1,l/2+3/2,l/2+2;j+l/2+4,j+l/2+9/2;-1)/Gamma(2j+l+8)+kGamma(l+1)HypergeometricPFQ(1,l/2+1/2,l/2+1;j+l/2+3, j+l/2+7/2;-1)/Gamma(2j+l+6))

Perhaps someone could programmatically check if there are any {i,j,k,l}, which would satisfy above?

Update #3:

Thanks to Jaume Oliver Lafont, at least one case, answering affirmatively to the last question, is identified: i=-1, j=-2, k=1, l=0

$$\pi = \int_{0}^{1}\frac{4}{1+x^2}\,\mathrm{d}x$$

Should there be an infinite number of such cases?

see Randall's answer at math.stackexchange.com/a/2198869/28343

Update #4

David Trimas looked into the option involving partial solution for my conjecture by setting j=l=0 and deriving formulas for "i" and "k" for j=l=0 condition.

The result is following:

(-1) ^ n * (Pi − A002485(n)/A002486(n)) = ((Abs(i)) * 2 ^ j) ^ (-1) * Int((x ^ l * (1 - x) ^ (2 * (j + 2)) * (k + (i + k) * x ^ 2 ))/(1 + x ^ 2 ), x=0 ...1)

holds true for any n>2 and j=l=0 when

i =(-1)^(n) * 3 * A002486(n)

k = (-1)^(n) * (47 * A002486(n) - 15 * A002485(n)

For example for n=3 where A002485(3)=22 and A002486(3)=7

i=(-1)^337=-21

k=(-1)^3*(477 - 1522)=1

and the quick check via Inverse Symbolic Calculator using Maple (also could be done using Mathematica) provides the confirmation.

i:=-21;j:=0;k:=1;l:=0;(abs(i)2^j)^(-1)int((x^l(1-x)^(2(j+2))*(k+(i+k)*x^2))/(1+x^2),x=0...1) = 22/7 - Pi

Substitution of

i =(-1)^(n) * 3 * A002486(n)

and

k = (-1)^(n) * (47 * A002486(n) - 15 * A002485(n)

into general conjectured formula with the condition that j=l=0 (under which above formulas for "i" and "k" were derived) confirmed those formulas for "i" and "k" validity by bringing lhs and rhs to be equal.

Though above finding is obviously circular in nature (it doesn't reveal direct dependency of "i" and and "k" on "n" but rather does it via A002485(n) and A002486(n) ) - it is in my view a step forward.

Update #5

As I mentioned above in Update #2, observation of all obtained solutions for my conjectured identity (which represents generalization of Stephen Lucas' experimentally obtained identities between Pi and its convergents) shows that abs(l-j)=2*m where "m" is some positive integer. Going further I make another conjecture re identity which relates Log(2) and its convergents and is analogous to the one which relates Pi and its convergents. Below is conjectured by me formula (expressed in Maple notations) for relating Log(2) (that is Ln(2)) with ALL of its convergents - those, which are described via A079942(n)/A079943(n) ratio where A079942(n) and A079943(n) are OEIS integer sequences.

(-1)^n*(Log(2) − A079942(n)/A079943(n))=(Abs(i)2^j)^(-1)Int((x^l(1-x)^(2(j+2))*(k+(i+k)*x^2))/(1+x^2),x=0...1)

where integer n>2 serves as the index for the terms in the OEIS A079942(n) and A079943(n) integer sequences, and {i,j,k,l} are some signed integer parameters (which are some implicit functions of “n” and to be found experimentally or otherwise for each value of “n”) , and abs(l - j) = 2*m + 1, where “m" is some positive integer.

Answer 4

Another solution is given by the integral

$$ 0 < \int_0^1 \frac{x^4(1-x)^8}{4(1+x^2)}dx = \pi -\frac{2419}{770} = \pi - \frac{333}{106}-\frac{1}{20405} $$

This proves the stricter condition $$ \pi > \frac{333}{106}+\frac{1}{20405} $$

which implies $$ \pi > \frac{333}{106} $$

Similarly, for the fourth convergent (formula (6) http://www.math.ucla.edu/~vsv/resource/general/Lucas.pdf)

$$0<\int_0^1 \frac{x^{10}(1-x)^8}{4(1+x^2)}dx=\frac{3849155}{1225224}-\pi=\frac{355}{113}-\frac{5}{138450312}-\pi$$

Therefore $$\pi<\frac{355}{113}-\frac{5}{138450312}$$ and $$\pi<\frac{355}{113}$$

Even for the first convergent $$0<2\int_0^1 \frac{x(1-x)^2}{(1+x^2)}dx=\pi-3$$

so $$\pi>3$$

(See https://math.stackexchange.com/a/1618454/134791 for a proof for $3<\pi<4$ that uses this integral)

A series proof that $\pi>\frac{333}{106}$ is given by

$$\frac{48}{371} \sum_{k=0}^\infty \frac{118720 k^2+762311 k+1409424}{(4 k+9) (4 k+11) (4 k+13) (4 k+15) (4 k+17) (4 k+19) (4 k+21) (4 k+23)} \\=\pi-\frac{333}{106}$$

Answer 5

From the relationship

$$\frac{333}{106}=\frac{377-2·22}{120-2·7}$$

and integrals $$\pi = \frac{22}{7} - \int_0^1 \frac{x^4(1-x)^4}{1+x^2}dx$$

and $$\pi = \frac{377}{120} - \frac{1}{2}\int_0^1 \frac{x^5(1-x)^6}{1+x^2}dx$$

we obtain

$$\int_0^1 \frac{14x^4(1-x)^4-60x^5(1-x)^6}{1+x^2}dx = \int_0^1 \frac{2x^4(1-x)^4(7-30x(1-x)^2)}{1+x^2}dx = \int_0^1 \frac{2x^4(1-x)^4(7-30x+60x^2-30x^3)}{1+x^2}dx = 106\pi-333$$

Therefore,

$$\frac{1}{53} \int_0^1 \frac{x^4(1-x)^4(7-30x+60x^2-30x^3)}{1+x^2}dx = \pi-\frac{333}{106}$$

The numerator can be shown to be nonnegative for $0\leq x\leq 1$.

In fact, we can get a smaller numerator that is still nonnegative, which leads to a closer approximation to $\pi$

$$\pi=\frac{21991}{7000} +\frac{1}{50} \int_0^1 \frac{x^4 (1 - x)^4 \left(4 - 27 x (1 - x)^2\right)}{\left(1 + x^2\right)} dx$$

Here the numerator has been adjusted to have a double zero in $(0,1)$ while not crossing the axis WA link

The resulting fraction is

$$\frac{21991}{7000} = \frac{22}{7}-\frac{9}{7000}$$

which allows for writing the double inequality

$$\frac{22}{7}-\frac{9}{7000}<\pi<\frac{22}{7}$$

or, equivalently, an upper bound for the error in Archimedes' approximation

$$\frac{22}{7}-\pi<\frac{9}{7000}$$

Answer 6

Let us consider the polynomial

$$P_n(x):=1-x^2+x^4-x^6+\cdots x^{2n-2}=\frac{x^{2n}+1}{x^2+1}.$$

We have

$$0<\int_0^1\left(\frac{P_n(x)}{x^2+1}-\frac1{x^2+1}\right)^2dx<\int_0^1\left(\frac{x^{2n}}{0+1}\right)^2dx=\frac1{4n+1},$$ and the integral can be made as small as desired.

On another hand, the remainder of the long division of $(P_n(x)-1)^2$ by $(x^2+1)^2$ is a binomial $ax^2+b$, with $a,b$ integer. Then

$$\int_0^1\left(\frac{P_n(x)}{x^2+1}-\frac1{x^2+1}\right)^2dx=\int_0^1\left(Q(x)+\frac{ax^2+b}{(x^2+1)^2}\right)dx.$$

As

$$\int_0^1\frac{ax^2+b}{(x^2+1)^2}dx=\frac{b-a}4+\frac{b+a}8\pi$$

we can get arbitrarily close rational approximations by a rational integral.

---

For example, with $P_2:=1-x^2+x^4$ we have

$$\frac{(x^4-x^2)^2}{(x^2+1)^2}=x^4-4x^2+8-\frac{12x^2+8}{(x^2+1)^2}$$ so that by integration

$$0<\frac15-\frac43+8-\frac{5\pi}2+1<\frac19$$ or

$$\frac{698}{225}<\pi<\frac{236}{75}.$$

Answer 7

Without restrictions on the integrals, the question is virtually meaningless, because you can take any known bounds on $\pi$ (such as terms of the Gregory series), say $a<\pi<b$ and write

$$\int_0^a dx<\pi<\int_0^bdx.$$

Or if the function must be non-trivial, use any convergent squeezing of $\dfrac4{x^2+1}$ and

$$\int_0^1f(x)\,dx<\int_0^1\frac{4\,dx}{x^2+1}=\pi<\int_0^1g(x)\,dx$$ that you can make as tight as wanted.

Or better, any convergent squeezing of $\pi$, and

$$\int_0^1f(x)\,dx<\pi<\int_0^1g(x)\,dx.$$

In my other answer, I showed constructively that that there are bracketings of $\pi$ with integrals of rational functions as tight as you want, even though they are pretty inefficient.

Answer 8

A short Maple script shows that among the Dalzell-type integrals, $$\alpha \int_0^1 \frac{x^i (1 - x)^j}{1 + x^2} \,dx,$$ the integrals minimal in the sense that there is no other example $(i', j')$ for which $i' \leq i$ and $j' \leq j$, along with their lower bounds for $\pi$ are: $$\begin{array}{crl} \hline \displaystyle\frac{1}{1024} \int_0^1 \frac{(1 - x)^{24}}{1 + x^2} \,dx & \frac{17\,941\,970\,723}{5\,711\,177\,472}\!\!\!\! &= \frac{333}{106} + \frac{13\,399\,231}{302\,692\,406\,016} \\ \displaystyle\frac{1}{128} \int_0^1 \frac{x (1 - x)^{18}}{1 + x^2} \,dx & \frac{13\,685\,795}{4\,356\,352}\!\!\!\! &= \frac{333}{106} + \frac{14\,527}{230\,886\,656} \\ \displaystyle\frac{1}{16} \int_0^1 \frac{x^2 (1 - x)^{12}}{1 + x^2} \,dx & \frac{2\,264\,177}{720\,720}\!\!\!\! &= \frac{333}{106} + \frac{1501}{38\,198\,160} \\ \displaystyle\frac{1}{4} \int_0^1 \frac{x^4 (1 - x)^8}{1 + x^2} \,dx & \frac{2\,419}{770}\!\!\!\! &= \frac{333}{106} + \frac{1}{20\,405} \\ \displaystyle\frac{1}{2} \int_0^1 \frac{x^7 (1 - x)^6}{1 + x^2} \,dx & \frac{174\,169}{55\,440}\!\!\!\! &= \frac{333}{106} + \frac{197}{2\,938\,320} \\ \displaystyle \int_0^1 \frac{x^{10} (1 - x)^4}{1 + x^2} \,dx & \frac{141\,511}{45\,045}\!\!\!\! &= \frac{333}{106} + \frac{181}{4\,774\,770} \\ 2\displaystyle \int_0^1 \frac{x^{29} (1 - x)^2}{1 + x^2} \,dx & \frac{457\,304\,942\,543}{145\,568\,097\,675}\!\!\!\! &= \frac{333}{106} + \frac{147\,383\,783}{15\,430\,218\,353\,550} \\ 4\displaystyle \int_0^1 \frac{x^{24\,036}}{1 + x^2} \,dx & \ast & \approx \frac{333}{106} + 1.11071 \cdot 10^{-8} \\ \hline \end{array}$$ The $\ast$ denotes a rational number that for space reasons cannot be conveniently typeset as a ratio of integers, but its value is $\pi - \Psi(\frac{24039}{4}) + \Psi(\frac{24037}{4})$, where $\Psi(x) = \frac{d}{dx} \log \Gamma(x) = \frac{\Gamma'(x)}{\Gamma(x)}$ is the digamma function. The fourth of these integrals (namely, $(i, j) = (4, 8)$) is the example Jaume Oliver Lafont mentions in their answer.

Notice in the above table that $\alpha = 2^{2 - \frac{j}{2}}$. All of these integrals can be used to give tighter bounds by generalizing the observation of Dalzell that ShreevatsaR mentioned in their answer: In general we have $$2^{2 - \frac{j}{2}} \int_0^1 \frac{x^i (1 - x)^j}{1 + x^2} \,dx \geq 2^{2 - \frac{j}{2}} \cdot \frac{1}{2} \int_0^1 x^i (1 - x)^j \,dx = 2^{1 - \frac{j}{2}} \mathrm{B}(i + 1, j + 1),$$ where $\mathrm{B}$ is the beta function, so for any rational approximant $q < \pi$ derived as above, in fact we have the improve bound $\pi > q + 2^{1 - \frac{j}{2}} \mathrm{B}(i + 1, j + 1)$. Accounting for this observation gives improvements for just the last two integrals in the table, replacing them with the $2\displaystyle \int_0^1 \frac{x^{\color{#bf0000}{17}} (1 - x)^2}{1 + x^2} \,dx$ and $4\displaystyle \int_0^1 \frac{x^{\color{#bf0000}{156}}}{1 + x^2} \,dx$, respectively.

D.P. Dalzell, "On $22/7$." J. London Math. Soc., 19 (75 Part 3): 133$-$134. doi:10.1112/jlms/19.75_part_3.133, MR 0013425, Zbl 0060.15306.

Answer 9

Here are a few more examples of the type user9413 gave, that is, integrals $$\alpha \int_0^1 \frac{x^i (1 - x)^j (a x^2 + b) \,dx}{1 + x^2}$$ with $i, j \in \Bbb Z_{\geq 0}$, $a > 0$, $b > 0$, and value $\pi - \frac{333}{106}.$ Some quick Maple computations show that the integrals of this type minimal in the sense that there is no other example $(i',j')$ for which $i' \leq i$ and $j' \leq j$, are $$\begin{array}{rl} \hline \displaystyle \frac{1}{2938320} \!\!\!\!& \displaystyle \int_0^1 \frac{(1 - x)^{12} (1\,501 x^2 + 185\,146) \,dx}{1 + x^2} \\ \displaystyle \frac{1}{1484} \!\!\!\!& \displaystyle \int_0^1 \frac{x^2 (1 - x)^8 (36 x^2 + 407) \,dx}{1 + x^2} \\ \displaystyle \frac{1}{530} \!\!\!\!& \displaystyle \int_0^1 \frac{x^5 (1 - x)^6 (197 x^2 + 462) \,dx}{1 + x^2} \\ \displaystyle \frac{1}{742} \!\!\!\!& \displaystyle \int_0^1 \frac{x^8 (1 - x)^4 (181 x^2 + 923) \,dx}{1 + x^2} \\ \displaystyle \frac{1}{2\,533\,697\,595} \!\!\!\!& \displaystyle \int_0^1 \frac{x^{27} (1 - x)^2 (294767566 x^2 + 5362162756) \,dx}{1 + x^2} \\ \hline \end{array}$$ The third entry ($(i, j) = (5, 6)$) is the example user9413 gave. The second entry ($(i, j) = (2, 8)$) has the numerator with the smallest possible degree ($10$) among integrals of this type.