Open Educational Resources

Finding Roots of Equations: Problems

  1. Let f:[-5,5]\rightarrow \mathbb{R} be such that f(x)=0.5+x\cos{x}. Use the graphical method to find the number and estimates of the roots of each of the equations f(x)=0 and f(x)=2.
  2. Compare the bisection method and the false position method in finding the roots of the two equations f(x)=0 and f(x)=0.1 in the interval [-1,1] where f(x) is:

        \[f(x)=2x^3-4x^2+3x\]


    Compare with the result using the Solve function and comment on the use of \varepsilon_r as a measure of the relative error.
  3. Locate the first nontrivial root of \sin{x}=x^2 where x is in radians. Use the graphical method and the bisection method with \varepsilon_s=0.002.
  4. Water is flowing in a trapezoidal channel at a rate of Q= 20 m^3/s. The critical depth y for such a channel must satisfy the equation:

        \[ 1-\frac{Q^2}{g A_c^3}B=0\]


    where g = 9.81 m/s^2, A_c is the cross sectional area (m^2), and B is the width of the channel at the surface (m). Assuming that the width and the cross‐sectional area can be related to the depth y by:

        \[\begin{split}B&=3+y\\A_c&=3y+\frac{y^2}{2}\end{split}\]


    Solve for the critical depth using the graphical method, the bisection method, and the false position method to find the critical depth in the interval [0,3]. Use \varepsilon_s=0.001. Discuss your results (difference in the methods, accuracy, number of iterations, etc.).

    (Note: The critical depth is the depth below which the flow of the fluid becomes relatively fast and affected by the upstream conditions (supercritical flow). Above that depth, the flow is subcritical, relatively slow and is controlled by the downstream conditions)

  5. Use the graphical method to find all real roots of the equation

        \[f(x)=-14-20x+19x^2-3x^3\]


    You are required to

     

    • Plot the function several times by using smaller ranges for the x axis to obtain the roots up to 3 significant digits.
    • Explain the steps you have taken to identify the number and estimates of the roots.
    • Explain the steps you have taken to obtain the roots to the required accuracy.
    • Compare your results to the solution obtained using the “Solve” function in Mathematica.
  6. Use the bisection method and the false position method to find all real roots of the equation

        \[f(x)=-14-20x+19x^2-3x^3\]


    Use \varepsilon_s=0.001.
  7. Consider the function:

        \[f(x)= -0.9 x^2+1.7x+2.5\]


    Determine a root for the expression f(x)=0

     

    • using fixed-point iteration
    • using the Newton-Raphson method

    Use an initial guess of x_0 = 5 and \varepsilon_s=0.01\%=0.0001. Compare your final answer with the Solve function in Mathematica.

  8. Consider the function:

        \[f(x)=7\sin{(x)} e^{-x}-1\]


    Find the lowest positive root for the expression f(x)=0

     

    • graphically
    • using the Newton-Raphson method
    • using the secant method
  9. Use the Newton-Raphson method to find the root of:

        \[f(x)=e^{-0.5x}(4-x)-2=0\]


    Employ initial guesses of 2, 6, and 8. Explain your results.
  10. The Manning equation can be written for a rectangular open channel as:

        \[Q=\frac{\sqrt{S}(BH)^\frac{5}{3}}{n(B+2H)^{\frac{2}{3}}}\]


    where Q is the fluid flow in m^3/s, S is the slope in m/m, H is the depth in m, and n is the Manning roughness coefficient. Develop a fixed-point iteration scheme to solve this equation for H given Q = 5, S =0.0002, B = 20, and n = 0.03. Consider \varepsilon_s=0.0005. Prove that your scheme converges for all initial guesses greater than or equal to zero.
  11. Why does the Babylonian method almost always converge?
  12. Consider the function

        \[f(x)=e^{-x}-x\]


    Compare the secant method and the Newton-Raphson method in finding the root of the equation f(x)=0.
  13. Solve the previous question using the Fixed-Point Iteration method
  14. Use the Fixed-Point Iteration method to find the roots of the following functions:

        \[f_1(x)=1-\frac{x^2}{4} \quad f_2(x)=2\sqrt{x-1}-x\]

  15. Engineers of many disciplines use thermodynamics extensively in their work. The following polynomial can be used to relate the zero-pressure specific heat of dry air c_p kJ/(kg K) to temperature (K):

        \[c_p = 0.99403 + 1.671\times10^{-4}T + 9.7215\times10^{-8}T^2 - 9.5828\times 10^{-11}T^3 + 1.9520\times10^{-14}T^4\]


    Plot c_p versus T in the range T\in [0, 1200]\ K​. Determine the temperature that corresponds to a specific heat of 1.1 kJ/(kg K) using the following methods. In each case use an approximate relative error tolerance of 0.0001% and report the number of iterations required to achieve that tolerance.
    – Bisection Method. Set up the initial brackets [0,1200] K.
    – Fixed-Point Iteration. Using an initial guess of T = 500 K.
    – Newton-Raphson Method. Using an initial guess of T = 500 K.
    – Secant Method. Use initial guesses of T = 540 K and T = 550K.
    Make a conclusion about the efficiency of the methods.

Leave a Reply

Your email address will not be published. Required fields are marked *