Numerical integration in Matlab

Hi I want to numerically integrate some function in Matlab using "integral". Say f_int(x). However, f_int can not handle vector inputs. This seems to cause "integral" troubles? Any ideas how to quickly fix this, w/o having to rewrite f_int(x)? Thanks a lot. M

Answer 1

[Andrei Bobrov]

Try:

f2_int = @(y)arrayfun(@(x)f_int(x),y);

out = integral(f2_int,a,b); % here a and b - double

[Mike Hosea]

Andrei's answer is the way I prefer to do it. However, in a pinch you can leverage the 'Arrayfun',true option to treat the problem as a vector-valued integral (just one that happens to return a vector of only one element)

out = integral(f_int,a,b,'ArrayValued',true);