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.14 1998/05/13 08:57:27 peterg
## Now uses simpar.h in place of odes.h
##
## Revision 1.13 1998/02/24 13:34:45 peterg
## Back under RCS
##
# Revision 1.12 1997/05/15 08:39:56 peterg
# Don't initialise states - now done in numpar file.
#
# Revision 1.11 1997/05/12 16:00:54 peterg
|
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
|
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
|
-
+
|
%Set up initial state
write "/* Set up initial state */"$
write " $1_state();"$
%Initialise main (Euler) integration loop
write "/* Initialise main (Euler) integration loop */"$
write " time = 0.0;"$
write " time = DT;"$
write " dt = DT/STEPFACTOR;"$
%% 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++)"$
|