12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
|
+
+
+
|
# Copyright (c) P.J.Gawthrop 1991, 1994, 1995, 1996, 1998
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.9 1999/10/26 23:47:58 peterg
## Now does cse.m file
##
## Revision 1.8 1999/07/26 06:07:19 peterg
## Reinstated removal of the odeo.p file
##
## Revision 1.7 1999/03/16 00:38:06 peterg
## Done some more tidying up
##
## Revision 1.6 1999/03/15 07:29:25 peterg
|
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
|
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
|
+
+
+
+
-
+
|
END;;
EOF
fi
if [ "$rep" = "cse" ]; then
touch mtt_odee
touch mtt_odee
fi
sed 's/mtt_matrix/mttdx/' <mtt_ode >> $1_$2.p
sed 's/mtt_matrix/mtte/' <mtt_odee >> $1_$2.p
sed 's/mtt_matrix/mtty/' <mtt_odeo >> $1_$2o.p
if [ "$rep" = "cse" ]; then
cat <<EOF >>$1_$2.p
## Solve E dx = dxe;;
mtt_solve(mttdx,mttE,mttdx,$Nx);;
EOF
fi
rm -f mtt_odeo mtt_ode mtt_odee
# Convert from Pascal syntax to Octave syntax
mtt_p2m<$1_$2o.p > $1_$2o.m
rm -f $1_odeo.p
fi
|