Numerical Methods for Engineers
Numerical Methods for Engineers
7th Edition
ISBN: 9780073397924
Author: Steven C. Chapra Dr., Raymond P. Canale
Publisher: McGraw-Hill Education
bartleby

Concept explainers

Differential Equations
Have you ever observed the movement of the satellites and rockets or how the planets go around the sun in their orbits? How will you determine their motion? Definitely they follow some scientific laws and these kinds of problems are framed as differentia…
bartleby

Videos

Textbook Question
Book Icon
Chapter 24, Problem 44P

A jet fighter's position on an aircraft carrier's runway was timed during landing:

t, s 0 0.52 1.04 1.75 2.37 3.25 3.83
x, m 153 185 208 249 261 271 273

where x is the distance from the end of the carrier. Estimate (a) velocity ( d x / d t ) and (b) acceleration ( d v / d t ) using numerical differentiation.

(a)

Expert Solution
Check Mark
To determine

To calculate: The velocity of the jet fighter landing on an aircraft carrier’s runaway using Numerical Integration.

Answer to Problem 44P

Solution:

The velocities at different points of time of a Jet fighter are given below,

t,s 0 0.52 1.04 1.75 2.37 3.25 3.83
x,m 153 185 210 249 261 271 273
v,m/s 68.26923 54.80769 50.97398 35.93855 16.05180 6.59273 0.303817

Explanation of Solution

Given Information:

The position of a Jet fighter on an aircraft carrier’s runway is,

t,s 0 0.52 1.04 1.75 2.37 3.25 3.83
x,m 153 185 210 249 261 271 273

Formula Used:

Second-order Lagrange interpolating polynomial.

f′(t)=f(ti−1)2t−ti−ti+1(ti−1−ti)(ti−1−ti+1)+f(ti)2t−ti−1−ti+1(ti−ti−1)(ti−ti+1)+f(ti+1)2t−ti−1−ti(ti+1−ti−1)(ti+1−ti)…… (1)

Here, t is the value at which the derivative is to be estimated.

Calculation:

Consider the function f(x) to be position (x) as shown in table above. So, f′(x) will represent the acceleration (v).

Apply second-order Lagrange interpolating polynomial to each set of three adjacent points for calculation of velocities at different points.

Calculate f′(0).

Substitute t=0 and i=1 in equation(1).

f′(0)=f(0)2(0)−t1−t2(t0−t1)(t0−t2)+f(t1)2(0)−t0−t2(t1−t0)(t1−t2)+f(t2)2(0)−t0−t1(t2−t0)(t2−t1)

Substitute the value of function from above table.

f′(0)=153(2(0)−0.52−1.04(0−0.52)(0−1.04))+185(2(0)−0−1.04(0.52−0)(0.52−1.04))+210(2(0)−0−0.52(1.04−0)(1.04−0.52))=68.26923

Calculate f′(0.52).

Substitute t=0.52 and i=2 in equation(1).

f′(0.52)=f(0.52)2(0.52)−t2−t3(t1−t2)(t1−t3)+f(t2)2(0.52)−t1−t3(t2−t1)(t2−t3)+f(t3)2(0.52)−t1−t2(t3−t1)(t3−t2)

Substitute the value of function from above table.

f′(0.52)=153(2(0.52)−0.52−1.04(0−0.52)(0−1.04))+185(2(0.52)−0−1.04(0.52−0)(0.52−1.04))+210(2(0.52)−0−0.52(1.04−0)(1.04−0.52))=54.80769

Calculate f′(1.04).

Substitute t=1.04 and i=3 in equation(1).

f′(1.04)=f(1.04)2(1.04)−t3−t4(t2−t3)(t2−t4)+f(t3)2(1.04)−t2−t4(t3−t2)(t3−t4)+f(t4)2(1.04)−t2−t3(t4−t2)(t4−t3)

Substitute the value of function from above table.

