Matlab - Problema Ode45

close all

clear

clc

k1=5;

k0=0.1;

k2=1;

S0=5;

E0=0.1;

f= @(t,n)[ k0*n(1)-k1*n(2)*(E0-k0); k1*n(2)*(E0-n(1))-(k0+k2)*n(1)  ];

 t= 0:20;

  [nsol,t]=ode45(f,t,[0 S0]);

    plot(t,nsol)

    xlabel('t');

    ylabel('C');

>> ode45caballero15