Numerical integration Matlab plot

The function which I want to integrate is:

fun = @(x) (1./(exp(0.00004781./(x))-1).*1./x.^4);

q = integral(fun,0,x1)

I want to make a semilogx plot of (q, x1) in matlab for values of x1 from 100*(10^(-9)) to 100*(10^(-6)). How can I do it using for loop?

Thanks in advance.

Raj Patel.

Answer 1

[Ameer Hamza]

You can do it using arrayfun (similar to for-loop, but compact)

x1 = logspace(-9, -6, 10);

fun = @(x) (1./(exp(0.00004781./(x))-1).*1./x.^4);
y = arrayfun(@(x1_) integral(fun, 0, x1_), x1);

semilogx(x1, int_fun);

[Raj Patel]

Thanks Ameer for your help. Appreciate your effort.

[Ameer Hamza]

I am glad to be of help!