11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
# Copyright (c) P.J.Gawthrop, 1996.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.4 2000/10/17 11:07:59 peterg
## *** empty log message ***
##
## 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
|
>
>
>
|
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.1 2000/10/17 12:55:29 peterg
## Initial revision
##
## Revision 1.4 2000/10/17 11:07:59 peterg
## *** empty log message ***
##
## 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
|
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
|
for t=T'
u = $1_input(x,zeros($Ny,1),t,par);
x = fsolve('f',x);
y = $1_odeo(x,u,t,par);
mtt_write(t,x,y,$Nx,$Ny); # Write it out
end;
save -ascii $1_odess.dat2 MTT_data
EOF
# Now invoke the standard error handling.
mtt_error mtt_error.txt
|
|
|
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
|
for t=T'
u = $1_input(x,zeros($Ny,1),t,par);
x = fsolve('f',x);
y = $1_odeo(x,u,t,par);
mtt_write(t,x,y,$Nx,$Ny); # Write it out
end;
save -text $1_odess.dat2 MTT_data
EOF
# Now invoke the standard error handling.
mtt_error mtt_error.txt
|