El OP pedido una respuesta mediante el análisis numérico. En primer lugar quisiera decir que incluso Mathematica y MAPLE no se dan exactamente las mismas respuestas para esta integral. ARCE da:
F(x) := ln(exp(2*x)+1)/(1+sin(2*x));
evalf(int(F(x),x=0..Pi/2),16);
1.830481056481415
En comparación con el valor en el comentario por Iuʇǝƃɹɐʇoɹ : $1.830481056481482$ .
A continuación, OK, aquí es un muy a la fuerza bruta (y , por tanto, un pedazo de la torta) programa:
programa numérico;
la función F(x : haga doble) : doble;
comenzar
F := ln(exp(2*x)+1)/(1+sin(2*x));
end;
la función inferior(N : integer) : doble;
var
k : integer;
x,dx,suma : doble;
comenzar
suma := 0;
dx := (Pi/2)/N;
para k := 0 hasta N-1 hacer
comenzar
x := (Pi/2)*k/N+dx/2;
{ Regla del punto medio }
suma := suma + F(x)*dx;
end;
inferior := suma;
end;
la función superior(N : integer) : doble;
var
k : integer;
x1,x2,dx,suma : doble;
comenzar
suma := 0;
dx := (Pi/2)/N;
x2 := 0;
para k := 1 hasta N hacer
comenzar
x1 := x2;
x2 := (Pi/2)*k/N;
{ Trapecio regla }
suma := suma + (F(x1)+F(x2))/2*dx;
end;
superior := suma;
end;
comenzar
Writeln(inferior(1000000));
Writeln(superior(1000000));
final.
Salida:
1.83048105648054 E+0000
1.83048105648321 E+0000
Ayuda a hacer un pequeño boceto de la función $F(x) = \ln(e^{2x}+1)/(1+\sin(2x))$ :
Debido a esta función de la conducta, la
Regla del punto medio
dar un límite inferior, mientras que el
La regla del trapecio
da un límite superior de la integral. Así, nos encontramos con la certeza de que:
$$
1.830481056480 <
\int_{0}^{\pi/2} \frac{\ln\left(e^{2x} + 1\right)}{1 + \sin2x}\mathrm dx
< 1.830481056484
$$
Los resultados con Mathematica y MAPLE están bien dentro de estos límites.