1NAME
2Reentry
3
4TITLE
5Stoer-Bulirsch, Kap. 7
6
7CONSTANTS
8R = 209
9beta = 4.26
10rho0 := 2.704e-3
11g := 0.32172e-3
12
13DYNAMIC_VARIABLES
14v
15gama
16xi
17
18CONTROL_VARIABLES
19u
20
21ABBREVIATIONS
22rho(~h) = rho0 * exp(- beta * h);
23c_W (u) :=  1.174 - 0.9 * cos(u);
24c_A (u) := 0.6 * sin(u);
25S := 53200 * m;
26h := xi * R;
27
28TIME_VARIABLE
29tt
30
31DIFFERENTIAL_EQUATIONS
32DOT(v) := -S*rho(h)*v^2*c_W(u)/(2*m)
33          - g * sin(gama) / (1 + xi)^2;
34DOT(gama) := S*rho(h)*v*c_A(u)/(2*m)
35              + v*cos(gama) / (R * (1 + xi))
36              - g*cos(gama)/(v * (1 + xi)^2);
37xi       v * sin (gama) / R;
38
39FUNCTIONAL
4010*v^3*sqrt(rho(h));
41
42END_POINT_IS_FREE
43
44BOUNDARIES
450;
46TIME_T;
47
48BOUNDARY_CONDITIONS
49v(0) - 0.36;
50gama(0) + 8.1*pi/180;
51xi(0)-4/R;
52v(TIME_T) - 0.27;
53gama(TIME_T);
54xi(TIME_T) -2.5/R;
55
56RESTRICTIONS
57v > 0;
58
59