Почему улетает в бесконечность

[Maik2013]

Здравствуйте, подскажите пожалуйста где может быт тут ошибку

$$ \dfrac{d^2T}{dx^2}=ua_{1}\dfrac{dT}{dx}+a_{2}(T-T_{0})-\dfrac{a_{3}}{T}\exp\left(-\dfrac{E}{RT}\right) $$

В М-файле с именем Rk.m:

Используется синтаксис Matlab M

function dT=RK(x,T)
alw=50; 
Rw=0.3; 
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; 
T0=300; 
E=126000; 
Q=15130000;
k0=50000000000;
u=0.0001;
v=1;
u0=v/u;
Al0=(2*alw)/Rw;
a1=(R1*(u0+1)*cp+R2*c2*u0)/(al1*la1+al2*la2);
a2=Al0/(al1*la1+al2*la2);
a3=(R1*Q*et0*k0)/(al1*la1+al2*la2);
dT(1,1) = T(2); % первая производная dT/dx
dT(2,1)=u*a1*T(2) + a2*(T(1)-T0) - (a3/T(1))*exp(-E/R/T(1));
end
 

и запуск расчетов:

Используется синтаксис Matlab M

[T,y]=ode45(@RK,[1 10], [300, 1200]);
plot(T,y)