%%%%%%%%%     MAS.864    %%%%%%%%%%%
%%% carlos rocha - rocha@mit.edu %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


N = 200;
init = zeros(2*(N+1),1);
init(1) = 1;

T = [0:0.1:100];
[t1,Y1] = ode45(@systemcl,T,init);
[t2,Y2] = ode45(@systemsimple,T,[1 0]);

figure(1);
hold on
for i=2:(N+1)
    plot(T,Y1(:,i))
end
plot(T,Y1(:,1),'r')
plot(T,Y2(:,1),'g')