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 commentaire

Walter Roberson
Walter Roberson le 6 Mai 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?

Connectez-vous pour commenter.

Réponses (1)

Alexander
Alexander le 7 Mai 2012

0 votes

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)

Community Treasure Hunt

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

Start Hunting!

Translated by