nr. 477 - 2011 - Institut for Natur, Systemer og Modeller (NSM)
nr. 477 - 2011 - Institut for Natur, Systemer og Modeller (NSM)
nr. 477 - 2011 - Institut for Natur, Systemer og Modeller (NSM)
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
C.3 The Function-File <strong>for</strong> Model C 171<br />
Mdot(5,1) = h*M(4)*M(2) - m*M(5);<br />
Mdot(6,1) = -(a0*(1+exp(-p2*t2))/(1+exp(p2*(t-t2))) + 0.001)*M(6)...<br />
- (4*10^7 *exp(-((t-9)/3)^2)/(1e9))*M(6) - (A*M(5)/(k_c + M(5)))*M(6) ...<br />
+ r0*(1+exp(-p1*t1))/(1+exp(p1*(t-t1)))*M(6);<br />
Simulation Code <strong>for</strong> Model C<br />
The following code was used to produce the simulations of Model C.<br />
clear all<br />
close all<br />
figure<br />
hold on<br />
M = [4.77e5,0,0,0,0,1e9];...<br />
tspan = [0;5000];<br />
options = odeset(’Reltol’,1e-20,’AbsTol’,1e-20);<br />
[t,M] = ode15s(@mareemodelfuld6C,tspan,M);<br />
Mtilde=M(:,6);<br />
Mtilde(find(Mtilde
Change language