ai =Beta-Ceta;
  bi =(Aeta/tau-Beta/h+2*Ceta/(h*h)+Deta);
  ci =Ceta/(h*h);
  fi =alw*T0-Aeta*Tn(i)/tau+ (T0*et0*Q*k0*ro*exp(-E/(R*T(i))))/T(i);

Rw=0.05; % где используется?
Tv=315;% где используется?

function [T,x] = funcgnew(t_end)
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
L=0.2;
% вместо A1...C4 необходимо поставить требуемые коэффициенты.
 A1=1;
 A3=1;
 B1=1;
A2=1;
B2=1;
C2=1;
A4=1;
B4=1;
C4=1;
 mf=500;
 
 
T=zeros(1,mf);
Ts=zeros(1,mf);
Tn=zeros(1,mf);
alfa=zeros(1,mf);
beta=zeros(1,mf);
 
N=500;% количество пространственных узлов
T0=1200;% - начальная температура
h=L/(N-1);% шаг по x
tau= t_end/100;%шаг по t
 
  for i=1: N 
    T(i) =T0;
  end;
  time =0;
while (time<t_end)
 
 time=time+tau;
alfa(1) =0;
beta(1)=1200;
   for i =1 :N 
   Tn(i)=T(i);
   end;
  while(1==1)
for i=1:N  
Ts(i)=T(i);
end
  for i=2:N-1 
 
  ai =A1/(h*h)-B1/h;
  bi =2*A2/(h*h)+B2/tau+C2/h;
  ci =A3/(h*h);
  fi =-A4*Tn(i)/tau-B4*exp(-1/(1*T(i)))/T(i)-T0-C4*T(i);
alfa(i) =ai/(bi-ci*alfa(i-1));
beta(i)=(ci*beta(i-1)-fi)/(bi-ci*alfa(i-1));
  end;
   T(N)=300;
         
  for i=N-1:-1:1  
T(i) =alfa(i)*T(i+1)+beta(i);
  end;
max=abs(T(1)-Ts(1));
  for i=2:N  
     if(max<abs(T(i)-Ts(i))) 
         max =abs(T(i)-Ts(i));
     end
  end
     if(max<=eps) 
         break;
     end
  
  end
  
end;
for i=1:N
    x(i)=h*(i-1);
end
T=T-273;
 
end

clear,clc
[T,x]=funcgnew(600);
plot(x,T)
grid on

function [T,x] = funcgnew(t_end)
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
L=0.2;
% вместо A1...C4 необходимо поставить требуемые коэффициенты.
alw=50; 
Rw=0.05; 
R1=0.3; 
R2=1500;
c2=660; 
la2=4; 
la1=0.084; 
al1=0.5; 
al2=0.5;
et0=0.077; 
R=8.314;
cp=1000; 
Tv=315;
E=126000; 
Q=1513000;
k0=50000000000;
u=0.001;
v=0.5;
u0=v/u;
Al0=(2*alw)/Rw;
Fi=(R2*c2)/(R1*cp);
k=(al1*la1+al2*la2)/(R2*c2+R1*cp);
A0=Al0/(R2*c2+R1*cp);
A1=(R1*k0*Q*et0)/(R2*c2+R1*cp);
A2=(v)/(1+Fi);
mf=500;
 
T=zeros(1,mf);
Ts=zeros(1,mf);
Tn=zeros(1,mf);
alfa=zeros(1,mf);
beta=zeros(1,mf);
N=500;% количество пространственных узлов
T0=1200;% - начальная температура
h=L/(N-1);% шаг по x
tau= t_end/100;%шаг по t
 
  for i=1: N 
    T(i) =T0;
  end;
  time =0;
while (time<t_end)
 
 time=time+tau;
alfa(1) =0;
beta(1)=1200;
   for i =1 :N 
   Tn(i)=T(i);
   end;
  while(1==1)
for i=1:N  
Ts(i)=T(i);
end
  for i=2:N-1 
   ai =k/(h*h);
  bi =1/tau+(2*k)/(h*h)+A2/h;
  ci =A2/h+(k)/(h*h);
  fi =-1/tau+A0*(T(i)-T0)-A1*exp(-E/R/T(i));
alfa(i) =ai/(bi-ci*alfa(i-1));
beta(i)=(ci*beta(i-1)-fi)/(bi-ci*alfa(i-1));
  end;
   T(N)=300;
         
  for i=N-1:-1:1  
T(i) =alfa(i)*T(i+1)+beta(i);
  end;
max=abs(T(1)-Ts(1));
  for i=2:N  
     if(max<abs(T(i)-Ts(i))) 
         max =abs(T(i)-Ts(i));
     end
  end
     if(max<=eps) 
         break;
     end
  end
end;
for i=1:N
    x(i)=h*(i-1);
end
T=T-273;
end