nr. 477 - 2011 - Institut for Natur, Systemer og Modeller (NSM)

C.1 Code for the Bifurcation Diagrams in Chapter 8 165
w_4,real(V_18),’.b’,w_4,real(V_19),’.b’,w_4,real(V_20),’.b’,
w_8,real(V_36),’.b’,w_8,real(V_37),’.b’,
w_8,real(V_38),’.b’,w_8,real(V_39),’.b’,w_8,real(V_40),’.b’,k,0,’black’)
% axis([0 3.5e-5 -4 4 ])
% Title(’Plot of eigenvalues for the Jacobian evaluated in the upper inflamed region’)
% xlabel(’f_1’)
% ylabel(’real(eig)’)
%
% subplot(2,1,2);plot(w,imag(V_1),’.b’,w,imag(V_2),’.b’,
w,imag(V_3),’.b’,w,imag(V_4),’.b’,w,imag(V_5),’.b’,
w_2,imag(V_6),’.b’,w_2,imag(V_7),’.b’,w_2,imag(V_8),’.b’,
w_2,imag(V_9),’.b’,w_2,imag(V_10),’.b’,w_3,imag(V_11),’.b’,
w_3,imag(V_12),’.b’,w_3,imag(V_13),’.b’,w_3,imag(V_14),’.b’,
w_3,imag(V_15),’.b’,w_4,imag(V_16),’.b’,w_4,imag(V_17),’.b’,
w_4,imag(V_18),’.b’,w_4,imag(V_19),’.b’,w_4,imag(V_20),’.b’,
w_8,imag(V_36),’.b’,w_8,imag(V_37),’.b’,w_8,imag(V_38),’.b’,
w_8,imag(V_39),’.b’,w_8,imag(V_40),’.b’,k,0,’black’)
% xlabel(’f_1’)
% ylabel(’imag(eig)’)
%
% --------------figure
subplot(2,1,1);plot(w_5,real(V_21),’.b’,w_5,real(V_22),’.b’,
w_5,real(V_23),’.b’,w_5,real(V_24),’.b’,
w_5,real(V_25),’.b’,w_6,real(V_26),’.b’,w_6,real(V_27),’.b’,
w_6,real(V_28),’.b’,w_6,real(V_29),’.b’,
w_6,real(V_30),’.b’,k,0,’black’)
axis([0 3.5e-5 -2 2 ])
Title(’Plot of eigenvalues for the Jacobian evaluated in the lower inflamed region ’)
xlabel(’f_1’)
ylabel(’real(eig)’)
subplot(2,1,2);plot(w_5,imag(V_21),’.b’,w_5,imag(V_22),’.b’,
w_5,imag(V_23),’.b’,w_5,imag(V_24),’.b’,
w_5,imag(V_25),’.b’,w_6,imag(V_26),’.b’,w_6,imag(V_27),’.b’,
w_6,imag(V_28),’.b’,w_6,imag(V_29),’.b’,
w_6,imag(V_30),’.b’,k,0,’black’)
xlabel(’f_1’)
ylabel(’imag(eig)’)
%--------------figure
plot(w_7,V_31,’.b’,w_7,V_32,’.b’,w_7,V_33,’.b’,