9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
|
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
|
+
+
+
-
-
-
-
+
+
-
+
-
+
|
# Copyright (c) P.J.Gawthrop 1999
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.3 2000/05/19 17:47:56 peterg
## Agument to state
##
## Revision 1.2 2000/05/11 13:43:14 peterg
## No change
##
## Revision 1.1 1999/03/15 23:27:16 peterg
## Initial revision
##
###############################################################
echo Creating $1_odes.dat2
$MATRIX <<EOF >octave_ode2odes.log 2>mtt_error.txt
global MTT_data
## Set up the simulation parameters
par = $1_numpar;
x_0 = $1_state(par);
simpar = $1_simpar;
$1_ode2odes(x_0,par,simpar);
if is_complex(MTT_data)
mtt_data = $1_ode2odes(x_0,par,simpar);
if is_complex(mtt_data)
mtt_error("octave_ode2odes: Simulated data is complex - something is wrong!");
MTT_data = real(MTT_data);
mtt_data = real(mtt_data);
endif;
save -ascii $1_odes.dat2 MTT_data
save -ascii $1_odes.dat2 mtt_data
EOF
mtt_error mtt_error.txt
|