12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
|
for i=1:length(theta)
args="";
for j=1:length(theta)
i_sensitivity=(j==i);
args = sprintf("%s%i %g ",args, i_sensitivity, theta(j));
endfor
args
command = sprintf("./%s_ode2odes.out %s > mtt_data.dat\n", system_name, args)
|
<
|
|
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
for i=1:length(theta)
args="";
for j=1:length(theta)
i_sensitivity=(j==i);
args = sprintf("%s%i %g ",args, i_sensitivity, theta(j));
endfor
command = sprintf("./%s_ode2odes.out %s > mtt_data.dat\n", system_name, args);
system(command);
## Retrieve data
load -force mtt_data.dat
t = mtt_data(:,1);
y = mtt_data(:,2);
y_theta = [y_theta mtt_data(:,3)];
|