16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
+
+
+
|
# Copyright (c) P.J.Gawthrop 1997.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
# Revision 1.11 1997/05/12 16:00:54 peterg
# Removed itime again.,
#
# Revision 1.10 1997/05/10 10:05:15 peterg
# Put _input in inner loop in front of call to _ode
#
## Revision 1.9 1997/05/10 09:54:34 peterg
## Moved _input call to after the inner integration loop.
##
## Revision 1.8 1997/05/10 08:12:23 peterg
|
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
|
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
|
-
-
+
+
|
write "/* Set up user-defined constants */"$
write " $1_numpar();"$
%Initialise main (Euler) integration loop
write "/* Initialise main (Euler) integration loop */"$
write " time = 0.0;"$
write " dt = DT/STEPFACTOR;"$
write " for (i=1; i<=MTTNX; i++)"$
write " x[i] = 0.0;"$
%% write " for (i=1; i<=MTTNX; i++)"$
%% write " x[i] = 0.0;"$
%Set up system inputs
write "/* Set up system inputs */"$
write " for (i=1; i<=", MTTNu, "; i++)"$
write " u[i] = 1.0;"$
write " $1_input(0.0);"$
|