40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
|
p = par; # Reset parameters
p(i) = 1; # Set sensitivity index to 1
[mtt_data] = ${sys}_ode2odes(x0,p,simpar);
if (i==index(1))
y = mtt_data(:,2:2:1+ny2);
endif
y_par = [y_par, mtt_data(:,3:2:2+ny2)];
x = [];
endfor
endif
endfunction
EOF
}
|
|
|
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
|
p = par; # Reset parameters
p(i) = 1; # Set sensitivity index to 1
[mtt_data] = ${sys}_ode2odes(x0,p,simpar);
if (i==index(1))
y = mtt_data(:,2:2:1+ny2);
endif
y_par = [y_par, mtt_data(:,3:2:2+ny2)];
x = mtt_data(:,3+ny2:2+ny2+nx2);
endfor
endif
endfunction
EOF
}
|