Integration by numerical methods

Hi, please someone can help how to compute the following integral numerically, since analytic integration is not possible for all values of "z". Here " z" lies in between "0" to "1".

Answer 1

[Alan Stevens]

Something like this:

z = 0.01:0.01:1;
I = zeros(1,numel(z));
for i = 1:numel(z)
I(i) = integral(@(x) sqrt(x)./(exp(x)./z(i) - 1),0,inf);
end
plot(z,I),grid

[P Rakesh Kumar Dora]

Can we get explicit form of the integral in terms of 'z' ?Like symbolic tool box gives explicit form.

[Alan Stevens]

Wolfram Alpha gives a result in terms of a function Li. MATLAB doesn't seem to have that function (unless it is there under another name!).

[P Rakesh Kumar Dora]

thanks