f′(1.04)=185(2(1.04)−1.04−1.75(0.52−1.04)(0.52−1.75))+210(2(1.04)−0.52−1.75(1.04−0.52)(1.04−1.75))+249(2(1.04)−0.52−1.04(1.75−0.52)(1.75−1.04))=50.97398

Calculate f′(1.75).

Substitute t=1.75 and i=1 in equation(1).

f′(1.75)=f(1.75)2(1.75)−t4−t5(t3−t4)(t3−t5)+f(t4)2(1.75)−t3−t5(t4−t3)(t4−t5)+f(t5)2(1.75)−t4−t3(t5−t3)(t5−t4)

Substitute the value of function from above table.

f′(1.75)=210(2(1.75)−1.75−2.37(1.04−1.75)(1.04−2.37))+249(2(1.75)−1.04−2.37(1.75−1.04)(1.75−2.37))+261(2(1.75)−1.04−1.75(2.37−1.04)(2.37−1.75))=35.93855

Calculate f′(2.37).

Substitute t=2.37 and i=5 in equation(1).

f′(2.37)=f(2.37)2(2.37)−t5−t6(t4−t5)(t4−t6)+f(t5)2(2.37)−t4−t6(t5−t6)(t5−t4)+f(t6)2(2.37)−t4−t5(t6−t4)(t6−t5)

Substitute the value of function from above table.

f′(2.37)=249(2(2.37)−2.37−3.25(1.75−2.37)(1.75−3.25))+261(2(2.37)−1.75−3.25(2.37−1.75)(2.37−3.25))+271(2(2.37)−1.75−2.37(3.25−1.75)(3.25−2.37))=16.05180

Calculate f′(3.25).

Substitute t=3.25 and i=6 in equation(1).

f′(3.25)=f(3.25)2(3.25)−t6−t7(t5−t6)(t5−t7)+f(t6)2(3.25)−t6−t7(t6−t5)(t6−t7)+f(t7)2(3.25)−t5−t6(t7−t6)(t7−t5)

Substitute the value of function from above table.

f′(3.25)=261(2(3.25)−3.25−3.83(2.37−3.25)(2.37−3.83))+271(2(3.25)−2.37−3.83(3.25−2.37)(3.25−3.83))+273(2(3.25)−2.37−3.25(3.83−2.37)(3.83−3.25))=6.59273

Calculate f′(3.83).

Substitute t=3.83 and i=7 in equation(1).

f′(3.83)=f(3.83)2(3.83)−t7−t8(t6−t7)(t6−t8)+f(t7)2(3.83)−t6−t8(t7−t6)(t7−t8)+f(t8)2(3.83)−t6−t8(t8−t6)(t8−t7)

Substitute the value of function from above table.

f′(3.83)=261(2(3.83)−3.25−3.83(2.37−3.25)(2.37−3.83))+271(2(3.83)−2.37−3.83(3.25−2.37)(3.25−3.83))+273(2(3.83)−2.37−3.25(3.83−2.37)(3.83−3.25))=0.303817

The velocities at different points of time of a Jet fighter are given below,

t,s 0 0.52 1.04 1.75 2.37 3.25 3.83
x,m 153 185 210 249 261 271 273
v,m/s 68.26923 54.80769 50.97398 35.93855 16.05180 6.59273 0.303817

(b)

Expert Solution
Check Mark
To determine

To calculate: The acceleration of the jet fighter landing on an aircraft carrier’s runaway using Numerical Integration.

Answer to Problem 44P

Solution:

The complete table of acceleration for different values of velocity is given below,

t,s x,m v=dxdt a=dvdt
0 153 68.26923 -35.14510
0.52 185 54.80769 -16.63004
1.04 210 50.97398 -13.20841
1.75 249 35.93855 -26.99478
2.37 261 16.05180 -23.26046
3.25 271 6.59273 -10.80560
3.83 273 0.303817 -10.88029

Explanation of Solution

