problem in this integration

1 view (last 30 days)
Nasir Qazi
Nasir Qazi on 6 May 2012
% 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
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?

Sign in to comment.

Answers (1)

Alexander
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)

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!