function [out1,out2,out3] = slode(t,y,flag)
% [out1,out2,out3] = slode(t,y,flag)
% Unstable linear ODE x'' - 9x' - 10x = 0
%                     x(0)=1,  x'(0)=-1
% Usage:
%   >> [tspan,y0,options] = ulode([],[],'init')
%   >> [t,y] = ode45('ulode',tspan,y0,options);
%   >> plot(t,y(:,1)), grid

switch flag
 case ''                      % Return dy/dt = f(t,y).
   out1 = f(t,y);
 case 'init'                 % Return default [tspan,y0,options].
   out1 = [0 10];
   out2 = [1; -1];
   out3 = odeset('Stats','on');
 otherwise
   error(['Unknown flag ''' flag '''.']);
 end
% -------------------------------------------------------------
function dydt = f(t,y)
 dydt = [0 1; 10 9] * y;
% -------------------------------------------------------------