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

C.2 Code for the Newton-Raphson Method 169
end
E=eig(jac);
%-----------------------------------------
%The code to "call" when solving the system
of differential equations by ode45.
function Mdot = eqdefT1D_full(M)
global g;
b = 0.09; %defines the parameter values
c = 0.1;
d = 0.5;
e = 1e-8;
f2=1e-5;
f1=1e-5
a = 5e4;
k = 0.4;
alpha=5e-9;
delta=25;
kc=1;
Amax = 2e7;
kb=(delta/alpha)*kc;
%DuCa model without apoptotic wave
Mdot = [a + (b + k) * M(2) - c * M(1) -
f1 * M(1) * M(3) - e * M(1) * ( M(1) + M(2) );
;f1 * M(1) * M(3) - k * M(2) - e * M(2) * ( M(1) + M(2) );
;(( Amax * M(5) )./( kc + M(5) )) - f1 * M(1) * M(3)
- f2 * M(2) * M(3) - d * M(3);
;d * M(3) - f1 * M(1) * M(4) - f2 * M(2)* M(4);
;alpha * M(4) * M(2) - delta * M(5)];
%------------------------------
%A sample code used for solving the system of differential
%equations by ode15s
clear all
close all
% initial values
M = [4e3 10e6 2e3 2e3 2e3];
% time span
tspan = [0 200];