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

ANHANG C
MATLAB-Programme
C.1. Programm 1. Berechnung stationärer
Zustände.
clear all; format long;
p=[0.42375, 0.36, 0.566154];
e=[0.4, 0.6, 0.6];
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];
for i=1:3
l(i)=0;
for j=1:3
l(i)=l(i)+s(i)/s(j)*c*b(j,i);
end;
end;
for k=1:3
ps1(k)=-2*p(k)*e(k)/K(k);
ps2(k)=-d*p(k)+p(k)*e(k)+c*e(k)/K(k)+l(k);
ps3(k)=d*c;
ps4(1)=-p(1)*b(1,2)-(s(1)/s(2))*p(2)*b(2,1);
ps4(2)=-p(2)*b(2,1)-(s(2)/s(1))*p(1)*b(1,2);
ps4(3)=-p(3)*b(3,1)-(s(3)/s(1))*p(1)*b(1,3);
ps5(1)=-p(1)*b(1,3)-(s(1)/s(3))*p(3)*b(3,1);
ps5(2)=-p(2)*b(2,3)-(s(2)/s(3))*p(3)*b(3,2);
ps5(3)=-p(3)*b(3,2)-(s(3)/s(2))*p(2)*b(2,3);
145