Given Information:

The velocities at different points of time of a Jet fighter are given below as calculated in Part (a).

t,s 0 0.52 1.04 1.75 2.37 3.25 3.83
x,m 153 185 210 249 261 271 273
v,m/s 68.26923 54.80769 50.97398 35.93855 16.05180 6.59273 0.303817

Formula Used:

Second-order Lagrange interpolating polynomial.

f′(t)=f(ti−1)2t−ti−ti+1(ti−1−ti)(ti−1−ti+1)+f(ti)2t−ti−1−ti+1(ti−ti−1)(ti−ti+1)+f(ti+1)2t−ti−1−ti(ti+1−ti−1)(ti+1−ti)…… (1)

Here, t is the value at which the derivative is to be estimated.

Calculation:

Consider the function f(x) to be velocity (v) as shown in table above. So, f′(x) will represent the acceleration (a).

Calculate the acceleration (a) at different points with respect to the velocities (v) given in the above table with the help of Equation (1).

Calculate f′(0).

Substitute t=0 and i=1 in equation(1).

f′(0)=f(0)2(0)−t1−t2(t0−t1)(t0−t2)+f(t1)2(0)−t0−t2(t1−t0)(t1−t2)+f(t2)2(0)−t0−t1(t2−t0)(t2−t1)

Substitute the value of function from above table.

f′(0)=68.26923(2(0)−0.52−1.04(0−0.52)(0−1.04))+54.80769(2(0)−0−1.04(0.52−0)(0.52−1.04))+50.97398(2(0)−0−0.52(1.04−0)(1.04−0.52))=−35.14510

Calculate f′(0.52).

Substitute t=0.52 and i=2 in equation(1).

f′(0.52)=f(0.52)2(0.52)−t2−t3(t1−t2)(t1−t3)+f(t2)2(0.52)−t1−t3(t2−t1)(t2−t3)+f(t2)2(0.52)−t1−t2(t3−t1)(t3−t2)

Substitute the value of function from above table.

f′(0.52)=68.26923(2(0.52)−0.52−1.04(0−0.52)(0−1.04))+54.80769(2(0.52)−0−1.04(0.52−0)(0.52−1.04))+50.97398(2(0.52)−0−0.52(1.04−0)(1.04−0.52))=−16.63004

Calculate f′(1.04).

Substitute t=1.04 and i=3 in equation(1).

f′(1.04)=f(1.04)2(1.04)−t3−t4(t2−t3)(t2−t4)+f(t1)2(1.04)−t2−t4(t3−t2)(t3−t4)+f(t4)2(1.04)−t2−t3(t4−t2)(t4−t3)

Substitute the value of function from above table.

f′(1.04)=54.80769(2(1.04)−1.04−1.75(0.52−1.04)(0.52−1.75))+50.97398(2(1.04)−0.52−1.75(1.04−0.52)(1.04−1.75))+35.93855(2(1.04)−0.52−1.04(1.75−0.52)(1.75−1.04))=−13.20841

Calculate f′(1.75).

Substitute t=1.75 and i=4 in equation(1).

f′(1.75)=f(1.75)2(1.75)−t4−t5(t3−t4)(t3−t5)+f(t1)2(1.75)−t3−t5(t4−t3)(t4−t5)+f(t2)2(1.75)−t4−t3(t5−t3)(t5−t4)

Substitute the value of function from above table.

f′(1.75)=50.97398(2(1.75)−1.75−2.37(1.04−1.75)(1.04−2.37))+35.93855(2(1.75)−1.04−2.37(1.75−1.04)(1.75−2.37))+16.05180(2(1.75)−1.04−1.75(2.37−1.04)(2.37−1.75))=−26.99478

Calculate f′(2.37).

Substitute t=2.37 and i=1 in equation(1).

f′(2.37)=f(2.37)2(2.37)−t1−t2(t0−t1)(t0−t2)+f(t1)2(2.37)−t0−t2(t1−t0)(t1−t2)+f(t2)2(2.37)−t0−t1(t2−t0)(t2−t1)

