We now discuss the numerical approximation of solutions of systems of differential equations. Our goal is to apply the methods of Sections2.4 through 2.6 to the initial value problem
x′=f(t,x),x(t0)=x0
x'=f(t,x),x(t0)=x0
(1)
for a system of m first-order differential equations. In (1) the independent variable is the scalar t, and
x=(x1,x2,…,xm)andf=(f1,f2,…,fm)
x=(x1,x2,…,xm)andf=(f1,f2,…,fm)
are vector-valued functions. If the component functions of ff and their first-order partial derivatives are all continuous in a neighborhood of the point (t0,x0)(t0,x0), then Theorems 3 and 4 of Appendix A guarantee the existence and uniqueness of a solution x=x(t)x=x(t) of (1) on some subinterval [of the t-axis] containing t0t0. With this assurance we can proceed to discuss the numerical approximation of this solution.
Beginning with step size h, we want to approximate the values of x(t)x(t) at the points t1,t2,t3,…,t1,t2,t3,…, where tn+1=tn+htn+1=tn+h for n≧0n≧0. Suppose that we have already computed the approximations
x1,x2,x3,…,xn
x1,x2,x3,…,xn
to the actual values
x(t1),x(t2),x(t3),…,x(tn)
x(t1),x(t2),x(t3),…,x(tn)
of the exact solution of the system in (1). We can then make the step from xnxn to the next approximation xn+1≈x(tn+1)xn+1≈x(tn+1) by any one of the methods of Sections2.4 through 2.6. Essentially all that is required is to write the iterative formula of the selected method in the vector notation of the present discussion.
Euler Methods for Systems
For example, the iterative formula of Euler’s method for systems is
xn+1=xn+hf(t,xn).
xn+1=xn+hf(t,xn).
(2)
To examine the case m=2m=2 of a pair of first-order differential equations, let us write
and the scalar components of the vector formula in (2) are
xn+1=xn+hf(tn,xn,yn),yn+1=yn+hg(tn,xm,yn).
xn+1yn+1==xnyn++hf(tn,xn,yn),hg(tn,xm,yn).
(4)
Note that each iterative formula in (4) has the form of a single Euler iteration, but with ynyn inserted like a parameter in the first formula (for xn+1xn+1) and with xnxn inserted like a parameter in the second formula (for yn+1yn+1). The generalization to the system in (3) of each of the other methods in Sections2.4 through 2.6 follows a similar pattern.
The improved Euler method for systems consists at each step of calculating first the predictor
un+1=xn+hf(tn,xn)
un+1=xn+hf(tn,xn)
(5)
and then the corrector
xn+1=xn+h2[f(tn,xn)+f(tn+1,un+1)].
xn+1=xn+h2[f(tn,xn)+f(tn+1,un+1)].
(6)
For the case of the two-dimensional initial value problem in (3), the scalar components of the formulas in (5) and (6) are
To describe in scalar notation the Runge–Kutta method for the two-dimensional initial value problem
x′=f(t,x,y),x(t0)=x0,y′=g(t,x,y),y(t0)=y0,
x'y'==f(t,x,y),g(t,x,y),x(t0)y(t0)==x0,y0,
(3)
let us write
x=[xy],f=[fg],andki=[FiGi].
x=[xy],f=[fg],andki=[FiGi].
Then the Runge–Kutta iterative formulas for the step from (xn,yn)(xn,yn) to the next approximation (xn+1,yn+1)≈(x(tn+1),y(tn+1))(xn+1,yn+1)≈(x(tn+1),y(tn+1)) are
G1,G2,G3,G1,G2,G3, and G4G4 are the similarly defined values of the function g.
Perhaps the most common application of the two-dimensional Runge–Kutta method is to the numerical solution of second-order initial value problems of the form
x″=g(t,x,x′),x(t0)=x0,x′(t0)=y0.
x′′=g(t,x,x'),x(t0)=x0,x'(t0)=y0.
(15)
If we introduce the auxiliary variable y=x′,y=x', then the problem in (15) translates into the two-dimensional first-order problem
x′=y,x(t0)=x0,y′=g(t,x,y),y(t0)=y0.
x'y'==y,g(t,x,y),x(t0)y(t0)==x0,y0.
(16)
This is a problem of the form in (3) with f(t,x,y)=y.f(t,x,y)=y.
If the functions f and g are not too complicated, then it is feasible to carry out manually a reasonable number of steps of the two-dimensional Runge–Kutta method described here. But the first operating electronic computers were constructed (during World War II) specifcally to implement methods similar to the Runge–Kutta method for the numerical computation of trajectories of artillery projectiles. The application material for this section lists TI-Nspire CX CAS and Python versions of Program RK2DIM that can be used with two-dimensional systems.
Example2
The exact solution of the initial value problem
x″=−x;x(0)=0,x′(0)=1
x''=−x;x(0)=0,x'(0)=1
(17)
is x(t)=sint.x(t)=sint. The substitution y=x′y=x' translates (17) into the two-dimensional problem
x′=y,x(0)=0;y′=−x,y(0)=1,
x'y'==y,−x,x(0)y(0)==0;1,
(18)
which has the form in (3) with f(t,x,y)=yf(t,x,y)=y and g(t,x,y)=−x.g(t,x,y)=−x. The table in Fig.7.7.1 shows the results produced for 0≦t≦50≦t≦5 (radians) using Program RK2DIM with step size h=0.05.h=0.05. The values shown for x=sint and y=costx=sint and y=cost are all accurate to five decimal places.
FIGURE7.7.1.
Runge–Kutta values (with h=0.05h=0.05) for the problem in Eq.(18).
t
x=sintx=sint
y=costy=cost
0.5
+0.47943+0.47943
+0.87758+0.87758
1.0
+0.84147+0.84147
+0.54030+0.54030
1.5
+0.99749+0.99749
+0.07074+0.07074
2.0
+0.90930+0.90930
−0.41615−0.41615
2.5
+0.59847+0.59847
−0.80114−0.80114
3.0
+0.14112+0.14112
−0.98999−0.98999
3.5
−0.35078−0.35078
−0.93646−0.93646
4.0
−0.75680−0.75680
−0.65364−0.65364
4.5
−0.97753−0.97753
−0.21080−0.21080
5.0
−0.95892−0.95892
+0.28366+0.28366
Example3
Lunar lander In Example4 of Section2.3 we considered a lunar lander that initially is falling freely toward the surface of the moon. Its retrorockets, when fired, provide a deceleration of T=4m/s2.T=4m/s2. We found that a soft touchdown at the lunar surface is achieved by igniting these retrorockets when the lander is at a height of 41,870 meters (just over 26 miles) above the surface and is then descending at the rate of 450 m/s.
Now we want to compute the descent time of the lunar lander. Let the distance x(t) of the lander from the center of the moon be measured in meters and measure time t in seconds. According to the analysis in Section2.3 (where we used r(t) in stead of x(t)), x(t) satisfes the initial value problem
where G≈6.6726×10−11 N⋅(m/kg)2G≈6.6726×10−11 N⋅(m/kg)2 is the universal gravitational constant and M=7.35×1022 kgM=7.35×1022 kg and R=1.74×106 mR=1.74×106 m are the mass and radius of the moon. We seek the value of t when x(t)=R=1,740,000.x(t)=R=1,740,000..
The problem in (19) is equivalent to the first-order system
The table in Fig.7.7.2 shows the result of a Runge–Kutta approximation with step size h=1h=1 (the indicated data agreeing with those obtained with step size h=2h=2). Evidently, touchdown on the lunar surface (x=1,740,000)(x=1,740,000) occurs at some time between t=180t=180 and t=190t=190 seconds. The table in Fig.7.7.3 shows a second Runge–Kutta approximation with t(0)=180,x(0)=1,740,059,y(0)=−16.83,t(0)=180,x(0)=1,740,059,y(0)=−16.83, and h=0.1.h=0.1. Now it is apparent that the lander’s time of descent to the lunar surface is very close to 187 seconds; that is, 3 min 7 s. (The final velocity terms in these two tables are positive because the lander would begin to ascend if its retrorockets were not turned off at touchdown.)
FIGURE7.7.2.
The lander’s descent to the lunar surface.
t(s)
x(m)
v(m/s)
0
1,781,870
−450.00−450.00
20
1,773,360
−401.04−401.04
40
1,765,826
−352.37−352.37
60
1,759,264
−303.95−303.95
80
1,753,667
−255.74−255.74
100
1,749,033
−207.73−207.73
120
1,745,357
−159.86−159.86
140
1,742,637
−112.11−112.11
160
1,740,872
−64.45−64.45
180
1,740,059
−16.83−16.83
200
1,740,199
30.77
FIGURE7.7.3.
Focusing on the lunar lander’s soft touchdown.
t(s)
x(m)
v(m/s)
180
1,740,059
−16.83−16.83
181
1,740,044
−14.45−14.45
182
1,740,030
−12.07−12.07
183
1,740,019
−9.69−9.69
184
1,740,011
−7.31−7.31
185
1,740,005
−4.93−4.93
186
1,740,001
−2.55−2.55
187
1,740,000
−0.17−0.17
188
1,740,001
2.21
189
1,740,004
4.59
190
1,740,010
6.97
Higher-Order Systems
As we saw in Section7.1, any system of higher-order differential equations can be replaced with an equivalent system of first-order differential equations. For example, consider the system
of four first-order equations in the unknown functions x1(t)=x(t),x2(t)=y(t),x3(t), and x4(t).x1(t)=x(t),x2(t)=y(t),x3(t), and x4(t). It would be a routine (if slightly tedious) matter to write a four-dimensional version of program RK2DIM for the purpose of solving such a system. But in a programming language that accommodates vectors, an n-dimensional Runge–Kutta program is scarcely more complicated than a one-dimensional program. For instance, the application material for this section lists the n-dimensional Matlab program rkn that closely resembles the one-dimensional program rk of Fig.2.6.11.
Example4
Batted baseball Suppose that a batted ball starts at x0=0,y0=0x0=0,y0=0 with initial velocity v0=160 ft/sv0=160 ft/s and with initial angle of inclination θ=30∘.θ=30∘. If air resistance is ignored, we find by the elementary methods of Section1.2 that the baseball travels a [horizontal] distance of 400√3 ft4003–√ ft (approximately 693 ft) in 5 s before striking the ground. Now suppose that in addition to a downward gravitational acceleration (g=32 ft/s2),(g=32 ft/s2), the baseball experiences an acceleration due to air resistance of (0.0025)v2(0.0025)v2 feet per second per second, directed opposite to its instantaneous direction of motion. Determine how far the baseball will travel horizontally under these conditions.
Solution
According to Problem30 of Section7.1, the equations of motion of the baseball are
d2xdt2=−cvdxdt,d2ydt2=−cvdydt−g
d2xdt2=−cvdxdt,d2ydt2=−cvdydt−g
(23)
where v=√(x′)2+(y′)2v=(x')2+(y')2−−−−−−−−−−√ is the speed of the ball, and where c=0.0025c=0.0025 and g=32g=32 in fps units. We convert to a first-order system as in (22) and thereby obtain the system
Note that x3(t) and x4(t)x3(t) and x4(t) are simply the x- and y-components of the baseball’s velocity vector, so v=√x23+x24.v=x23+x24−−−−−−√. We proceed to apply the Runge–Kutta method to investigate the motion of the batted baseball described by the initial value problem in (24) and (25), first taking c=0c=0 to ignore air resistance and then using c=0.0025c=0.0025 to take air resistance into account.
Without Air Resistance:Figure7.7.4 shows the numerical results obtained when a Runge–Kutta program such as rkn is applied with step size h=0.1h=0.1 and with c=0c=0 (no air resistance). For convenience in interpreting the results, the printed output at each selected step consists of the horizontal and vertical coordinates x and y of the baseball, its velocity v, and the angle of inclination αα of its velocity vector (in degrees measured from the horizontal). These results agree with the exact solution when c=0.c=0. The ball travels a horizontal distance of 400√3≈692.82 ft4003–√≈692.82 ft in exactly 5 s, having reached a maximum height of 100 ft after 2.5 s. Note also that the ball strikes the ground at the same angle and with the same speed as its initial angle and speed.
FIGURE7.7.4.
The batted baseball with no air resistance (c=0)(c=0).
t
x
y
v
αα
0.0
0.00
0.00
160.00
+30+30
0.5
69.28
36.00
152.63
+25+25
1.0
138.56
64.00
146.64
+19+19
1.5
207.85
84.00
142.21
+13+13
2.0
277.13
96.00
139.48
+7+7
2.5
346.41
100.00
138.56
+0+0
3.0
415.69
96.00
139.48
−7−7
3.5
484.97
84.00
142.21
−13−13
4.0
554.26
64.00
146.64
−19−19
4.5
623.54
36.00
152.63
−25−25
5.0
692.82
0.00
160.00
−30−30
FIGURE7.7.5.
The batted baseball with air resistance (c=0.0025)(c=0.0025).
t
x
y
v
αα
0.0
0.00
0.00
160.00
+30+30
0.5
63.25
32.74
127.18
+24+24
1.0
117.11
53.20
104.86
+17+17
1.5
164.32
63.60
89.72
+8+8
2.0
206.48
65.30
80.17
−3−3
2.5
244.61
59.22
75.22
−15−15
3.0
279.29
46.05
73.99
−27−27
3.5
310.91
26.41
75.47
−37−37
4.0
339.67
0.91
78.66
−46−46
FIGURE7.7.6.
The batted ball’s apex and its impact with the ground.
t
x
y
v
αα
1.5
164.32
63.60
89.72
+8+8
1.6
173.11
64.60
87.40
+5+5
1.7
181.72
65.26
85.29
+3+3
1.8
190.15
65.60
83.39
+1+1
← Apex
1.9
198.40
65.61
81.68
−1−1
2.0
206.48
65.30
80.17
−3−3
⋮⋮
⋮⋮
⋮⋮
⋮⋮
⋮⋮
3.8
328.50
11.77
77.24
−42−42
3.9
334.14
6.45
77.93
−44−44
4.0
339.67
0.91
78.66
−46−46
← Impact
4.1
345.10
−4.84−4.84
79.43
−47−47
4.2
350.41
−10.79−10.79
80.22
−49−49
With Air Resistance:Figure7.7.5 shows the results obtained with the fairly realistic value of c=0.0025c=0.0025 for the air resistance for a batted baseball. To within a hundredth of a foot in either direction, the same results are obtained with step sizes h=0.05h=0.05 and h=0.025.h=0.025. We now see that with air resistance the ball travels a distance well under 400 ft in just over 4 s. The more refined data in Fig.7.7.6 show that the ball travels horizontally only about 340 ft and that its maximum height is only about 66 ft. As illustrated in Fig.7.7.7, air resistance has converted a massive home run into a routine fly ball (if hit straightaway to center field). Note also that when the ball strikes the ground, it has slightly under half its initial speed (only about 79 ft/s) and is falling at a steeper angle (about 46∘46∘). Every baseball fan has observed empirically these aspects of the trajectory of a fly ball.
Variable Step Size Methods
The Runge–Kutta method for a large system requires an appreciable amount of computational labor, even when a computer is employed. Therefore, just as the step size h should not be so large that the resulting error in the solution is unacceptable, h ought not to be so small that too many steps are needed, hence requiring an unacceptable amount of computation. Thus the practical numerical solution of differential equations involves a tradeoff between accuracy and efficiency.
To facilitate this tradeoff, modern variable step size methods vary the step size h as the solution process proceeds. Large steps are taken in regions where the dependent variables are changing slowly; smaller steps are taken when these variables are changing rapidly, in order to prevent large errors.
An adaptable or variable step size Runge–Kutta method employs both a pre-assigned minimum error tolerance MinTol and a maximum error tolerance MaxTol to attempt to ensure that the error made in the typical step from xnto xn+1xnto xn+1 is neither too large (and hence inaccurate) nor too small (and hence in efficient). A fairly simple scheme for doing this may be outlined as follows:
Having reached xnxn with a Runge–Kutta step of length tn−tn−1=h, let x(1)tn−tn−1=h, let x(1) denote the result of a further Runge–Kutta step of length h and let x(2)x(2) denote the result of two successive Runge–Kutta steps each of length h/2.
On the grounds that x(2)x(2) should be a more accurate approximation to x(tn+h)x(tn+h) than is x(1),x(1), take
Err=|x(1)−x(2)|
Err=|x(1)−x(2)|
as an estimate of the error in x(1).x(1).
If MinTol≦Err≦MaxTol, then let xn+1=x(1),tn+1=tn+h,MinTol≦Err≦MaxTol, then let xn+1=x(1),tn+1=tn+h, and proceed to the next step.
If Err<MaxTol,Err<MaxTol, then the error is too small! Hence let xn+1=x(1),tn+1=tn+h,xn+1=x(1),tn+1=tn+h, but double the step size to 2h before making the next step.
If Err>MaxTol,Err>MaxTol, then the error is too large. Hence reject x(1)x(1) and start afresh at xnxn with the halved step size h/2.
The detailed implementation of such a scheme can be complicated. For a much more complete but readable discussion of adaptive Runge–Kutta methods, see Section 17.2 of William H. Press et al., Numerical Recipes: The Art of Scientific Computing (Cambridge: Cambridge: University Press, 2007).
Several widely available scientific computing packages (such as Maple, Mathematica, and Matlab) include sophisticated variable step size programs that will accommodate an essentially arbitrary number of simultaneous differential equations. Such a general-purpose program might be used, for example, to model numerically the major components of the solar system:the sun and the nine (known) major planets. If mimi denotes the mass and ri=(xi,yi,zi)ri=(xi,yi,zi) denotes the position vector of the ith one of these 10 bodies, then—by Newton’s laws—the equation of motion of mimi is
mir″=∑j≠iGmimj(rij)3(rj−ri),
mir''=∑j≠iGmimj(rij)3(rj−ri),
(26)
where rij=|rj−ri|rij=|rj−ri| denotes the distance between mi and mj.mi and mj. For each i=1,2,…,i=1,2,…, 10, the summation in Eq.(26) is over all values of j≠ij≠i from 1 to 10. The 10 vector equations in (26) constitute a system of 30 second-order scalar equations, and the equivalent first-order system consists of 60 differential equations in the coordinates and velocity components of the 10 major bodies in the solar system. Mathematical models that involve this many (or more) differential equations—and that require sophisticated software and hardware for their numerical analysis—are quite common in science, engineering, and applied technology.
Earth–Moon Satellite Orbits
For an example of a program whose efficient solution requires adaptive step size methods, we consider an Apollo satellite in orbit about the Earth E and Moon M. Figure7.7.8 shows an x1x2x1x2-coordinate system whose origin lies at the center of mass of the Earth and the Moon and which rotates at the rate of one revolution per “moon month”of approximately τ=27.32τ=27.32 days, so the Earth and Moon remain fixed in their positions on the x1x1-axis. If we take as unit distance the distance (about 384,000 kilometers, assumed constant) between the Earth and Moon centers, then their coordinates are E(−μ,0)E(−μ,0) and M(1−μ,0),M(1−μ,0), where μ=mM/(mE+mM)μ=mM/(mE+mM) in terms of the Earth mass mEmE and Moon mass mMmM. If we take the total mass mE+mMmE+mM as the unit of mass and τ/(2π)≈4.35τ/(2π)≈4.35 days as the unit of time, then the gravitational constant is G=1G=1 in Eq.(26), and the equations of motion of the satellite position S(x1,x2)S(x1,x2) are
where rErE and rMrM denote the satellite’s distance to the Earth and Moon (indicated in Fig.7.7.8). The initial two terms on the right-hand side of each equation result from the rotation of the coordinate system. In the system of units described here, the lunar mass is approximately mM=0:012277471.mM=0:012277471. The second-order system in (27) can be converted to an equivalent first-order system (of four differential equations) by substituting
x′1=x3,x′2=x4,so thatx″1=x′3,x″2=x′4.
x'1=x3,x'2=x4,so thatx′′1=x'3,x′′2=x'4.
Suppose that the satellite initially is in a clockwise circular orbit of radius about 2400 kilometers about the Moon. At its farthest point from the Earth (x1=0.994)(x1=0.994) it is “launched” into Earth–Moon orbit with initial velocity v0.v0. The corresponding initial conditions are
x1(0)=0.994,x2(0)=0,x3(0)=0,x4(0)=−v0.
x1(0)=0.994,x2(0)=0,x3(0)=0,x4(0)=−v0.
An adaptive step size method (ode45) in the Matlab software system was used to solve numerically the system in (27). The orbits in Figs.7.7.9 and 7.7.10 were obtained with
v0=2.031732629557andv0=2.001585106379,
v0=2.031732629557andv0=2.001585106379,
respectively. [In the system of units used here, the unit of velocity is approximately 3680 km/h.] In each case a closed but multilooped periodic trajectory about the Earth and the Moon—a so-called bus orbit—is obtained, but a relatively small change in the initial velocity changes the number of loops! For more information, see NASA Contractor Report CR-61139, “Study of the Methods for the Numerical Solution of Ordinary Differential Equations,” prepared by O. B. Francis, Jr. et al. for the NASA–George C. Marshall Space Flight Center, June 7, 1966.
So-called Moon–Earth “bus orbits” are periodic—that is, are closed trajectories traversed repeatedly by the satellite—only in a rotating x1x2x1x2-coordinate system as discussed above. The satellite of Fig.7.7.9 traverses its closed orbit and returns to rendezvous with the Moon about 48.4 days after its insertion into orbit. Figures7.7.11 and 7.7.12 illustrate the motion of the same satellite—but in an ordinary nonrotating xy-coordinate system centered at the Earth, in which the Moon encircles the Earth counterclockwise in a near-circular orbit, completing one revolution in about 27.3 days. The Moon starts at point S, and after 48.4 days it has completed a bit over 1.75 revolutions about the Earth and reaches the point R at which its rendezvous with the satellite occurs. Figure7.7.11 shows the positions of Moon and satellite a day and a half after the satellite’s insertion into its orbit, each traveling in a generally counterclockwise direction around the Earth. Figure7.7.12 shows their positions a day and a half before their rendezvous at point R, the satellite meanwhile having encircled the Earth about 2.5 times in an orbit that (in the indicated xy-coordinate system) appears to resemble a slowly varying ellipse.
7.7 Problems
A hand-held calculator will suffice for Problems 1 through 8. In each problem an initial value problem and its exact solution are given. Approximate the values ofx(0.2) andy(0.2) in three ways: (a) by the Euler method with two steps of sizeh=0.1h=0.1; (b) by the improved Euler method with a single step of sizeh=0.2h=0.2; and (c) by the Runge–Kutta method with a single step of sizeh=0.2.h=0.2.Compare the approximate values with the actual valuesx(0.2) andy(0.2).
A computer will be required for the remaining problems in this section. In Problems 9 through 12, an initial value problem and its exact solution are given. In each of these four problems, use the Runge–Kutta method with step sizesh=0.1h=0.1andh=0.05h=0.05to approximate to five decimal places the valuesx(1) andy(1) Compare the approximations with the actual values.
Crossbow bolt Suppose that a crossbow bolt is shot straight upward with initial velocity 288 ft/s. If its deceleration due to air resistance is (0.04)v, then its height x(t) satisfies the initial value problem
x″=−32−(0.04)x′;x(0)=0,x′(0)=288.
x′′=−32−(0.04)x';x(0)=0,x'(0)=288.
Find the maximum height that the bolt attains and the time required for it to reach this height.
Repeat Problem13, but assume instead that the deceleration of the bolt due to air resistance is (0.0002)v2.(0.0002)v2.
Height of projectile Suppose that a projectile is fired straight upward with initial velocity v0 from the surface of the earth. If air resistance is not a factor, then its height x(t) at time t satisfies the initial value problem
d2xdt2=−gR2(x+R)2;x(0)=0,x′(0)=v0.
d2xdt2=−gR2(x+R)2;x(0)=0,x'(0)=v0.
Use the values g=32.15 ft/s2≈0.006089 mi/s2 for the gravitational acceleration of the earth at its surface and R=3960 mi as the radius of the earth. If v0=1 mi/s, find the maximum height attained by the projectile and its time of ascent to this height.
Batted Baseball
Problems 16 through 18 deal with the batted baseball of Example4, having initial velocity 160 ft/s and air resistance coefficientc=0.0025.
Find the range—the horizontal distance the ball travels before it hits the ground—and its total time of flight with initial inclination angles 40∘,45∘, and 50∘.
Find (to the nearest degree) the initial inclination that maximizes the range. If there were no air resistance it would be exactly 45∘, but your answer should be less than 45∘.
Find (to the nearest half degree) the initial inclination angle greater than 45∘ for which the range is 300 ft.
Home run Home run Find the initial velocity of a baseball hit by Babe Ruth (with c=0.0025 and initial inclination 40∘) if it hit the bleachers at a point 50 ft high and 500 horizontal feet from home plate.
Crossbow bolt Crossbow bolt Consider the crossbow bolt of Problem14, fired with the same initial velocity of 288 ft/s and with the air resistance deceleration (0.0002)v2 directed opposite its direction of motion. Suppose that this bolt is fired from ground level at an initial angle of 45∘. Find how high vertically and how far horizontally it goes, and how long it remains in the air.
Artillery projectile Artillery projectile Suppose that an artillery projectile is fired from ground level with initial velocity 3000 ft/s and initial inclination angle 40∘. Assume that its air resistance deceleration is (0.0001)v2. (a) What is the range of the projectile and what is its total time of flight? What is its speed at impact with the ground? (b) What is the maximum altitude of the projectile, and when is that altitude attained? (c) You will find that the projectile is still losing speed at the apex of its trajectory. What is the minimum speed that it attains during its descent?