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