MATLAB: Numerical integration Matlab plot

matlab semilogx 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.

Best Answer

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);
Related Question