这段程序的积分函数怎么改才能用对

load('q7.txt');
h=6.63*10^-34;
c=3*10^10;
k=1.38*10^-23;
t=296;
E=4576.2256;
v0=5000.495078;
N=0.66;
P=1;
X=0.3;
T=297;
GAMA0=0.0596;
GAMA1=0.063;
GAMA=GAMA0*(t/T)^N;
GAMA2=GAMA1*(t/T)^N;
ln2=log(2);
q=q7(:,2);
Q=q(T);
S=0.637*q(296)/Q*(t/T)*exp(-h*c*E/k*(1/T-1/t))*(1-exp(-h*c*v0/(k*T)))*(1-exp(-h*c*v0/(k*t)))^-1;
v=5000:0.01:5010;
for i=1:1001
vD=7.1623*10^-7*v0*sqrt(T/44);
vC=P*X*2*GAMA+P*(1-X)*2*GAMA2;
w=2*sqrt(ln2)*(v(i)-v0)/vD;
%y=2*v*sqrt(ln2)/vD;
a=sqrt(ln2)*vC/vD;
Faiv=@(y) 10*a./pi*(exp(-y.^2)/(a.^2+(w-y).^2));
FAIV(i)=integral(Faiv,0,5);
end