ci = [0.000, 0.3177, 0.0530, 0.5959]; %% estos son los valores estacionarios [t,x] = ode45 ('HH',[0,40],ci); v=x(:,1); % variable v n=x(:,2); % variable w m=x(:,3); % variable w h=x(:,4); % variable w figure plot(t,n,t,m,t,h) % Grafico de la evolucion temporal legend ('n','m','h'); figure plot(t,v) axis([0 40 -20 120])