/*! This file is auto-generated */ .wp-block-button__link{color:#fff;background-color:#32373c;border-radius:9999px;box-shadow:none;text-decoration:none;padding:calc(.667em + 2px) calc(1.333em + 2px);font-size:1.125em}.wp-block-file__button{background:#32373c;color:#fff;text-decoration:none} Q15E Use the fourth-order Runge–Kut... [FREE SOLUTION] | 91Ó°ÊÓ

91Ó°ÊÓ

Use the fourth-order Runge–Kutta subroutine with h= 0.1 to approximate the solution to\({\bf{y' = cos}}\;{\bf{5y - x,y(0) = 0}}\),at the points x= 0, 0.1, 0.2, . . ., 3.0. Use your answers to make a rough sketch of the solution on\(\left[ {{\bf{0,3}}} \right]\).

Short Answer

Expert verified

\({{\rm{x}}_{\rm{n}}}\)

\({{\rm{y}}_{\rm{n}}}\)

0.5

0.21462

1.0

0.13890

1.5

-0.02668

2.0

-0.81879

2.5

-1.69491

3.0

-2.99510

The rough sketch is given below.

Step by step solution

01

Find the values of \({{\bf{k}}_{\bf{i}}}{\bf{.i = 1,2,3,4}}\)

Using the improved 4th order Runge-Kutta subroutine.

Since \({\bf{f(x,y) = cos}}\;{\bf{5y - x}}\) and \({\bf{x = }}{{\bf{x}}_{\bf{0}}}{\bf{ = 0,y = }}{{\bf{y}}_{\bf{o}}}{\bf{ = 0}}\) and h = 0.1, M = 30

\(\begin{array}{l}{{\bf{k}}_{\bf{1}}}{\bf{ = h}}{\rm{f}}{\bf{(x,y) = 0}}{\bf{.1(cos5y - x)}}\\{{\bf{k}}_{\bf{2}}}{\bf{ = hf}}\left( {{\bf{x + }}\frac{{\bf{h}}}{{\bf{2}}}{\bf{,y + }}\frac{{{{\bf{k}}_{\bf{1}}}}}{{\bf{2}}}} \right){\bf{ = 0}}{\bf{.1}}\left( {{\bf{cos}}\left( {{\bf{5y + }}\frac{{{{\bf{k}}_{\bf{1}}}}}{{\bf{2}}}} \right){\bf{ - (x + 0}}{\bf{.05)}}} \right)\\{{\bf{k}}_{\bf{3}}}{\bf{ = hf}}\left( {{\bf{x + }}\frac{{\bf{h}}}{{\bf{2}}}{\bf{,y + }}\frac{{{{\bf{k}}_{\bf{2}}}}}{{\bf{2}}}} \right){\bf{ = 0}}{\bf{.1}}\left( {{\bf{cos}}\left( {{\bf{5y + }}\frac{{{{\bf{k}}_{\bf{2}}}}}{{\bf{2}}}} \right){\bf{ - (x + 0}}{\bf{.05)}}} \right)\\{{\bf{k}}_{\bf{4}}}{\bf{ = hf}}\left( {{\bf{x + h,y + }}{{\bf{k}}_{\bf{3}}}} \right){\bf{ = 0}}{\bf{.1(y + }}{{\bf{k}}_{\bf{3}}}{\bf{)(cos(5(y + }}{{\bf{k}}_{\bf{3}}}{\bf{) - (x + 0}}{\bf{.05))}}\\{{\bf{k}}_{\bf{1}}}{\bf{ = h(x,y) = 0}}{\bf{.1}}\\{{\bf{k}}_{\bf{2}}}{\bf{ = hf}}\left( {{\bf{x + }}\frac{{\bf{h}}}{{\bf{2}}}{\bf{,y + }}\frac{{{{\bf{k}}_{\bf{1}}}}}{{\bf{2}}}} \right){\bf{ = 0}}{\bf{.091891}}\\{{\bf{k}}_{\bf{3}}}{\bf{ = hf}}\left( {{\bf{x + }}\frac{{\bf{h}}}{{\bf{2}}}{\bf{,y + }}\frac{{{{\bf{k}}_{\bf{2}}}}}{{\bf{2}}}} \right){\bf{ = 0}}{\bf{.092373}}\\{{\bf{k}}_{\bf{4}}}{\bf{ = hf}}\left( {{\bf{x + h,y + }}{{\bf{k}}_{\bf{3}}}} \right){\bf{ = 0}}{\bf{.079522}}\end{array}\)

02

Find the values of x and y

\(\begin{array}{c}{\bf{x = 0 + 0}}{\bf{.1 = 0}}{\bf{.1}}\\{\bf{y = 0 + }}\frac{{\bf{1}}}{{\bf{6}}}\left( {{{\bf{k}}_{\bf{1}}}{\bf{ + 2}}{{\bf{k}}_{\bf{2}}}{\bf{ + 2}}{{\bf{k}}_{\bf{3}}}{\bf{ + }}{{\bf{k}}_{\bf{4}}}} \right)\\{\bf{ = 0}}{\bf{.091342}}\end{array}\)

The solution of the given IVP at x=0.1 is approx. 0.091.

03

Evaluate the other values of x and y

x

y

x

y

x

y

0.2

0.157

1.2

0.087

2.2

-1.173

0.3

0.195

1.3

0.055

2.3

-1.300

0.4

0.212

1.4

0.019

2.4

-1.454

0.5

0.215

1.5

-0.027

2.5

-1.695

0.6

0.208

1.6

-0.086

2.6

-2.037

0.7

0.196

1.7

-0.170

2.7

-2.309

0.8

0.180

1.8

-0.306

2.8

-2.501

0.9

0.161

1.9

-0.535

2.9

-2.698

1

0.139

2

-0.819

3

-2.995

1.1

0.114

2.1

-0.029

04

Plot the Graph

Hence the solution is

\({{\rm{x}}_{\rm{n}}}\)

\({{\rm{y}}_{\rm{n}}}\)

0.5

0.21462

1.0

0.13890

1.5

-0.02668

2.0

-0.81879

2.5

-1.69491

3.0

-2.99510

Unlock Step-by-Step Solutions & Ace Your Exams!

  • Full Textbook Solutions

    Get detailed explanations and key concepts

  • Unlimited Al creation

    Al flashcards, explanations, exams and more...

  • Ads-free access

    To over 500 millions flashcards

  • Money-back guarantee

    We refund you if you fail your exam.

Over 30 million students worldwide already upgrade their learning with 91Ó°ÊÓ!

One App. One Place for Learning.

All the tools & learning materials you need for study success - in one app.

Get started for free

Most popular questions from this chapter

Use the fourth-order Runge–Kutta subroutine with h = 0.25 to approximate the solution to the initial value problemy'=1-y,y(0)=0, at x = 1. Compare this approximation with the one obtained in Problem 6 using the Taylor method of order 4.

By experimenting with the improved Euler’s method subroutine, find the maximum value over the interval [0,2]of the solution to the initial value problemy'=sin(x+y),y(0)=2 Where does this maximum value occur? Give answers to two decimal places.

Sailboats A and B each have a mass of 60 kg and cross the starting line at the same time on the first leg of a race. Each has an initial velocity of 2 m/sec. The wind applies a constant force of 650 N to each boat, and the force due to water resistance is proportional to the velocity of the boat. For sailboat A the proportionality constants arebefore planing when the velocity is less than 5 m/sec andwhen the velocity is above 5 m/sec. For sailboat B the proportionality constants arebefore planing when the velocity is less than 6 m/sec andwhen the velocity is above. If the first leg of the race is 500 m long, which sailboat will be leading at the end of the first leg?

Show that when the trapezoid scheme given in formula (8) is used to approximate the solutionf(x)=exofy'=y,y(0)=1 , at x = 1, then we get yn+1=1+h21-h2yn,n = 0, 1, 2, . . . , which leads to the approximation (1+h21-h2)1hfor the constant e.Compute this approximation for h= 1,10-1,10-2,10-3,and10-4and compare your results with those in Tables 3.4 and 3.5.

Escape Velocity.According to Newton’s law of gravitation, the attractive force between two objects varies inversely as the square of the distances between them. That is,Fg=GM1M2/r2whereM1andM2are the masses of the objects, ris the distance between them (center to center), Fgis the attractive force, and Gis the constant of proportionality. Consider a projectile of

constant mass mbeing fired vertically from Earth (see Figure 3.12). Let trepresent time and v the velocity of the projectile.

  1. Show that the motion of the projectile, under Earth’s gravitational force, is governed by the equationdvdt=-gR2r2, where ris the distance between the projectile and the center of Earth, Ris the radius of Earth, Mis the mass of Earth, andg=GM/R2.
  2. Use the fact that dr /dt = v to obtainvdvdt=-gR2r2
  3. If the projectile leaves Earth’s surface with velocity vo, show thatv2=2gR2r+vo2-2gR
  4. Use the result of part (c) to show that the velocity of the projectile remains positive if and only ifvo2-2gR>0. The velocityve=2gRis called the escape velocityof Earth.
  5. If g= 9.81 m/sec 2and R= 6370 km for Earth, what is Earth’s escape velocity?
  6. If the acceleration due to gravity for the moon is gm= g/6 and the radius of the moon is Rm= 1738 km, what is the escape velocity of the moon?
See all solutions

Recommended explanations on Math Textbooks

View all explanations

What do you think about this solution?

We value your feedback to improve our textbook solutions.

Study anywhere. Anytime. Across all devices.