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

C.4. Programm 4. Das Integral-Maximumprinzip.
y(2)*0.42375*0.000375 ...
*y(1))/(2*0.00125))-0.4*y(2)+2*0.4*y(2)*y(1)+ ...
((y(2)*0.42375*0.000375*y(1)-y(2)^2*0.000375^2 ...
*y(1)*exp(0.06*y(3)))/(2*0.00125))
1];
C.4. Programm 4. Das Integral-Maximumprinzip.
clear all;
options = [];
sol = bvpinit(linspace(0,1,50),[1 1 1]);
sol = bvp4c(@odesys4,@bc4,sol,options);
x = [sol.x];
y1 = [sol.y(1,:)]; y2 = [sol.y(2,:)]; y3 = [sol.y(3,:)];
y=[y1;y2;y3];
x=linspace(0,1,50);
f=((1-0.000375*38.519/0.4)*0.25* ...
exp((0.4-0.000375*38.519).*x)./ ...
((1-0.000375*38.519/0.4)+0.25.* ...
(exp((0.4-0.000375*38.519).*x)-1)));
z=0.42375.*f.*exp(-0.06.*x)-2*0.00125*38.519.* ...
exp(-0.06.*x)-y2.*0.000375.*f;
l = trapz(x,z),
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function dydx = odesys4(x,y,n)
% ODE function
dydx = [ 0.4*y(1)*(1-y(1))-0.000375*(38.519*y(1))
(-0.42375*38.519*exp(-0.06*y(3))-0.4*y(2)+ ...
2*0.4*y(2)*y(1)+ 0.000375*38.519*y(2))
1];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function res = bc4(ya,yb)
% Boundary conditions
res = [ ya(1)-0.25
yb(2)-206.46
ya(3)];
C.5. Programm 5. Logarithmisches Modell.
clear all
options = [];
149