function f=te5(t,x)
f=zeros(3,1);
f(1)=-0.04*x(1)+10^4*x(2)*x(3);
f(2)=0.04*x(1)-10^4*x(2)*x(3)-3*10^7*x(2)^2;
f(3)=3*10^7*x(2)^2;

