11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
|
+
+
+
|
# Copyright (c) P.J.Gawthrop, 1996.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.3 1998/05/19 19:27:04 peterg
## Zapped all the parameter stuff
##
## Revision 1.2 1996/09/13 19:40:51 peter
## Fixed problem with default paramaters.
##
## Revision 1.1 1996/09/12 19:26:57 peter
## Initial revision
##
## Revision 1.4 1996/08/16 13:04:46 peter
|
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
|
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
|
+
+
+
+
+
+
+
|
$MATRIX << EOF > ode2odess_m.log 2>mtt_error.txt
%Read in parameters
$1_numpar;
[nx,ny,nu,nz,nyz] = $1_def;
%Read in simulation parameters
$1_simpar;
T = [0:DT:LAST];
t=0; %Just in case it appears in the parameter list.
%Defaults
if exist('T')==0
T=[0:0.1:1]
end;
if exist('x0')==0
|