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

164 matlab Code
M=M’;
%
end
s_7(r)=M(end,2,:);
Q_7(r,:)=M(end,:);
V_31(r)=E(1);
V_32(r)=E(2);
V_33(r)=E(3);
V_34(r)=E(4);
V_35(r)=E(5);
w_7(1)=[];
%--------------------------------------------------------------------
k=w(1):0.0005e-05:w_7(end);
%--------------------------------------------------------------------
%--------------------------------------------------------------------
figure
plot(w,s,’-b’,w_2,s_2,’-b’,w_3,s_3,’--r’,w_4,s_4,’--r’,w_5,s_5,’--r’,...
w_6,s_6,’--r’,w_7,s_7,’-g’,w_8,s_8,’--r’)%
axis([0 4e-5 -0.1e6 6e6])
xlabel(’f_1’)
ylabel(’M_a(t)*’)
%ini=[ini_1;ini_2;ini_3;ini_4;ini_5;ini_6;ini_7;ini_8];
%legend([repmat(’ini = ’,size(ini,1),1) num2str(ini)])
Title(’Numerical bifurcation diagram’)
%-----------------------------------------------------
% figure
%
% subplot(2,1,1);plot(w,real(V_1),’.b’,w,real(V_2),’.b’,w,
real(V_3),’.b’,w,real(V_4),’.b’,w,real(V_5),’.b’,
w_2,real(V_6),’.b’,w_2,real(V_7),’.b’,w_2,real(V_8),’.b’,
w_2,real(V_9),’.b’,w_2,real(V_10),’.b’,w_3,real(V_11),
’.b’,w_3,real(V_12),’.b’,w_3,real(V_13),’.b’,w_3,real(V_14),’.b’,
w_3,real(V_15),’.b’,w_4,real(V_16),’.b’,w_4,real(V_17),’.b’,