25
26
27
28
29
30
31
32
33
34
35
36
37
38
|
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
|
+
+
+
|
######################################
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.14 2002/08/20 16:14:35 gawthrop
## Include Q in documentation
##
## Revision 1.13 2002/08/20 15:43:45 gawthrop
## Works with ident DIY rep
##
## Revision 1.12 2002/06/11 11:25:25 gawthrop
## No longer delay the simulated data.
##
## Revision 1.11 2002/05/20 13:32:36 gawthrop
|
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
|
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
|
-
+
|
##
###############################################################
## Copyright (C) 1999,2000,2001,2002 by Peter J. Gawthrop
## Simulation command
sim_command = sprintf("%s_ssim(x_0,par,simpar,u,i_s)", system_name);
sim_command = sprintf("%s_ssim(x_0,par,simpar,u,i_s);", system_name);
## Extract indices
i_t = free(:,1); # Parameters
i_s = free(:,2)'; # Sensitivities
if nargin<9
extras.criterion = 1e-5;
|