I veru much thankful to you for your valuable inputs in the above code. its works fine for me.
Based on the above inputs I am trying to attempt one research problem. please find the attached files including base paper and problem file which i wish to solve. I attempted momentum eqn (eqn 9 in problem figure).
I dont have any errors in the code but have few issues.
In maple code in eqn 6 (f1 to f6 ) has 1/P^3.. terms are involed, that is denominator with 'p' powers.
for m from 0 to 5 do
f[m+1] := unapply(simplify(int(int(int(R[m](x), x), x), x)), x):
As I am integrating with respect to x, why we are getting p powers in the denominator. (Is it possible to implement with out using unapply command-Initially I wrote without using unapply, you added unapply.)
I also wish to bring to ur notice that fig profiles in the code are not satisfying the boundary conditions. I humbly request you that please fix the bugs which produce figure profile as shown in the basepaper. I wrote only for momentum, how to implement for both eqations (equation 9 and 10 subject to 11-see problem.pdf).
please help me sir, As this is my research problem. If you could help, I will be greatful to you.
looking forward your kind response.
Thank you in advance