Run this:
ClearAll["Global`*"];
b[x_] := (a + e)^c[x];
v1 = Normal[Series[Integrate[b[x], x], {e, 0, 1}]]
v2 = Integrate[Normal[Series[b[x], {e, 0, 1}]], x]
to get this:
The first one is clearly wrong but the second one seems correct.
What am I missing or this is a bug in Mathematica?
I would not ask here but I need to perform some complicated transformations of integrals and this seems killing what I need to do. I am running 11.1.0.0 version of Mathematica. Does this work correctly in the most recent version?
Integrate[b[x], x]
expression to return? $\endgroup$Integrate[b[x], x]
to return, well,Integrate[b[x], x]
. The transformations above are "trivial" for a human. Mathematica seems to fail in a Taylor expansion there. $\endgroup$