close all
clear all


n_elements = 4
num_state_vars = 2*n_elements;


%Solver Tolerances
rtol = 1e-6;
atol = 1e-5*ones(num_state_vars,1);
options = odeset('RelTol',rtol,'AbsTol',atol);

%Initial Conditions
init_conds = zeros(1,num_state_vars);

%Time range
t = [0 1];

[T,Y] = ode45(@rhs,t,init_conds,options);

figure
plot(T,Y)
