clear;
N=1000000;I_0=1;eta=3;gamma=0.05;
sim('prj2_model2')
figure; 
plot(J.time,J.signals.values,'g',I.time,I.signals.values,'k--',Q.time,Q.signals.values,'b-^')  %%% use this when using format ”structure with time”.
% plot(J,'g')
% hold on; plot(I,'k--')
% plot(Q,'b-^')% use these three plots when using format ”array”
legend('J(t)','I(t)','Q(t)')
xlabel('t')
ylabel('number');
hold off