21.07.2013 Views

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

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

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

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

C.1 Code <strong>for</strong> the Bifurcation Diagrams in Chapter 8 163<br />

1.601221771499518e-002 ];<br />

<strong>for</strong> r=1:n;<br />

w_6(r+1)=w_6(r)+dw_6;<br />

g=w_6(end);<br />

[M]=newtonRaphson2(@eqdefT1D_full,ini_6);<br />

M=M’;<br />

%<br />

end<br />

s_6(r)=M(end,2,:);<br />

Q_6(r,:)=M(end,:);<br />

V_26(r)=E(1);<br />

V_27(r)=E(2);<br />

V_28(r)=E(3);<br />

V_29(r)=E(4);<br />

V_30(r)=E(5);<br />

w_6(1)=[];<br />

% %---------------------------------------------------------------------<br />

n=80;<br />

w_7=0;<br />

dw_7=0.05e-05;<br />

s_7=zeros(1,n);<br />

Q_6=zeros(n,5);<br />

ini_7 =[4.772255750458808e+005 0 5.600225683160789e-007 0 0];<br />

<strong>for</strong> r=1:n;<br />

w_7(r+1)=w_7(r)+dw_7;<br />

g=w_7(end);<br />

[M]=newtonRaphson2(@eqdefT1D_full,ini_7);

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!