Matlab - Matlab ecuaciones diferenciales

function ode45ecuaciones4

 %Programa de ecuacion diferencial

 options = odeset('RelTol',1e-7,'AbsTol',[1e-7 1e-7 1e-7 1e-7 1e-7  ]);

 [T,X]=ode45(@Eq,[0 1e-17],[1 0 0 0 0],options);

 figure(gcf)

 plot(T,X(:,1),  T,X(:,2),  T,X(:,3) ,  T,X(:,4) , T,X(:,5),'linewidth',2  )

 axis([0  1e-18 abs(min(min(X(:,1:2))))  0.55  ])

 legend('x(1)','x(2)','x(3)','x(4)','x(5)')

 grid on

end

 function dx = Eq(t,x)

  R=8.314; %ctte de los fases

 To= 298; %temperatura inicial Kelvin

 h=5; %velocidad de calentamiento Kelvin/s

t=100;

 %Constantes de velocidad para la descomposicion de biomasa

 K1=(4.38E9)*exp(152700/(R*(To+(h*t))));%Descomposicion de la biomasa en gas

 K2=(1.08E10)*exp(148000/(R*(To+(h*t))));%Descomposicion de la biomasa en bio-oil

 K3=(3.375E10)*exp(111700/(R*(To+(h*t))));%Descomposicion de la biomasa en solido-intermediario

 K4=(1.38E10)*exp(161000/(R*(To+(h*t))));%Descomposicion de solido intermediario en carbonizado

 K5=(1.00E5)*exp(108000/(R*(To+(h*t))));%Descomposicion de bio-oil en carbonizado

 K6=(4.28E6)*exp(108000/(R*(To+(h*t))));%Descomposicion de bio-oil en gas

dx = zeros(5,1);    % a column vector

dx(1) = -(K1+K2+K3)*x(1);

dx(2) = K3*x(1)-K4*x(2);

dx(3) = K4*x(2)+K5*x(5);

dx(4)=K1*x(1)+K6*x(5);

dx(5)=K2*x(1)-(K6+K5)*x(5);

 end

>> ode45ecuaciones4