arrayfun applied to integral

Hello!

I would like to vectorialize a subroutine calculating an integral.

The main entry is an array that includes, in row, parameters of the function to integrate and the boundaries for the integration.

It works fine when called into a for loop, but I failed to vectorize it using arrayfun. Below is the code.

Any help to speed up this code would be most the welcome!!!

Many thanks

Olivier

function hV=change_potential_by_history_NoDOV3(t,t0,tpeak,w,past,slops,extrems)

%Matlab2022

parfor i=1:length(t)

int(i)=integrate_PSP(t(i),t0(i),tpeak(i),w(i),past(i));

end

%I WANT SOMETHING LIKE

%int=arrayfun(@integrate_PSP, t,t0,tpeak,w,past)

%

hV=sum

end

function int=integrate_PSP(t,t0,tpeak,w,past)

if t>t0

tr=t-t0;

if tr-past>0

int=integral(@(tr) w.*tr.*exp(-tr./tpeak),tr-past,tr)/past;

else

int=integral(@(tr) w.*tr.*exp(-tr./tpeak),0,tr)/tr;

end

else

int=0;

end

end

[Torsten]

Hint:

syms x a b
f = a*x*exp(-x/b)

f = 

F = int(f)

F = 

[Olivier]

Thank!