Using Python write a 4th-order Runge-Kutta method to numerically solve the following forcedamped harmonic oscillator equationx = −2(w₁x - w²2x + Fmxcwith initial conditions x(0)= 1 and r(0)= 0. In addition, use ( = 1/10, wo = 10and w = wo for plot 1 and w = 2wn for plot 2. Note, you will be solving the equation twice,one for a driving force equal to sin wot and the second with a driving force equal to sin 2wot.h = 0.1 and 0 ≤t≤40Check for step size convergence and once achieved, plot both solutions.

Detailed explanation: Solving this problem using a 4th-order Runge-Kutta method in Python, we need…

Answered: Using Python write a 4th-order…

C++ for Engineers and Scientists

4th Edition

ISBN:9781133187844

Author:Bronson, Gary J.

Publisher:Bronson, Gary J.

Chapter5: Repetition Statements

Section5.7: Do While Loops

Problem 6E: (Numerical analysis) Here’s a challenging problem for those who know a little calculus. The...

See similar textbooks

Similar questions

(Numerical analysis) Here’s a challenging problem for those who know a little calculus. The Newton-Raphson method can be used to find the roots of any equation y(x)=0. In this method, the (i+1)stapproximation,xi+1,toarootofy(x)=0 is given in terms of the ith approximation, xi, by the following formula, where y’ denotes the derivative of y(x) with respect to x: xi+1=xiy(xi)/y(xi) For example, if y(x)=3x2+2x2,theny(x)=6x+2 , and the roots are found by making a reasonable guess for a first approximation x1 and iterating by using this equation: xi+1=xi(3xi2+2xi2)/(6xi+2) a. Using the Newton-Raphson method, find the two roots of the equation 3x2+2x2=0. (Hint: There’s one positive root and one negative root.) b. Extend the program written for Exercise 6a so that it finds the roots of any function y(x)=0, when the function for y(x) and the derivative of y(x) are placed in the code.
(Mechanics) The deflection at any point along the centerline of a cantilevered beam, such as the one used for a balcony (see Figure 5.15), when a load is distributed evenly along the beam is given by this formula: d=wx224EI(x2+6l24lx) d is the deflection at location x (ft). xisthedistancefromthesecuredend( ft).wistheweightplacedattheendofthebeam( lbs/ft).listhebeamlength( ft). Eisthemodulesofelasticity( lbs/f t 2 ).Iisthesecondmomentofinertia( f t 4 ). For the beam shown in Figure 5.15, the second moment of inertia is determined as follows: l=bh312 b is the beam’s base. h is the beam’s height. Using these formulas, write, compile, and run a C++ program that determines and displays a table of the deflection for a cantilevered pine beam at half-foot increments along its length, using the following data: w=200lbs/ftl=3ftE=187.2106lb/ft2b=.2fth=.3ft
(Conversion) An object’s polar moment of inertia, J, represents its resistance to twisting. For a cylinder, this moment of inertia is given by this formula: J=mr2/2+m( l 2 +3r 2 )/12misthecylindersmass( kg).listhecylinderslength(m).risthecylindersradius(m). Using this formula, determine the units for the cylinder’s polar moment of inertia.
(Thermodynamics) The work, W, performed by a single piston in an engine can be determined by this formula: W=Fd F is the force provided by the piston in Newtons. d is the distance the piston moves in meters. a. Determine the units of W by calculating the units resulting from the right side of the formula. Check that your answer corresponds to the units for work listed in Table 1.1. b. Determine the work performed by a piston that provides a force of 1000 N over a distance of 15 centimeters.
(Statics) An annulus is a cylindrical rod with a hollow center, as shown in Figure 6.7. Its second moment of inertia is given by this formula: I4(r24r14) I is the second moment of inertia (m4). r2 is the outer radius (m). r1 is the inner radius (m). a. Using this formula, write a function called annulusMoment ( ) that accepts two double-precision numbers as parameters (one for the outer radius and one for the inner radius), calculates the corresponding second moment of inertia, and displays the result. b. Include the function written in Exercise 5a in a working program. Make sure your function is called from main(). Test the function by passing various data to it.
(Practice) Determine names for functions that do the following: a. Find the average of a set of numbers. b. Find the area of a rectangle. c. Find the minimum value in a set of numbers. d. Find the density of a steel door. e. Sort a set of numbers from lowest to highest
(Civil eng.) Modify the program written for Exercise 9 to determine the maximum load that can be placed at the end of an 8-foot I-beam, shown in Figure 2.21, so that the stress on the fixed end is 20,000lbs/in2. Use the fact that this beam’s rectangular moment of inertia is 21.4 in4 and the value of c is 3 in.
(Electrical eng.) You’ve been asked to write a C++ program to calculate the total resistance of a series circuit. In this circuit, the total resistance is the sum of all individual resistance values. The circuit consists of a number of 56-ohm, 33-ohm, and 15-ohm resistors. a. For this programming problem, how many outputs are required? b. How many inputs does this problem have? c. Determine a formula for converting input items into output items. The number of 56-ohm resistors is m, the number of 33-ohm resistors is n, and the number of 15-ohm resistors is p. d. Test the formula written for Exercise 1c using the following sample data: m=17,n=24,andp=12.
(Program) Write a program that tests the effectiveness of the rand() library function. Start by initializing 10 counters, such as zerocount, onecount, twocount, and so forth, to 0. Then generate a large number of pseudorandom integers between 0 and 9. Each time 0 occurs, increment zerocount; when 1 occurs, increment onecount; and so on. Finally, display the number of 0s, 1s, 2s, and so on that occurred and the percentage of time they occurred.
Find the error(s) in the following code: (6)
Determine the value of each of the following expressions. (For decimal numbers, round your answer to two decimal places.) (1) abs(-18) fabs(20.5) fabs(-87.2) pow(4, 2.0) pow(8.4, 3.5) sqrt(7.84) sqrt (196.0) sqrt (38.44)* pow(2.4, 2) / fabs(-3.2) floor(27.37) ceil(19.2) floor(12.45) + ceil(6.7) floor (-8.9) + ceil (3.45) floor(9.6) / ceil(3.7) pow(-4.0, 6.0) pow(10, -2.0) pow(9.2, 1.0 / 2)
(Practice) For the following correct algebraic expressions and corresponding incorrect C++ expressions, find the errors and write corrected C++ expressions: Algebra C++ Expression a.(2)(3)+(4)(5)(2)(3)+(4)(5) b. 6+1826+18/2 c. 4.512.23.1S4.5/12.23.1 d. 4.6(3.0+14.9)4.6(3.0+14.9) e. (12.1+18.9)(15.33.8)(12.1+18.9)(15.33.8)

SEE MORE QUESTIONS

Recommended textbooks for you

C++ for Engineers and Scientists

Computer Science

ISBN:

9781133187844

Author:

Bronson, Gary J.

Publisher:

Course Technology Ptr

Operations Research : Applications and Algorithms

Computer Science

ISBN:

9780534380588

Author:

Wayne L. Winston

Publisher:

Brooks Cole

C++ Programming: From Problem Analysis to Program…

Computer Science

ISBN:

9781337102087

Author:

D. S. Malik

Publisher:

Cengage Learning

C++ for Engineers and Scientists

Computer Science

ISBN:

9781133187844

Author:

Bronson, Gary J.

Publisher:

Course Technology Ptr

Operations Research : Applications and Algorithms

Computer Science

ISBN:

9780534380588

Author:

Wayne L. Winston

Publisher:

Brooks Cole

C++ Programming: From Problem Analysis to Program…

Computer Science

ISBN:

9781337102087

Author:

D. S. Malik

Publisher:

Cengage Learning

Systems Architecture

Computer Science

ISBN:

9781305080195

Author:

Stephen D. Burd

Publisher:

Cengage Learning

Programming Logic & Design Comprehensive

Computer Science

ISBN:

9781337669405

Author:

FARRELL

Publisher:

Cengage

EBK JAVA PROGRAMMING

Computer Science

ISBN:

9781337671385

Author:

FARRELL

Publisher:

CENGAGE LEARNING - CONSIGNMENT

SEE MORE TEXTBOOKS