disser1.pdf (2006 KB) - Ernst-Moritz-Arndt-Universität Greifswald

C.2. Programm 2. Berechnung der Adjungierten.
p=[0.42375, 0.36, 0.566154];
K=[1, 1.2, 1.3];
b=[0, 0.0167, 0.0154; 0.0104, 0, 0.0064; 0.0096, 0.0064,
0];
c=0.125; d=0.06;
s=[0.000375, 0.001333, 0.001231];
x=[0.5792, 0.7071, 0.6920];
w=[0.1253, 0.4094, 3.9333];
t=20; z=0.0001;
for i=20/z+1:-1:w(3)/z+1;
l(i,1)=t;
for j=1:3
l(i,j+1)=(p(j)-c/x(j))*exp(-d*t)/s(j);
end;
t=t-z;
end;
% plot(l(:,1),l(:,2:4));
a1=0;a2=176.7985;a3=219.7581;
x2=0.7071;x3=0.6920;
y0=[l(w(3)/z+1,2:4) w(3) 0.5792];
[t,y]=ode45(@fua5,w(3):-0.0001:w(2),y0);
% plot(t,y(:,1:3));
size(y);
for j=1:size(y,1)-1
l(w(3)/z+1-j,2:4)=y(1+j,1:3);
l(w(3)/z+1-j,1)=y(1+j,4);
end;
% plot(l(:,1),l(:,2:4));
a1=0;a2=176.7985;a3=900;
x2=0.7076;
y0=[l(ceil(w(2)/z+1),2:4) w(2) y(size(y,1),5) 0.692];
[t,y]=ode45(@fua55,w(2):-0.0001:w(1),y0);
% plot(t,y(:,1:3));
size(y);
for j=1:size(y,1)-1
l(ceil(w(2)/z+1)-j,2:4)=y(1+j,1:3);
l(ceil(w(2)/z+1)-j,1)=y(1+j,4);
end;
% plot(l(:,1),l(:,2:4));
a1=0;a2=900;a3=900;
y0=[l(w(1)/z+1,2:4) w(1) y(size(y,1),5) y(size(y,1),6)
0.7071];
[t,y]=ode45(@fua56,w(1):-0.0001:0,y0);
147