function dTdt
clear all;clc
format long

[t T]=ode15s(@odes,[0 1e-12],[300 300 1])
figure(1)
plot(t,T(:,1),'b',t,T(:,2),'r',t,T(:,3),'k'),axis([0 max(t) min(min(T))-10 300+10]);




function f=odes(t,y)

Te=y(1);
T1=y(2);
T2=y(3);
Ce=88*Te;
C1=2.7e16;

J=44;
R=0.6;
b=12.5e-9;
tp=0.1e-12;
G=2.7e16;

Q=-sqrt(4*log(2)/pi)*(1-R)*J/b/tp*exp(  -4*log(2)*((t-tp)/tp)^2  );

f1=-G*(Te-T1)+Q;
f2=G*(Te-T1);
K=20e12;
TC=800;
f3=K*T2*T1/TC*(1-T2*coth(T2*TC/Te));

f=[f1/Ce;f2/C1;f3];