Substitute the value of function from above table.

f′(2.37)=35.93855(2(2.37)−2.37−3.25(1.75−2.37)(1.75−3.25))+16.05180(2(2.37)−1.75−3.25(2.37−1.75)(2.37−3.25))+6.59273(2(2.37)−1.75−2.37(3.25−1.75)(3.25−2.37))=−23.26046

Calculate f′(3.25).

Substitute t=3.25 and i=1 in equation(1).

f′(3.25)=f(3.25)2(3.25)−t6−t7(t5−t6)(t5−t7)+f(t6)2(3.25)−t6−t7(t6−t5)(t6−t7)+f(t7)2(3.25)−t5−t6(t7−t6)(t7−t5)

Substitute the value of function from above table.

f′(3.25)=16.05180(2(3.25)−3.25−3.83(2.37−3.25)(2.37−3.83))+6.59273(2(3.25)−2.37−3.83(3.25−2.37)(3.25−3.83))+0.303817(2(3.25)−2.37−3.25(3.83−2.37)(3.83−3.25))=−10.80560

Calculate f′(3.83).

Substitute t=3.83 and i=7 in equation(1).

f′(3.83)=f(3.83)2(3.83)−t7−t8(t6−t7)(t6−t8)+f(t7)2(3.83)−t6−t8(t7−t6)(t7−t8)+f(t8)2(3.83)−t6−t8(t8−t6)(t8−t7)

Substitute the value of function from above table.

f′(3.83)=16.05180(2(3.83)−3.25−3.83(2.37−3.25)(2.37−3.83))+6.59273(2(3.83)−2.37−3.83(3.25−2.37)(3.25−3.83))+0.303817(2(3.83)−2.37−3.25(3.83−2.37)(3.83−3.25))=−10.88029

Want to see more full solutions like this?

Subscribe now to access step-by-step solutions to millions of textbook problems written by subject matter experts!
Previous
Chapter 24, Problem 43P
Next
Chapter 24, Problem 45P
Students have asked these similar questions
Q2) A: Find the region where ODEs has no limit cycle: x = y + x³ y=x+y+y³ 6
Q3)A: Given H(x,y)=x2-x+ y²as a first integral of an ODEs, find this ODES corresponding to H(x,y) and show the phase portrait by using Hartman theorem and by drawing graph of H(x,y)-e. Discuss the stability of critical points of the corresponding ODEs.
Q/ Write Example is First integral but not Conservation system.

Chapter 24 Solutions

Numerical Methods for Engineers

