Skip to main content

what are numerical methods and fixed points

 



Unveiling the Magic of Successive Substitution: The Simplest Numerical Method for Root Finding 🌱



Are you tired of scratching your head over complex numerical methods for finding roots? Fear not, for in the vast world of mathematics, simplicity often reigns supreme. Enter the humble yet powerful technique known as successive substitution, also referred to as the fixed-point method. 🎩✨



For example, if 𝑓 is defined on the real numbers by:



𝑓(𝑥) = 𝑥² - 3𝑥 + 4,



then 2 is a fixed point of 𝑓, because 𝑓(2) = 2. 🔒✨


### Unraveling the Mystery


Successive substitution is a numerical method used to find the roots of a function by iteratively applying a transformation until convergence is achieved. It's like solving a puzzle one piece at a time, gradually getting closer to the solution with each step.


### How It Works


Let's dive into a simple Python example to illustrate the magic of successive substitution:


```python

def fixed_point_iteration(g, x0, tol=1e-6, max_iter=100):

    """Fixed-point iteration to find the root of g(x)=x."""

    x = x0

    for _ in range(max_iter):

        x_new = g_xnew(x)

        if abs(x_new - x) < tol:

            return x_new

        x = x_new

    return None


# Define your function g(x)

#p=rt/(v-b)-a/v2  

#p(v-b)=rt-a(v-b)/v2

#pv2=v2rt/(v-b)-a

#isolate one v =f(v) for eithwr format                                

def g_xnew(x):

    return 


# Initial guess

x0 = 1.5


# Apply fixed-point iteration

root = fixed_point_iteration(g, x0)

print("Approximate root:", root)

```


In this example, we want to find the v using the function. We start with an initial guess of \( x_0 = 1.5 \) and iteratively update \( x \) until convergence.


### The Beauty of Simplicity


Successive substitution shines in its simplicity. With just a few lines of code, you can embark on a journey to uncover the roots of any function. Its ease of implementation makes it an excellent choice for quick prototyping and solving problems on the fly.


### Caveats and Considerations


While successive substitution is straightforward, it's essential to choose a suitable initial guess and ensure that the function \( g(x) \) satisfies the conditions for convergence. Additionally, the method may converge slowly for certain functions or fail to converge altogether if approached haphazardly.


Successive substitution may fail to converge when the derivative of the function being iterated does not satisfy certain conditions, particularly around the fixed point.


One such condition for convergence is that the derivative of the function at the fixed point must be less than 1 in absolute value. This condition ensures that successive iterations move closer to the fixed point rather than diverging away from it. When the absolute value of the derivative at the fixed point is greater than 1, the iterations may oscillate or diverge instead of converging.


This phenomenon can be understood intuitively: if the slope of the function at the fixed point is greater than 1, the iterations will amplify rather than dampen, leading to divergence. Conversely, if the slope is less than 1, the iterations will gradually approach the fixed point, resulting in convergence.


In summary, successive substitution may not converge when the derivative of the function at the fixed point exceeds 1 in absolute value, leading to divergence instead of convergence. This highlights the importance of analyzing the derivative of the function to ensure convergence in numerical methods.


### Conclusion


In a world where complexity often clouds our vision, the simplicity of successive substitution offers a beacon of hope for numerical problem solvers. With its straightforward approach and elegant solutions, it proves that sometimes, less truly is more. So why not give it a try and unlock the hidden roots of your mathematical challenges? 🌟


Let's embrace the beauty of simplicity and embark on a journey of numerical exploration with successive substitution at our side. Happy coding! 💻🔍


Interactive example:

https://colab.research.google.com/drive/1Sg1K-LAsSwznV1YzfL7JUaG8Df3t5hs_?usp=sharing

Comments

Popular posts from this blog

yer a wizard - making your own custom %%cellmagics for colab notebooks

    ## Dive into the Magic of Jupyter %%cellmagic! ✨📚 Hey there, fellow data enthusiasts! 👋 Today, let's dive into the fascinating world of Jupyter's `%%cellmagic` 🪄. This little-known feature can supercharge your Jupyter Notebook workflow! 🚀💡 ### What's `%%cellmagic`? 🧙‍♂️✨ In the Jupyter ecosystem, `%magic` and `%%cellmagic` commands add special functionalities to your notebook cells. Think of them as magical commands that can transform how your cells behave! 🌟 For example, `%%time` can measure the execution time of a cell. But what if you want to create your own custom magic? That's where `%%cellmagic` shines! 💥 ### Example: Create Your Own Cell Magic 🛠️🔮 Let's say you want to create a custom magic that processes a cell in a specific way. Here's a simple example to get you started: ```python from IPython.core.magic import (Magics, magics_class, cell_magic, needs_local_scope) @magics_class class class_mycellmagic(Magics):     @needs_local_scope ...

x=? or how can i make a random variable in python ?

 **Unleashing the Power of Randomness in Python/Numpy for Simple Game Structures! 🎲🔀🃏** Welcome, fellow programmers, game enthusiasts, and curious minds! Today, we embark on an exciting journey into the realm of randomness within Python and Numpy. Whether you're a seasoned coder or a newbie explorer, buckle up as we uncover the magic of random functions and how they can breathe life into simple game structures. 🚀 **1. Uniform Randomness:** 🎲 Ah, the beauty of unpredictability! With Python's `random` module or Numpy's `numpy.random` package, we can effortlessly generate uniformly distributed random numbers. This feature is ideal for scenarios like rolling dice, selecting random players, or determining the movement of objects in a game world. ```python import random # Roll a fair six-sided die roll_result = random.randint(1, 6) print("You rolled:", roll_result) ``` **2. List Choice:** 🔀 In the realm of games, sometimes decisions need to be made from a pool of ...

how to do the linear regression in python??

  📊 **Unlocking the Power of Linear Regression with Python's SciPy Library!** 📈 Hey there, data enthusiasts! Today, we're diving into the world of linear regression using Python's powerful SciPy library. Strap in as we explore how to perform linear regression, calculate the coefficient of determination (R-squared), and unleash the potential of your data with just a few lines of code! ### 📊 What is Linear Regression? Linear regression is a fundamental statistical technique used to model the relationship between two variables. It's like fitting a straight line to a scatter plot of data points, allowing us to make predictions and understand the underlying relationship between the variables. ### 💻 Let's Get Coding! First things first, fire up your Python environment and make sure you have SciPy installed. If not, a quick `pip install scipy` should do the trick. Once that's done, import the necessary libraries: ```python from scipy.stats import linregress ``` Now...