problem in this integration

2 次查看(过去 30 天)
Nasir Qazi
Nasir Qazi 2012-5-6
% 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 个评论
Walter Roberson
Walter Roberson 2012-5-6
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?

请先登录,再进行评论。

回答(1 个)

Alexander
Alexander 2012-5-7
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