problem in this integration
1 view (last 30 days)
Show older comments
% Unable to integrate this , Please help
syms T T0
b = 0.141;cp0 = 38.12;h0 = 1684;R= 8.314;Act_W = 0.98;D_Vm = 3.0;DCP_0 = 1263; DH_0 = 1389;
cp = cp0 + b(T-T0);
h = h0 + int((cp,T,T0,T));
DCP_0/R*T + int((h/R*T^2),T,T0,T) + int((D_Vm/R*T),P,0,P) - log(Act_w)
1 Comment
Walter Roberson
on 6 May 2012
What problem are you observing? Are you getting an error message? Which of the three integrals are you having difficulty with?
Did you happen to notice that for int((D_Vm/R*T),P,0,P) you are integrating with respect to a variable that does not appear in the expression being integrated?
Answers (1)
Alexander
on 7 May 2012
Not sure what you want to do. Could it be, that you meant this:
syms T T0 P
b = 0.141;cp0 = 38.12;h0 = 1684;R= 8.314;Act_W = 0.98;D_Vm = 3.0;DCP_0 = 1263; DH_0 = 1389;
cp = cp0 + b*(T-T0);
h = h0 + int(cp,T,T0,T);
DCP_0/R*T + int((h/R*T^2),T,T0,T) + int((D_Vm/R*T),P,0,P) - log(Act_W)
I did these changes to your code:
- syms T T0 P instead of syms T T0
- b*(T-T0) insteat of b(T-T0)
- int(cp,T,T0,T) instead of int((cp,T,T0,T))
- log(Act_W) instead of log(Act_w)
0 Comments
See Also
Categories
Find more on Symbolic Math Toolbox in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!