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.15 2003/05/08 18:46:45 gawthrop
## Fixed range bug by using isfield in place of exists
##
## Revision 1.14 2000/05/20 15:35:18 peterg
## Parameterised version
##
## Revision 1.13 1998/09/03 09:12:32 peterg
## Removed debugging lines
##
## Revision 1.12 1998/09/02 12:02:14 peterg
|
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
|
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
|
-
+
-
+
-
+
-
+
|
# Simulation parameters
simpar=$1_simpar;
[nx,ny,nu,nz,nyz] = $1_def;
%Defaults
if (struct_contains(simpar,"wmin")==0)
if (isfield(simpar,"wmin")==0)
simpar.wmin = -1;
end;
if (struct_contains(simpar,"wmax")==0)
if (isfield(simpar,"wmax")==0)
simpar.wmax = 2
end;
if (struct_contains(simpar,"wsteps")==0)
if (isfield(simpar,"wsteps")==0)
simpar.wsteps = 100;
end;
W = logspace(simpar.wmin,simpar.wmax,simpar.wsteps)';
if (struct_contains(simpar,"mttinput")==0)
if (isfield(simpar,"mttinput")==0)
mttinput = 1;
end;
mttu0 = zeros(nu,1);
mttu0(mttinput) = 1;
[n,m]=size(W);
if m>n
|