How to do integration using MATLAB

how to get integral of (u^4*exp(u)/(exp(u)-1)^2) du with limit 0 to T (where T is temperature) in MATLAB. I want analytic solution in T where T value will be assigned later. I tried the following: >> u=sym('u')

u =

u

>> int((u^4*exp(u)/(exp(u)-1)^2)) Warning: Explicit integral could not be found.

ans =

int((u^4*exp(u))/(exp(u) - 1)^2, u)

>>

[Mahesha MG]

Is there any way in symbolic math? I tried the following. >> syms u theta t >> int(u^4*exp(u)/(exp(u)-1)^2,0,theta/t) Warning: Explicit integral could not be found.

ans =

int((u^4*exp(u))/(exp(u) - 1)^2, u = 0..theta/t)

Answer 1

[Pourya Alinezhad]

instead of using symbolic math you can use numerical integration .

function answer=myfunction(T)
   %* Labatto quadrature –
    q =quadl('(u.^4.*exp(u))./((exp(u)-1).^2)',0, T)
   %you can use quad() and trapz() functions instead
  end

[Mahesha MG]

Thank you... It is working.