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.3 1997/05/01 11:15:33 peterg
# Back under RCS
#
# Revision 1.2 1997/03/20 14:36:56 peterg
# Includes the sympar.h file
#
## Revision 1.1 1997/01/21 22:54:54 peterg
## Initial revision
##
###############################################################
|
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
|
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
|
-
-
-
-
+
+
+
+
-
-
+
+
|
write "#include ""$1_odes.h"" "$
%External (global) variable list
write "#include ""$1_sympar.h"" "$
write "/* Declare standard arrays */"$
write "double y[", MTTNy+1, "]; /* $1_ode output */"$
write "double dx[", MTTNx+1, "]; /* $1_ode state derivative */"$
write "double x[", MTTNx+1, "]; /* $1_ode state */"$
write "double u[", MTTNu+1, "]; /* $1_ode input */"$
write "float y[", MTTNy+1, "]; /* $1_ode output */"$
write "float dx[", MTTNx+1, "]; /* $1_ode state derivative */"$
write "float x[", MTTNx+1, "]; /* $1_ode state */"$
write "float u[", MTTNu+1, "]; /* $1_ode input */"$
write "/* Files */ "$
write " FILE *fopen(), *fps, *fpso;"$
write "main()"$
write " "$
write "{"$
write "/* Counters etc*/ "$
write " double time;"$
write " double dt;"$
write " float time;"$
write " float dt;"$
write " int i;"$
write " int k;"$
write "/*functions */ "$
write " extern $1_numpar();"$
%Open the output files
|