integration of number array with variable

T(ez) is a number array. is there anyone who can help me to do this integration?

[Torsten]

Use trapz.

[Abu Nazir]

I do have calculated 600 value of T(Ez) corresponding to Ez number array. EF1, EC1, EF2 and EC2 600 numerical value which I have.

I want to vary T - temperature and calculate the J.

Answer 1

[Torsten]

Ez=...;
TEz = ...;

q=...;
mstar=...;
h=...;
kB=...;
    
T=...;
J = zeros(size(T));
for i = 1:numel(T)
    Y = TEz.*log((1+exp((EF1-Ez-EC1)/(kB*T(i))))./(1+exp((EF2-Ez-EC2)/(kB*T(i)))));
    J(i) = q*mstar/(2*pi^2*h^3)*kB*T(i)*trapz(Ez,Y);
end

plot(T,J)