problem in this integration

% 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

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)

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)

Asked:

on 6 May 2012

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!