How to integrate a function u^3/(exp(u)-1) in the range of [0 0.033]??.

The code is as follows:

syms x

fun=x^3/(exp(x)-1)

defIntVal=int(fun,[0 0.033])

I have used integral command also in place of int

the defIntVal is giving a very complex sym type variable. How to get value from this variable???

[Bjorn Gustavsson]

Wouldn't a numerical integration be good enough?

[Raj Kishor]

That is what I am looking for. Please send me the final evaluated expression.

Answer 1

[Bjorn Gustavsson]

if numerical integration is good enough:

Val = integral(@(x) x.^3./(exp(x)-1),0,0.033)
% for which I get: 1.1831e-05

HTH