clear all;

%----- input here -----%
%%%%% Set the system parameters here %%%%%
% global s1 s2 s3 ...sn (n : number of parameters)
% here parameters are (a, b, c), so write as below.
global s1 s2 s3 

%----- input here -----%
%%%%% Initial values of parameter %%%%%
% a=1, b=0, c=1
s1=1;
s2=0;
s3=1;

%----- input here -----%
%%%%% Initial values of x %%%%%
initial=[1 1 0 0 0 0 0 0];

%%%%% Solve (n+np) ODE equations %%%%%
% time step : [0,10]
% Initial values : initial
[t,Y] = ode45(@getf,[0 10],initial);

%%%%% Plot the results %%%%%
figure(1);
plot(t,Y(:,3),'-',t,Y(:,5),t,Y(:,7));
figure(2);
plot(t,Y(:,4),'-',t,Y(:,6),t,Y(:,8));