求助一个求积分的程序

suntht

要写程序得到如下表达式的结果，T随时间t的改变
程序总是报错，请大家指教
alpha=9.003e-5;k=220.0;

Pi=3.141592654;
beta=1.0e3;
b=1.0e-3;
B=1.53e-8;
I0=1.3e6;

A1=0.93*I0*B/k;

v=1.0;
x=1.0e-3;
y=b;
  % 固定x,改变t  
  
fid=fopen('Gn,x=1.0e-3.txt','w');

for i=linspace(1e-5,1e-2);
    t=i;
   
Tall=0.0;
  for j=linspace(1,7)
       n=j;
    gamma=n*Pi/b;
   lm=alpha*gamma^2;
   
   fn=B*b^2/(B^2*b^2+n*n*Pi*Pi)*(1.0-(-1)^n*exp(-B*b));
syms tau  
g1=exp(-beta*tau)/sqrt(t-tau)*exp(-lm*(t-tau))*exp(-(x-v*tau)^2/4.0/alpha/(t-tau));
gn=int(g1,tau,0,t);
gn1=vpa(gn,6);
gn2=eval(gn1);

T1=A1*alpha/b/sqrt(Pi*alpha)*cos(n*Pi/b*y)*fn*gn2;
Tall=Tall+T1;

end
     fprintf(fid,'%10.7f  %20.15f\n',t,Tall);
  end
  sta=fclose (fid);