︙ | | |
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
+
+
+
|
# Copyright (c) P.J.Gawthrop 1991, 1994, 1995, 1996
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.13 1998/05/21 16:19:54 peterg
## Modified to include explicit algebraic loop solution
##
## Revision 1.12 1998/05/21 12:55:48 peterg
## Put in algebraic equation stuff
##
## Revision 1.11 1998/05/21 08:05:23 peterg
## Back under RCS
##
## Revision 1.10 1998/04/14 07:25:02 peterg
|
︙ | | |
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
|
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
|
-
+
|
cat $1_ode.m4 >> $1_ode.m
cat <<EOF >> $1_ode.m
% Solve the algebraic equations (if any)
if nyz>0
global xx tt;
xx = x; tt=t;
mttui = fsolve('$1_odea',mttui);
MTTui = fsolve('$1_odea',mttui);
else
mttui = [];
end
EOF
cat $1_ode.m3 >> $1_ode.m
|
︙ | | |
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
|
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
|
-
+
|
# Create the odea.m function
cat <<EOF > $1_odea.m
function mttyz = $1_odea(mttui);
% mttyz = $1_odea(mttui);
%Algebraic equations in Octave form for system $1;
%File $1_odea.m;
%Generated by MTT on `date`;
%Generated by mtt on `date`;
EOF
# Create the globals
sympar2global_txt2m $1 >> $1_odea.m
cat <<EOF >> $1_odea.m
|
︙ | | |
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
|
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
|
-
+
|
# Create the odeo.m function
cat <<EOF > $1_odeo.m
function mtty = $1_odeo(x,t);
% mtty = $1_odeo(x,t);
%Algebraic equations in Octave form for system $1;
%File $1_odeo.m;
%Generated by MTT on `date`;
%Generated by MTT on `date`;
EOF
# Create the globals
sympar2global_txt2m $1 >> $1_odeo.m
#Common bit
|
︙ | | |