Error using integral (line 85). A and B must be floating-point scalars.

I=.3*pi:.1*pi:.8*pi;
plot(I,sin(I)./I+integral(@(t) (1-cos(t))./t, 0, I)

Error using integral (line 85)
A and B must be floating-point scalars.

Answer 1

[Walter Roberson]

integral() is not designed to integrate to a list of different values.

You can use arrayfun() to integrate to each of the values in turn, having it do all the previous work all over again. Or you can break it up into intervals and integrate over each of them and use cumsum() to get the list of values.

targets = [0 I];
cumsum( arrayfun( @(IDX) integral(@(t) (1-cos(t))./t, targets(IDX), targets(IDX+1)), 1:length(targets)-1 ) )

[Oleg Primachenko]

Thank!