10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
# Copyright (c) P.J.Gawthrop, 1991, 1994, 1998
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.1 1998/07/19 16:59:56 peterg
## Initial revision
##
###############################################################
# Language
if [ -n "$2" ]; then
|
>
>
>
|
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
# Copyright (c) P.J.Gawthrop, 1991, 1994, 1998
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.2 1998/07/30 10:49:02 peterg
## Put [] round output of zeros.
##
## Revision 1.1 1998/07/19 16:59:56 peterg
## Initial revision
##
###############################################################
# Language
if [ -n "$2" ]; then
|
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
|
BEGIN
write "matrix mtta(", mttnx, ",", mttnx, ");";
write "matrix mttb(", mttnx, ",", mttnu, ");";
write "matrix mttc(", mttny, ",", mttnx, ");";
write "matrix mttd(", mttny, ",", mttnu, ");";
END;
IF (lang = m) THEN
BEGIN
write "[mtta] = zeros(", mttnx, ",", mttnx, ");";
write "[mttb] = zeros(", mttnx, ",", mttnu, ");";
write "[mttc] = zeros(", mttny, ",", mttnx, ");";
write "[mttd] = zeros(", mttny, ",", mttnu, ");";
end;
|
<
<
<
<
<
<
<
<
|
71
72
73
74
75
76
77
78
79
80
81
82
83
84
|
BEGIN
write "matrix mtta(", mttnx, ",", mttnx, ");";
write "matrix mttb(", mttnx, ",", mttnu, ");";
write "matrix mttc(", mttny, ",", mttnx, ");";
write "matrix mttd(", mttny, ",", mttnu, ");";
END;
% find MTTA : the A matrix
FOR j := 1:MTTNx DO
BEGIN
xj := MTTX(j,1);
FOR i := 1:MTTNx DO
|