I am trying to solve the indefinite integral with user-defined rules. the integral expression and details are shown below
where the P is a function of a,b and c. and I would like Mathematica to follow what I did. I keep struggling for a very long time, but I could only make the very beginning part work, the code I had right now is
TT = \[Integral]\[Integral]\[Integral]D[P[x1, x2, x3, t]*x1, x1]*x1^r*
x2^s*x3^n \[DifferentialD]x1 \[DifferentialD]x2 \
\[DifferentialD]x3
ExpandAll[TT]
I tried using the ExpandAll to expand out the expression into 2 different parts, but it does not work.
any help from anyone will be very grateful.
Sincerely, Li

rule1 = \!\(\*SuperscriptBox[\(P\), TagBox[ RowBox[{"(", RowBox[{"1", ",", "0", ",", "0", ",", "0"}], ")"}], Derivative], MultilineFunction->None]\)[x1, x2, x3, t]*\[DifferentialD]x1 -> \[DifferentialD]p;, any help from anyone is very grateful – xiaofu li May 04 '20 at 22:09