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

158 matlab Code
w(r+1)=w(r)+dw;
g=w(end);
[M]=newtonRaphson2(@eqdefT1D_full,ini_1);
M=M’;
end
s(r)=M(end,2,:);
Q(r,:)=M(end,:);
V_1(r)=E(1);
V_2(r)=E(2);
V_3(r)=E(3);
V_4(r)=E(4);
V_5(r)=E(5);
w(1)=[];
%
% %--------------------------------------------------------------------
%
n=855;
w_2=8.6e-6;
dw_2=0.002e-05;
s_2=zeros(1,n);
Q_2=zeros(n,5);
V=zeros(n,5);
ini_2 =[1.177753840657023e+006 5.017144005834548e+006...
1.967867544648059e+005 1.588274974061828e+003...
1.593754564310941e+000]
for r=1:n;
w_2(r+1)=w_2(r)+dw_2;
g=w_2(end);
[M]=newtonRaphson2(@eqdefT1D_full,ini_2);
M=M’;
Q_2(r,:)=M(end,:);