Ch. 24 - Perform the same computation as Sec. 24.1, but...Ch. 24 - 24.2 Repeat Prob. 24.1, but use Romberg...Ch. 24 - 24.3 Repeat Prob. 24.1, but use a two- and a...Ch. 24 - 24.4 Integration provides a means to compute how...Ch. 24 - Use numerical integration to compute how much mass...Ch. 24 - 24.6 Fick’s first diffusion law states...Ch. 24 - The following data were collected when a large oil...Ch. 24 - 24.8 You are interested in measuring the fluid...Ch. 24 - Prob. 10PCh. 24 - 24.11 Glaucoma is the second leading cause of...
Ch. 24 - One of your colleagues has designed a new...Ch. 24 - Video an giography is used to measure blood flow...Ch. 24 - 24.14 Perform the same computation as in Sec....Ch. 24 - Perform the same computation as in Sec. 24.2, but...Ch. 24 - 24.16 As in Sec. 24.2, compute F using the...Ch. 24 - Stream cross-sectional areas (A) are required for...Ch. 24 - 24.18 As described in Prob. 24.17, the...Ch. 24 - 24.21 A transportation engineering study requires...Ch. 24 - 24.22 A wind force distributed against the side of...Ch. 24 - 24.23 Water exerts pressure on the upstream ...Ch. 24 - 24.24 To estimate the size of a new dam, you have...Ch. 24 - The data listed in the following table gives...Ch. 24 - The heat flux q is the quantity of heat flowing...Ch. 24 - 24.27 The horizontal surface area of a lake at a...Ch. 24 - 24.28 Perform the same computation as in Sec....Ch. 24 - 24.29 Repeat Prob. 24.28, but use five...Ch. 24 - Repeat Prob. 24.28, but use Romberg integration to...Ch. 24 - Faradays law characterizes the voltage drop across...Ch. 24 - 24.32 Based on Faraday’s law (Prob. 24.31), use...Ch. 24 - Suppose that the current through a resistor is...Ch. 24 - If a capacitor initially holds no charge, the...Ch. 24 - 24.35 Perform the same computation as in Sec....Ch. 24 - 24.36 Repeat Prob. 24.35, but use (a) Simpson’s ...Ch. 24 - 24.37 Compute work as described in Sec. 24.4, but...Ch. 24 - As was done in Sec. 24.4, determine the work...Ch. 24 - 24.39 The work done on an object is equal to the...Ch. 24 - The rate of cooling of a body (Fig. P24.40) can be...Ch. 24 - 24.41 A rod subject to an axial load (Fig....Ch. 24 - If the velocity distribution of a fluid flowing...Ch. 24 - 24.43 Using the following data, calculate the work...Ch. 24 - 24.44 A jet fighter’s position on an aircraft...Ch. 24 - 24.45 Employ the multiple-application Simpson’s...Ch. 24 - The upward velocity of a rocket can be computed by...Ch. 24 - Referring to the data from Problem 20.61, find the...Ch. 24 - Fully developed flow moving through a 40-cm...Ch. 24 - Fully developed flow of a Bingham plasticfluid...Ch. 24 - 24.50 The enthalpy of a real gas is a ...Ch. 24 - Given the data below, find the isothermal work...Ch. 24 - 24.52 The Rosin-Rammler-Bennet (RRB) equation is...Ch. 24 - For fluid flow over a surface, the heat flux to...Ch. 24 - The pressure gradient for laminar flow through a...Ch. 24 - 24.55 Velocity data for air are collected at...
Knowledge Booster
Background pattern image
Advanced Math
Learn more about
Need a deep-dive on the concept behind this application? Look no further. Learn more about this topic, advanced-math and related others by exploring similar questions and additional content below.
Similar questions
SEE MORE QUESTIONS
Recommended textbooks for you
Text book image
Algebra: Structure And Method, Book 1
Algebra
ISBN:9780395977224
Author:Richard G. Brown, Mary P. Dolciani, Robert H. Sorgenfrey, William L. Cole
Publisher:McDougal Littell
Text book image
Trigonometry (MindTap Course List)
Trigonometry
ISBN:9781337278461
Author:Ron Larson
Publisher:Cengage Learning
Text book image
Algebra & Trigonometry with Analytic Geometry
Algebra
ISBN:9781133382119
Author:Swokowski
Publisher:Cengage
Text book image
College Algebra
Algebra
ISBN:9781337282291
Author:Ron Larson
Publisher:Cengage Learning
SEE MORE TEXTBOOKS
01 - What Is A Differential Equation in Calculus? Learn to Solve Ordinary Differential Equations.; Author: Math and Science;https://www.youtube.com/watch?v=K80YEHQpx9g;License: Standard YouTube License, CC-BY
Higher Order Differential Equation with constant coefficient (GATE) (Part 1) l GATE 2018; Author: GATE Lectures by Dishank;https://www.youtube.com/watch?v=ODxP7BbqAjA;License: Standard YouTube License, CC-BY
Solution of Differential Equations and Initial Value Problems; Author: Jefril Amboy;https://www.youtube.com/watch?v=Q68sk7XS-dc;License: Standard YouTube License, CC-BY